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PREFACE 


This volume is the first of four which make up the Multisatellite Attitude 
Detcrminatioii/Optical Aspect Bias Determination (MSAD/OABIA^ System 
Description and Operating Guide. The volumes are 

Volume 1 - Introduction and Analysis 
Volume 2 - System Description 
Volume 3 - Operating Guide 

Volume 4 - Program Listing and Sample Execution 

This volume contains an introductory exposition of the MSAD/OABIAS System 
and describes the analytic basis for the OABIAS subsystem. This includes a 
detailed discussion of the recursive estimator algorithm, each of flie 12 state 
vector elements, and the 8 observation models used. 

Volume 2 de^ribes the system flow and the components of the MSAD/OABIAS 
System. The table language description in this volume provides detailed infor- 
mation relating the operational displays on the IBM 2250 display device to 
specific COMMON areas and subroutines within the MSAD/OABIAS System. 

Volume 3 contains a con^lete description of all MSAD/OABIAS NAMEDCST 
control parameters, a description ar't sample of all piinted output unique to 
OABIAS and of each IBM 2250 graphics display, an explanation of and user 
response for all error messages generated by the MSAD/OABIAS Sj'stem, 
and a listing of Uie Job Control Language (JCL) required to operate the system. 

Volume 4 contains the program listing with supplementary output and line 
printer plots of all IBM 2250 displays occurring during a sample execution 
of die program. Hiis volume preserves, in source form, the MSAD/OABIAS 
system as it is presented in tliis document. 
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This document dcscritjcs the MuUisatcllite Attitude Dctermination/C^tical 
Aspect Bias Deterroioation (AISAD/OABIA^ System, designed to determine 
spin axis orientation and biases in the alignment or performance of optical or 
infrared horizon sensors and Sun sensors used for spacecraft attitude deter- 
mination. MSAD/OABIAS uses any combination of eight observation models to 
process data from a single onboard horizon sensor and Sun sensor to determine 
simultaneously the two components of the attitude of the spacecraft, the initial 
plmse of the Sun sensor, the spin rate, seven sensor biases, and the orbital 
in-track error associated mth the spacecraft ephemeris information supplied 
to the system, hi addition, the MSAD/OABIAS System pro\ddes a data simu- 
lator for system and performance testing, an independent deterministic attitude 
system for preprocessing and independent testily of biases determined, and a 
multipurpose data prediction and comparison system. 

MSAD/OABIAS has extensive capabilities for an interactive graphics mode and 
makes use of the Graphics Executive Support System (GESS^, formerly known 
as the Multisatellite Attitude Determination System (MSAD) services. MSAD/ 
OABIAS is a multisatellite system capable of supporting, in its present form, 

3 

the Small Scientific Satellite (S ), the Interplanetary Monitoring Platform (MP), 
the Atmosphere Explorer (AE), and the Synchronous Meteorological Satellite 
(SMS) missions or an^' similar missions using optical or infrared horizon 
scanners and providing attitude data that can be read by the MSAD/OABIAS 
System. 
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SECTION 1 - INTRO^yUCTION 
1.1 MOTI\^ATTON FOR BIAS DETERMINATION 

In the simplcsl, spin-stobilixed satellite attitude determination problem, the 
state vector consi;>ts of only t^vo clemcn<^s A>'.ich define the orientation of the 
angular momentum vector in inertial space. In practice it is known that the 
presence of biases caused by the misalignment of sensors, or by misspecified 
hardware characteristics, \vill cause -‘naccuracies or the complete breakdown 
of the simple attitude determination computations. This can lead to inaccu-* 
racies in attitude control if the spacecraft spin axis orientation cannot be deter- 
mined to Mthin control system uncertainties, or worse, it can lead to attitude 
determination blackouts during which computations yield no solution for attitude. 
Therefore, in more realistic analyses, additional elements, e. g. , sensor 
mounting angles which affect the accuracy of an attitude estimate, are included 
in the state vector computation to provide improved knowledge of the values 
of parameters. 

Accurate bias dcterminahc»a will permit more accurate spacecraft attitude 
determination and control liian could be accomplished purely on the basis of 
prclaunc't measurements. Highly accurate prelaunch measurements not only 
are difficult and expensive to obtain but also may be invalidated by postlaunch 
changes in the spacecraft due to th< rmal or mechanical shocks. Also, space- 
craft d 3 mamic imbalance effects which may appear as sensor biases are ex- 
pensive to remove completely before launch and may chaiige due to disci*ete 
events in the mission profile, such as apogee motor firing or boom aeployment. 

The Multisatellito Attitude Dctcrniinalioa/Optical Aspect Bias Determination 
(MSaD/OAIjIAS) System in Us present form provides the necessary bias de- 
termination and spin axis attitude refinement for spin-stabilized spacecraft 
nith sensor confirurations which include Sun sensors and horizon sensors. 
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The sensor observables are accurately mcdclcd and are used in a recursive 
least- squares filter technique to estimate the elements of a state vector, which 
include the tw’O comijonents of the attitude; Sun sensor azimuth, elevation and 
plane tilt; Earth sensor azimuth and elevation; Eax*th sensor-triggering thresh- 
old; and fixed time delays on sensor triggerings. In addition, a state vector 
component is used to correct the in- track orbital position, '..hich is the most 
likely orbital error as well as tlic largest single orbital uncei*tainty. 

The interactive graphics capabilities programmed into MSAD/OAB*^AS provide 
for a synergistic collaboration between analyst and computer. Generally, it is 
not feasible to solve for spin axis attitude, sensor misalignments, and orbital 
position simultaneously. Analyst intervention is required tx> select the data 
span and to choose the variables to be solved using that data span. The com- 
puter program can best evaluate the numerical parameters and the fit to the 
data. Analyst scrutiny is required further to detect the possibility of unmodcled 
effects in the data, which may arise from the use of an early orbit estimate or 
from anomalous sensor performance. 

MSAD/OABIAS is, therefore, an interactive data processing system for de- 
termining the maximum information content of a spacecraft attitude sensor 
data set. 

1.2 HISTORICAL BACKGROUND 

Prior to the development of the MSAD/OABIAS System, the Optical Aspect 
Attitude Determination System (OASYS) u'as used to evaluate biases. During 
mission support, biases were obsci*ved through large dispersions bct'> ecn 
single-frame deterministic attitude solutions or through large spans ol data 
which yielded no solution. Several trial and cri'or techniques were dcv'elopcd 
to adjust bias parameter on input to OASYS to reduce the attitude solution dis- 
persion and to improve the quality of the fit to the data. This latter adjustment 
was accomplished by displaying the raw data superimposed on a computer 
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pi'cdiction of the data based on the system-determined attitudes nod system- 
estimated biases. In addition, it was possible to determine fixed-angle biases 
on the Sun angle, the nadir angle, and the dihedral angles calculatea for input 
to the differential correction subsystem of OASYS. Biases which could not be 
treated as fLxcd-angle biases due to the Earth sensor-triggering tiu'cshold (such 
as, the sensor-mounting angle bias or the bias oh the apparent radius -of the 
Earth) could be determined only by the manual trial and error techniques. Al- 
though it was possible to develop considerable analyst expertise in recogpuzing 
What biases could be applied to fit the data, these techniques had several draw- 
backs. 

One drawback was that because the nadir angles or dihedral angles are computed 
through nonlinear transformation of measurements, the biases of these angles 
v.'ere not necessarily physically meaningful. Also, if more dian one bias was 
present, manual trial and error techniques were too slow for real-time sui^rt 
and led to unce<.'tainties about the uniqueness of the bias solution determined. 

MSAD/OABIAS xma designed to address the above difflculties. It uses as input 
the true measiu'cments received in the telemetry, and its interactive graphics 
control structure is designed to speed operational and analytical processing. 
MSAD/OABIAS evolved from a study of analytical techniques for recursive 
least-squares Kalman filters (Reference 1). The preliminary design (Refer- 
ence 2) added a bias determination subsystem to the existing OASYS at the point 
at wluch the differential correction subsystem xvas invoked. 

MSAD/OABL\S is a multisntcllitc system and has been used on data from tlic 

3 

Interplanetary Monitoxang Platform (BMP), the Small Scientific Satellite (S ), 
the Radio Astronomy Explox*cr-B (RAE-B), the Atmosphere Explorer-C (AE-C) 
and the Synchronous Meteorological Satellite-A and -B (SMS-A, -B). The user 
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simply supplies MSAD/OAlilAS with a sensor-measurement data set in a 
standard format.^ 

1.3 SYSTEM OVERVIEW 

All interactive gx'aphics display sendees and the interactive control of program 
flow in MSAD/OABIAS are accomplished through the Graphics Executive Support 
System (GESS). This system may be considered an overall executive control 
program or a sidxsystem of programs which supply graphics support to MSAO/ 
OABIAS at specified points in the program flow. In addition to GESS, five main 
functional subs^'stems exist in MSAD/OABIAS. These are 

1. Data selection and adjustment subsystem (OADRIV) 

2. Data simulator subsystem (ODAP) 

3. Deterministic altitixle determination processor sub^stem (OASY^ 

4. Recursive least-squares filter subsystem (OABIAI^ 

5. Data prediction and solution comparison subsj'stem (PLOTOQ 

Also, an auxiliary' routine (AECOPY) copies AE-C spacecraft data in a form 
readable by the MSAD/OABIAS System. Figure 1-1 shows these major func- 
tional subs\*stems in a progpram baseline diagpram hierarchy. OAMAIN and 
OPMAIN are drivers and are shown for the sake of accuracy. . 

In the standard processing sequence, the GESS executive calls flie system 
driver OADRIV. OADRIV either obtains the data to be processed from an 
OABIAS measurement data set or causes the data simulator ODAP to be in- 
voked. In the case of AE-C data, AECOPY must be invoked by OADRIV to 
produce the standard OABIAS measurement data set. OADRIV then calls the 
dolci'ministic processor OASYS (through OPMAIN). and the deterministic so- 
lution is used as an a priori estimate in the rccursiw least-squares Alter 


1 -■ 

Due to the exigencies of the AE-C mission, this sUindard format w;;s violated, 

necessitating the development of a utility interface routine. 
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Figure 1-1. MSAD/OABIAS Functional Subsystems 
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subsystem OAHIAS. The user may call the data prediction and graphical solu- 
tion comparison subs^tem PLOTOC fi*om OABIAS or asynehronoi.sly from 
OASYS. (Refer to Volume III for a thorough presentation of the mechanics of 
the graphical processing options.') 

MSAD/OAGIAS has multiple (^tions for accessii^ ephemcris data which arc 
required by all of the subsystems. Cations arc available through the multi- 
satellite routine EPIIKMG to read spacecraft ephemcrides as well as lunar and 
solar ephemcrides. The position of the spacecraft can be obtained from cither 
the Definitive Orbit Determination System (DOD^ file (Reference 3) or the 
Goddard Trajectory Determination System (GTDS) file (Reference 4). It also 
can be internally generated through a simple orbit generator (ORBGEN). The 
positions of the Sun and the Moon can be obtained from a Jet Propulsion Lab- 
oratory (JPL) ephemcris file (Reference 5) or a SUNRD file (Reference 6). 
These positions likewise can be generated internally by using SUNl for the 
Son and ORBGEN for the Moon, hitemally-genr rated ephemcrides are useful 
for simulation purposes. The user is not constrained to use the same ephem- 
eris source h\ simulating and processing data. 

1.3.1 OADRn" 

OADRIV is the main control subroutine for the MSAD/OABIAS System. Opera- 
tionally, it serves as the data evaluation, selection, and adjustment subsystem. 
Although the rest c-f the MSAD/OABIAS Sj'stem is limited to 200 frames of data, 
if core storai,c is not limited, OADRIV can accept over 1200 frames. This 
enabh rapid and detailed data evaluation and selection by the operator. The 
OADRIV capacity is limited to approximately 1200 frames because the storage 
capacity of the IBM 2250 display device buffer is limited. 

OADRIV provides a variety of options for sifting the data, including periodic 
reading of the data by cither time or frame number. Also, individual data 
po.ntj> in either plot or character displays may be flagged, noise or biases 
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may be added to the data (for c\*aluation or testing purposes) , and the appro- 
priate interval may be scleetcd for processing by the attitude and bias deter- 
mination subsystems. 

The capability of the user to edit the data entered in th"* determination systems 
is crucial. Because of the \*ariablc quality of attitude data and the subtle data 
rejection decisions which must be made based on Imowledge of the hardware 
or an understanding of spurious events (such as boom reflections), operator 
intervention in the execution of the program is required. Rejection of spurious 
data could be done by machine, but imdoubtedly, this would require an unde- 
terminable amount of additional storage. 

1.3.2 ODAP - 

ODAP, the data simulator subsystem in MSAD/OABIAS, is a modified graphics 
version of the program described in Reference 7. It allows all system capabil- 
ities to be exercised in the simulation mode, which is essential to any system 
as large and complex as MSAD/OABIAS. An attached simulator not only facil- 
itates system testing but also makes it feasible. The simulator is also used 
for prelaunch analysis and simulations. Data can be simulated for any mission 
conditions, and noise and biases can be applied and passed through to the atti- 
tude processors. This simulation-processing sequence can be repeated as 
rtotscssary during a single execution of the system. The data simulator is not 
normally used during mission support and can be overlayed when not in use. 

1.3.3 OASYS 

The OASYS subsystem, an interactive graphics version of the program de- 
scribed in Rcfci'cncc 8, serves as a deterministic attitude determination 
processor within MSAD/OABIAS. Each data frame consists of the Sun angle, 
the Sun crossing time, and the Earth-in and Earth-out triggering times, and 
can be used to compute attitude by four geometric mcihods — Earth-in, Earth- 
out, Earth width, and mid- scan dihedral angle. The average of these single 
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frame solutions provides an initial estimate to the reeursivc least-squares filter 
in OAniAS. 

In addition graphical plots of the framc-bj'-frame altitude solutions obtained 
through the different methods give insight into what biases ma 3 '^ be present. 

The Earth width method, for example, is sensitive to a bias on the apparent 
angular radius of the Earth and, in the event of sueh a bias, yields results 
whieh are different from the mid-sean dihedral angle method whieh is not sen- 
sitive to this bins. As a check on tlie biases determined, the biases can be 
entered in OASYS and the solutions observed graphically'. The correct solu- 
tion will have the effect of reducing the attitude dispersion in the single-frame 
solutions and minimizing the difference bettveen the curves for the different 
solution methods. 

1.3.4 OABIAS 

OABIAS is the subsy'stem which determines the biases and from which the sys- 
tem as a whole takes its name. It is a recursive least-squares filter designed 
to determine spin axis attitude, sensor biases, and in-track orbit error (Refer- 
ence 2). The 12-component state vector is composed of the following elements: 

• X inertial component of the spin axis^ ^ 

• Y inertial component of the spin axis^ 

• Spin rate 

• Phase of I'otation measured from the projection of the sunline in 
the spacecraft spin plane 

• Seven sensor bias parameters 

• In-track orbiUil error 

^This component is converted to right ascension and declination of the spin 
axis in all output displays. 
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The OAIJLNS subsystem uses eight models (functions) of the obscn'ablcs which 
are dependent on the state vector elements in the filtering scheme. The pro- 
gram provides graphical displays of the state vector elements on a point-by- 
point basis as well as displays of statistical parameters; c.g. , residuals and 
correlation coefficients. The user has control of all input parameters from 
the graphics device and can make multiple passes through the data to achieve 
an optimal solution. (See Section 1.4 for details of the filter and the observa- 
tion models.) 

1.3.5 PLOTOC 

The PLOTOC subsystem provides plots of the predicted and observed Earth 
widths and rotation angles. PLOTOC will display up to three predicted plo^s 
simultaneously. Tipically this data will be displayed with the deterministic 
solution input to OABIAS and the OABIAS solution. The degree of fit to the 
data can then be determined visually. 

Parameters can be varied manually to see the effect on the observed fit. 
Although these parameters — which include height of the atmosphere, oblate- 
ness coefficients, and orbital elements — are not in the state vector, their 
potential effect on the solution must be gauged. The PLOTOC plots clearly 
show unmodeled effects. PLOTOC provides a check on the validity of the 
OABIAS solutions and, because all GESS-generated plots can be routed to 
hardcopy CalComp plots, a way to document the results. 

1.4 OABIAS RECURSIVE LEAST-SQUARES FILTER 

1.4.1 Introduction 

The OABIAS I'ccursivc least-squares filler is based on eight distinct obser- 
V'alion models. They are 

• Model 1 — Sun angle model ^ 

• Model 2 — .Sun sighting time model 
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• Model 3 — Nadir vector projection model 

• Model 4 — Horizc i crossing time model 

. • Model 5 — Sun to Earth-in and Sun to Earth-out dihedral angle model 

- • \ 

• Model G — Earth width model 
0 Model 7 — Small target model 

0 Model 8 — Sun to Earth mid-se; \ dihedral angle model 

(See Section 1.4.2 for individual descriptions.) Any combination of these 
models may be used. However, statistically independent models should be 
chosen to obtain valid statistical uncertainties. For example, if the obser- 
vations available arc the Sun angle, the Sun sighting time, and the times at 
which the horizon sensor acquires and loses the disk of the Earth, a valid 
choice of statisticallj' independent models would be Mod 's' 6, and 8. The 
addition of any other models, such as Model 5, to process the same observa- 
tions would yield unrealistically low uncertainties. 

Each obsen'ation model is formulated as a function of a subset of the 12 state 
vector components which include two components of the spin axis attitude, the 
spin rate, the initial phase of the X axis (Sun sensor), and the following eight 
biases (seven sensor biases and one orbit parameter): 

% 

1. Horizon sensor mounting angle bias (equivalent to a bias in the 
elevation of the sensor relative to the spin plane) 

2. Azimuth bias of the horizon sensor relative to the Sun sensor on 
horizon-in triggering 

3. Azimuth bias of the horizon sensor relative to the Sun sensor on 
horizon-out triggering 

4. Bias on the ajigular radius of the central body (equivalent to a 
sensor-triggering level bias) 

i* 

5. Sun angle bias 
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C. Sun sensor plane tilt 

7. Panoramic attitude sensor plane tilt (used for RAE-13 mission) 

8. Orbital in-track error 

It is not possible to determine all of the state vector elements from an}' one 
observation model. The dependence of the observation models on the state 
vector components is summarized in Table 1-1, 

Within OADIAS, data is processed one frame at a time. lii order to facilitate 
analysis^ options are available to update the state vector after 

• . Each observation is processed 

• Each frame is processed 

o A specified number of frames is processed 

the option is taken to update the state vector on a frame-by-frame basis, an 
acMitional option is available to iter^;te through all data frames since the last 
update. In thi. ivay, any desired degree of recursiveness is proWded so that 
even a batch processing differential correction can be obtained. 

1.4.2 O>scrvation Models 

Of the eight OABIAS observation models, two are associated with the Sun and 
six with the central body (eith'^r the Earth or the Moon). The formulation of 
these observables is summarized here. (See Section 3 of this volume for ad- 
ditional details.) 

• Model 1— Sun angle model. The observable for Model 1 is the Sun 
angle reported with the data. The predicted value is determined 
by the position of the Sun (fi*om either a solar cphcmcris or an 
algebraic routine) relative to the spacecraft, the spin axis attitude, 
the Sun sensor plane tilt, and ttic Sun angle bias. 
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Model 2 — Sun sighting lime model. The obscrvntion for Model 2 is 
the projection of the Sun vector onto n vector normal to the plane 
containing the spin axis and Sun sensor slit. The expression for 
this projection is evaluated at the Sun sighting time. Since the Sun 
must lie in the plane containing the spin axis and the Sun sensor 
slit when it is sighted, the observed value is al\v.ays zero. 

Model 3 — Nadir veetor projection model. The observation for 
Model 3 is the projection of a unit vector along the optical axis of 
the horizon sensor onto the nadir vector. The observation is eval- 
uated at the time of a horizon crossing. 

Model 4 — Horizon crossing time model. The observation for this 
model is the projection of the horizon vector onto the normal to the 
plane containing the spin axis and the optical axis of the horizon 
sensor. It is evaluated at the time of a horizon crossing. There- 
fore, as in the case of Model 2, the observed value is zero. 

Model 5 — Sun to Earth-in and Sun to Earth-out dihedral angle model. 
The two observables for Model 5 are the rotation angles from the 
Sun to Earth-in crossing and from the Sun to Earth-out crossing. 
These observables are available directly from the data. 

Model G — Earth width model. The observable for Model 6 is the 
Earth width dihedral angle, defined as the Earth-out rotation angle 
minus the Earth-in rotation angle. This model docs not depend on 
a Sun sighting. Therefore, it may be used when the Sun is not vis- 
ible or when a Sun sighting reference time is not available from the 
data (as with the AE-(!J wheel-mounted hoiizon sensors). 

Model 7 — Small target model. This model is used when the ecnlral 
body is treated as a point source; c. g. , the RAE-n mission. The 
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observable is the sensor mounting angle; i. e, , the angle between 
the spm axis of the spacecraft and the nadir vector. The nadir vec- 
tor is assumed to be along the line of sight of the horizon sensor. 

• Model 8 — Sun to Earth mid-scan dihedral angle model. The ob- 

servable for this model is the dihcdi'al angle from the Sun to the 
mid- point between the Earth-in crossing and the Earth-* rossing. 
For a spherical Earth this would be the same as the dih- angle 
from the Sun to the center of the Earth. However, for an oblate 
Earth the perpendicular bisector of the great circle between the 
Earth-in crossing and the Earth-out crossing generally does not 
go through the Earth's center. 

Models 6 and 8 are statistically independent whereas the two parts of Models 4 
and 5 are not. Therefore, Models 1, 6, and 8 generally would be prefer ed 
for processing over Models 1 and 5 or Models 1 and 4. However, Models 0 
and 8 depend on both horizon crossings being available, whereas Models 4 and 5 
treat the t^vo ho*'izon crossings separately. Thus, Modv*)! 4 or 5 would be used 
when only a single horizon crossing was available; e. g. , a sensor operating in 
the visible range and triggering on one horizon crossing and the terminator. 
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SECTION 2 - SENSOU DESCIUPTIONS 

This section describes the Sun sonsors and the body-fixed horizon oc lectors 
which usually arc used onboard satellites supported by the MSAD/OARIAS Sys- 
- tem. Brief descriptions of the wheel- mounted horizon scanner used by the 
AE-C satcliilc and the panoramic attitude scanner (PAS) used on the RAE-B 
satellite also are included. MSAD/OABIA3 supported both the AE-C and the 
RAE-B missions. The RAE-B PAS application necessitated OABIAS modifi- 
cations which arc described throughout this document. 

The geometry of the Sun-sensing and horizon-sensing operations is of primary 
significance to the MSAO/OABIAS System and, therefore, is emphasized. De- 
tailed physical and internal characteristics of currently available sensors are 
less important to MSAD/OABIAS and, hence, are not described in this docu- 
ment. 

2 . 1 SUN SENSOR 

Figure 2-1 shows the main geometric features of t'le Sun sensor system which 
has been assumed in the design of MSAD/OABIAS. The sensor's ffeld-of-view 
(FOV) is fan-shaped anu is indicated in Figure 2-1 by the heavy line. The FOV 
is centered about the face axis which is at an angle 4 ;^ from the satellite 
equatorial plane. A common value for the total FOV angle t is * 28 degrees. 
The sensor is mounted such that its nominal FOV plane contains .he satellite's 

Ai /V 

spin .axis S. (In Figure 2-1, S is showai lying in the FOV plane.) Sun sensor 

systems commonly employ two sensors moumed with their face axis x on 

S8 

opposite sides of the satellite's r*-uatorial plane. This arrangement permits 
full coverage of tlic celestial sphere with Iwx) 128-dcgrce sensors and provides 
rcdundaii'’'- over the salcllilc's equatorial region. 

The lot:- .1 of the satellite causes the Sun sensor's FOV to sweep out a wide 
belt on the celestial sphere. In Figure 2-1, this Is the area between the two 
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hca\y ckishcd lines. The sensor sees the Sun once per spin period if the Sun 
lies within tais belt. Olhcrudsc, it does not see the Sun at all. The Sun sighting 
points oecur when the rotating FOV plane crosses the Sun vector U . 


Each time tlic Siui is sighted, the Sun sensor S 3 ’Stcm provides MSAD/OADIAS 

with two pieces of information — the time t at wiiich the sighting occurred and 

s 

the elevation angle f of the Sun abo\'C the sensor face axis x at t . The 
° ^s ss s 

oidmard sensor itself normally docs not establish t^ . It provides a reference 

pulse at each sighting which is used elseuiiere in the spacecraft. Normally, 

t is a time tag which is placed on the telemetry data on the gpround. The sat- 
s 

ellite spin rate m is computed from the time intervals between successive 

Sun sighting pulses. The elevation angle ^ is used to determine the desired 

s 

angle p between the satellite spin axis S and the Sun vector U . Neglecting 
sensor mounting alignment errors, = 90® ~ ~ ^ ' 

2.2 HORIZON DETECTOR 


.Figure 2-2 show's the main geometric features of the horizon detector system 
which has been assumed in the design of MSAO/OABIAS. The sensor's FOV 
is narrow and usually' is either circular or square, hi practice, the FOV diam- 
eter or side is commonly in the 1- to 1. 5-degree range. A hypothetical unit 
vector L is located in the center of the FOV; £ is fixed in the satellite. The 

/V 

sensor mounting angle y bctw'cen the satellite's spin axis S and L is tailored 
to mission requirements. Satellites which employ horizon scanners commonly' 
include two or more units with different y aisles. This provides redundancy 
and permits inci*cascd covcj^age throughout the mission since the band swopt 
out bv the horizon scanner is small and w'ill miss the central body in some 
spacecraft orientations. 

The rotation of the satellite causes L to sweep out a small circle on the ccles- 
tial sphere, indicated in Figure 2-1 by the dashed line. The horizon detector 
system provides output signals (1) at *'in-timcs'' t^^^ , when L crosses the 
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Figure 2-2. Horizon Detection Geometry 
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boundary from the sky to the central body (Earth or Moon), aiKl (?• at "out- 
times" tjjQ , when L crosses the boundary from the central bodj'^ to the skj<’. 
There will be one in-crossing and one out-crossing in each spacecraft revolu- 
tion during those portions of the mission when the central body lies within the 
FOV loci On the celestial sphere (assuming intersection with only one central 
bodjr). The outputs of the horizon scanner system used by MSAD/OABIAS are 
the horizon crossing times SlI and tjj^ . In most systems, the horizon 

crossing parameters measured by the onboard equipment actually are t - t 

HJl 8 

and t,,^ - t . Ground (derations convert these measurements to the inputs 
HU S 

^ and tjjQ required by MSAD/OABIAS. 

Horizon detectors arc desigpicd to be sensitive to either visible light or to in- 
frared radiation. For MSAD/OABIAS operations, the significant difference 
between the two tjqies of sensors is that visible light sensors detect planetary 
disk terminators, while infrared sensors detect only the true planetary disk 
boundary. The OASYS portion of the MSAD/OABIAS System is capable of de- 
tecting and rejecting terminator crossing^. 

2.3 AE-C WHEEL-MOUNTED HORIZON SENSORS 

The AE-C spacecraft has one horizon sensor mounted on its body and two sen- 
sors that cffectivclj' are mounted on the spacecraft's momentum wheel, hi 
fact, the two w'hcel sensors arc mounted on the body of the spacecraft with their 
field of view' nominally parallel to the spin axis of the momentum wheel. Each 
of these sensors looks into n mirror mounted o.. the momentum wheel so that 
the motion of the wheel carries the sensor scan about the celestial sphere. 

Within OABIAS, the modeling of the wheel sensors is identical to that of the 
body sensor. Since there is no wheel-mounted Sun sensor, only Sun ai^lc data 
(from the body-mounted Sun sensor) and Earth width data arc available for 
attitude determination. The modeling of the wheel-sensor biases is the same 
as for the body-sensor biases. Ilowxwor, because the bolometer associated 
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with the wheel-mounted sensors is mounted on the body of the spacecraft, an 
additional physical misalignment is possible. Specifically, the bolometer coiild 
be mounted off-axis of the bod^' ot the spacecraft and, therefore, be misaligned 
with tlie axis of the wheel. This would cause a sinusoidal oscillation in the 
wheel-sensor. Earth v.7dih data as tlie spin of the spacecraft carried the wheel 
sensor axis in a small circle about the spin axis. The modeling appropriate 
to a misalignment of the sensor axis is discussed in detail in Reference 9. 

^.4 RAE-B PANORAMIC ATTITUDE SENSOR' (PA^ 

Figure 2-3 shows the geometry of the PAS which was used on the RAE-B satel- 
lite. ha this figure, and Ay are small misalignment aisles which will be 
discussed in detail in Section 3. 

One difference between the PAS system and the usual horizon detector is of 

primary significance to OABIAS; namely, the mounting angle y oi the ^AS 

s 

line-of-sight vector L relative to the satellite's spin axis S is not constant, 
bstead, y can be stepped in a predetermined manner. 

On RAE-B, the PAS sensor system supported the translunar phase of the mis- 
sion during which the central body — the Moon — was sufficiently small to be 
approximated as a point source for attitude determination purposes. The 

migle y was advanced 0. 7 degree per spin revolution. These steps were 
s 

actuated by reference pulses from a Sun sensor. The y advance signal was 

s 

Inhibilcd upon scanner acquisition of a central bod}'. The angle y at which 

s 

the central body was acquired was the primal}' information which PAS provided. 

The time variation of y on the PAS system made the addition of the misalign- 

s 

mcnl angle to the OABIAS mathematical representation of the horizon 
detector necessary. Wlicn modeling the more common body-fixed horizon 
detectors, is not needed. IVo new obset*valion models — 7 and 8 — were 
added to OABIAS for use with P.AS. 
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NOTE I: AND Ay INDICATE MISALIGNMENT ANGLES. 

NOTE 2: VARIES IN 0.7 DEGREE STEPS. 


Figure 2-3* RAE-B PAS Geometry 
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3.1 INTRODUCTION 

• ♦ • T 

Tills section contains the matlicmatical algorithms used in the recursive esti- 
mator subsi'stem (OABIAS) of tlic MSAD/OABIAS System. The log^c used in 
OASYS, the deterministic attitude determination subsystem, is sununarized 
in Section 3. 2. Section 3. 3 describes the general recursive estimation tech- 
nique and its application to attitude determination problems. The description 
covers the following three topics: 

1. The recursive and batch processing approaches with a comparison 
of their merits and drawbacks. 

2. The fundamental recursive estimation algorithm on which OABIAS 
is based. (See Appendix B for the derivation of this algorithm. ) 

3. The selection or computation of the observation w»eighting factors. 

The general features of the implementation of the estimation algorithm in 
OABIAS are discussed in Section 3. 4. This section describes the modifications 
made to the algorithm presented in Section 3. 3, denotes the main inputs which 
the program requires, and presents an introduction to the observables and state 
vector elements which OABIAS employs. Also discussed are the options in 
OABIAS w'hich enable the user to influence its performance as a recursive esti- 
mator, i.c. , the state vector updating frequency selection and the iterative 
operation options. 

The remaining portions of Section 3 describe the OABIAS mathematics in more 
detail. Section 3. 5 discusses the satellite and sensor geometry assumed and 
employed in ihe OaBIAS subsystem. This seetion includes coordinate frames, 
transform ;it ion matrices, and the mathematical representation of the sensors. 

It also presents precise definitions of the elements x in the state vector X . 
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Each of the eight observation models available in OABIAS arc discussed in 
Section 3. 6. (See Appendix A for the equations of the observation partial de- 
rivatives.) Sections 3.7 through 3.9 cover three additional topics: the com- 

putation of the central body angular radius p , the horizon crossing vector D , 

c 

and the weighting factors \v. . Finally, Section 3. 10 summarizes the Section 3 
material. 

The notation used throughout Section 3 is as follows. General 3x1 Cartesian 
vectors are represented by an overhead arrow (^) , and 3x1 Cartesian vec- 
tors of unit length are represented by an overhead caret ( V) . The magpiitude 
of a vector normally is designated by encasing it with two vertical lines (|^|) ; 
however, a simpler representation (V) is used in places where this can be 
done without ambiguity or confusion. Coordinate frame resolution is signified 

-=*CrI 

by superscripts (V ) . 


More general (n x i) coliunn vectors are designated by an underline . 

Matrices are signified by upper case letters with no underlines or overhead 

T 

symbols (P) . Superscript T designates the transpose of a matrix (G ) or 

T 

of a column vector (X ) . Where it is necessary to dis<^inguish betu^een ma- 
trices and scalars, the scalars are designated by lower case letters. For 
example, the elements of matrix (A) are signified by lower case letters with 
subscripts (a^.) . The dimensions of vectors and matrices are indicated be- 
neath the symbols the first time they appear in equations. 


A state vector estimate is indicated by an overhead caret (X) to distinguish 
it from the true state vector (X) . This should cause no confusion with the 

/V 

unit vector notation (V) , because the correct meaning should be obvious from 
the context of the equations. 
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The notation ilisitinj:;;uishcs between variables and algebraic expressions which 
arc used to compute tlie values of these variables. For example, in the cqu 'tion 

- • 

= Y^ (M. X, t) 

Xp » 2$ » t) is a known algebraic expression. The numerical values of the 

elements of Y arc computed by inserting values of ^ , X , and t into it. 
c 

Parentheses ( ) arc used in Section 3 almost exclusively in this manner. Also, 
brackets [] usually denote matrices and braces [}, vectors. However, all 
three — parentheses, brackets, and braces — are used to enclose scalars. 
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3.2 OASYS— I)i: ri:ilMINIS riC AITITITDE nETKRMINATION SUBSYSTEM 


The analysis required for an overall understanding of MSAD/OABIAS is sum- 
marized in this section. (See Referoneo 8 for a detailed description of the 
mathematical techniques used in OASYS. ) 

Attitude computations arc based on hori<^on crossings only; therefore, OASYS 
must identifi* and reject terminator crossings. The attitude calculation is 
basca on cither a single-horizon crossing method or a double-horizon crossing 
method. OASYS can process horizon crossings occurring on either the Earth 
or the Moon, and can determine which central body is bc’ng observed, if nec- 
essary. . In most cases the system mil function without the benefit of an initial 
attitude estimate. 


3.2.1 Input to Deterministic Process 

Each input frame of preprocessed telemetry is assumed to contain the following 
information from a single spin pcr.od of the spacecraft: 

o The Sun angle, or angle between the spin axis and the Sun vector 
o The Sun sighting time 

• The horizon sensor "on” time 

• The horizon sensor "off" time 

« The spacecraft spin rate 

OASYS processes each frame of data to identify and reject any horizon sensor 
triggering which occurred at a terminator crossing. For each remaining hori- 
zon sensor triggering, the system computes up to two possible dctcrndnistic 
attitudes. After a block of tlata has been processed in this fashion, a block- 
averaging technique is used to select from each pair of attitudes a single atliludc, 
such that all chosen allitu<les remain approximately constant throughout the 
block, '^hc chosen attitude vectors arc then averaged. 
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3. 2, 2 Design Assumptions 
The following assumptions are necessar 3 ' for deterministic processing: 

1. The spacecraft attitude is assueicd to be approximatcl 3 ' constant 
during one spin period. Note, however, that the spacecraft position 
is not assumed to be constant during one spin period. 

2. The spacecraft spin rate is assumed to be approximately constant 
during one spin period. 

3. The spacecraft is assumed to be sufficiently close to the Earth so 
that the vector from the Earth to the Sun is approximately parallel 
to the vector from the spacecraft to the Sun. 

The quality of the input sensor data and the accuracy of the initial attitude esti- 
mate determine the system logic used dm*ing any particular execution. For 
example, the status of the terminator must be checked for each scan. OAS^' S 
must handle the following cases: when the Earth is fully sunlit, when the ter- 
minator is visible but is not present in the scan, and W'heu the terminator is 
intersected by the sensor scan. In all cases the logic is based on the compu- 
tation of attitude from a single-horizon crossing, with either one attitude 
computation or two, depending on whether the terminator was present in the 
scan. The ambiguities involved in calculating attitude from a single-horizon 
crossing arc eliminated by logic spanning successive data points. The as- 
sumptions arc made tlvit the satellite attitude i.^ changing slowly and Uiat suc- 
cessive nadir angle calculations should 3 'icld an unambiguous attitude. The 
presence of the terminator in the scan is dctcrmiricd b 3 » a recursive technique 
utilizing data predictor modides. Singular conditions, which yield ambiguous 
results for attitude over a short time span, exist both in the case of a sunlit 
Earth and in the ease of the presence of the terminator, hi the latter case, 
an a priori attitude must be used to resolve the ambiguit 3 % 
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3.2.3 Dotcrininistic T^oiyic 

3. 2. 3. 1 Single Frame Processing 
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Attitude computation is chai*actcrizcd b}* three levels: the level of a single 
telemetry frame containing one or two valid horizon crossing triggering times, 
the higher level on which a full block of processed frames is averaged by the 
system, and the lower level on which the single triggerings are processed. 


For each central bod^' under consideration (Earth, Moon) , the following func- 
tions are pci'formcd at the intermediate, single telemetry frame level: 

1. The lighting conditions on the central body are computed from 
ephemeris data. If the central body is dark or not visible, this 
triggering is rejected. 

If a terminator is visible, the terminator flag for this frame is 
set. If data are not to be included in attitude calculations while 
fl terminator is visible, the triggering is rejected; otherwise, 
processing continues "s for the sunlit case. 

2. The attitude '^termination routine, ATTOET, is called and as 
many as two attitudes are computed for this crossing. 

3. If ATTOET produced no solution, there is no possible attitude con- 
sistent ^vith the assumption that this triggering resulted from a sun- 
lit horizon crossing on this central body; therefore, the crossing 

is rejected. 

4. If the central body is fully sunlit, the processing of this crossing 
is complete. If a terminator is visible, it must be determined 
whether this triggering resulted from a terminator crossing or 

a horizon crossing. 


f 
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If a tci'minator is visible, it must be determined whether f* 's intersected by 
the spaceex'aft sensor scan, i.c, , whether a particular triggering of the sensor 
was in fact a hoi'ison crossing or a terminator crossing. The procedure for 
• this determination is a recursive one, utilizing ODAP (data simulator subsys- 
* i) (Reference 7) modules as tools. (Sec Reference 5 for background analysis 
: ODAP. ' 

The routine TERCIIK is called once for each possible attitude computed in sub- 
X itine ATTDET, TERCHK calls the ODAP module with the computed attitude 
and determines whether a scan of the central body with this attitude would have 
produced a sunlit horizon crossing for the in or out triggering, as required. If 
the computed attitude is not consistent with the assumption that this triggering 
occurred at a sunlit horizon crossing, this attitude is rejected. 

Note that it is possible that this test wiP fail to reject a terminator crossing 
when the attitude computed from a terminator crossing is so far from the true 
attitude that a scan with the erroneous attitude would gdve a sunlit horizon 
crossing at this triggering. VVlien this occurs, one attitude is consistent with 
the assumption that tliis was a terminator crossing, and a second attitude is 
consistent with the assumption that the triggering was a sunlit horizon crossing. 
Therefore, there is no deterministic procedure for recognizing this problem. 

However, when the problem occurs, the resulting computed attitude generally 
has a large error. Since the error is large, the erroneous attitude is easily 
recognized and rejected in the block-averaging module, SPINAV, assuming 
either an a priori altitude or a large block of data is available. 

3. 2. 3. 3 Single-Horizon Crossing Computation 

Single-ljorizon crossing events are processed in the module ATTDET. As many 
as two possible alUludes arc calculated for each event, each with corresponding 
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nadir angles and dihednu .ingles. ATTDET is the key module in the inlerfranic 
processing. (See Reference 8 for a complete discussion of this computation*) 

3.2.3. 4 Block- Ayci'aging Procedure 

Once each input .elenictiy frame has been processed singly, the best estimate 
of the attitude must be computed based on the single frume results. Each input 
telemcti'y frame results in two output measurement frames, each of which may 
contain zero, one, or two attitudes. The ambiguities \vhich could not be re- 
solved on a single frame basis now can be eliminated if the block of data 
large enough. 

The module SPINAV processes a block of output frames to resolve the ami. ki - 
lties and determine the awrage attitude. For each output frame containing twv 
possible attitudes, SPINAV must determine which of the attitudes, if any, is 
valid. Note that in some cases both attitudes must be rejected as err^'ieous. 

The asrumpti^’-n Inherent to SPINAV is that attitude should remain apprcvlmately 
constant over the duration a block of data. Therefore, the ambiguities should 
be resolved in such a w’ay that the selected attitudes are as nearly constant as 
possible. The following procedure for resolvi;.g the ambiguities applies to the 
most general case, in which each output frame contains two attitudes and no 
a priori attitude is available. 

The first attitude from the first output frame is selected as a trial attitude. In 
each other output frame the pair of attitudes is examined, and the attitude from 
that pair which is closer to the trial attitude (in degrees of arc-length) is se- 
lected. One attitude is selected in this manner from c.*:ch output frame. 

The set of altitudes thus selected Is avemged using SPNAVl. The attitudes 
arc expressed as unit vectors, each component of the vectors is averaged, and 
the resulting vector is unitized. (If data weighting is employed, a weighted 
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average is computed, with each weight equal to l/(attitudc uncertaintj' in de- 
2 

grccs of are-length) “.) The resulting unit vector is the average attitude for 
tho set of attitudes selected. 

• r 

However, the set nia\' include cri*oucous attitudes, resulting from noisy or 
biased d:ita, or from terminator erossings wliich could not be rejected. There- 
fore, a residual edit must be performed within SPNAVl. First, the standard 
deviation of the set of attituckis is computed as follows: Let the residual p. for 
each attitude be defined as the angle in degrees between that attitude vector and 
the average attitude vector. Then the weighted standard deviation a is com- 
puted as follows: 



where o). denotes the weight for the ith frame and the summations are taken 
over all frames in the block. 

The resulting value for a is the standard deviation of the set of attitudes, in 
degfrees of arc-length. Next, each attitude in the set which has a residual p^ 
gpreater than na , where n is an input parameter, is rejected. The default 
value of n is 3 . The attitudes not rejected in this process are then re- 
averaged. A new standard dcvia<^ion is computed and the entire procedure is 
repeated until no further rejections occur. 

Then SPNAVl returns witli the block-average attitude, the standaitl dc\iation, 
and Uic number of frames included in the final average M . Next, SPINA V' 
computes a got dness-of-fil jviramelcr, which is intended to measure the amount 
of scatter in the sot of selected attitudes. The goodness-of-fit piiramctcr is 
equal to a/M . Note that this parameter takes into account both the final stand- 
ard deviation and the number of frames inoluded in its computation. 
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The above procedui'c describes the steps followed using the first attitude from 
the first fniine as a trial attitude. This procedure is repeated usii^ each atti- 
tude from each frame as a trial attitude. A running compel rison is maintained 
on the goodness-of-fit parameter, and at each step, the block-average attitude 
associated with the best fit is saved. 


The average attitude associated with tlie best fit is refined further bj'-thc follow- 
ing iterative technique. Using the average attitude to resolve the ambiguities, 
the program selects one attitude from each pair, as above. SPNAVl Is called 
to awrage the selected attitudes. The a\'erage obtained by SPNAVl is used as 
a new average and the sequence is repeated. The process terminates ivhcn the 
set of attitudes selected remains identical for two successive iterations. Con- 
vergence normally occurs ia two or three iterations. 

The following discussion is intended to clarify the theory behind the block- 
averaging procedure and explain its strengths and weaknesses. 

Consider a block of N output frames, with each frame having two attitude 

N 

solutions. The ambiguities can be resolved in 2 different ways, if one 

attitude is selected from each pair. If a possibiliiy exists that neither attitude 

N 

in the pair is correct, then there are 3 different ways to resolve the am- 

big^tics. A goodness-of-fit parameter can be defined which measures the 

amount of scatter in the selected attitudes and defines the "correct" way for 

..esolving the ambiguities as the way which minimizes this parameter. An ex- 

N N 

haustivc search of the 2 (or 3 ) choices ir the only procedure guaranteed 
to 5 'icld this correct solution. Clearly such a procedure is impractical unless 
N is very small. (Note that if N - 2 , the procedure is quite practical. In 
fact, it corresponds to the case of two horizon crossings, cacli yielding two 
attitudes, in tli:it case, a check of the four possilile pairings will reveal whicli 
pair of altitudes is optimal. Wlicn N - 2 , the method in SPINAV reduces to 
the same situation.) If N is large, some simplifying assumption must be 
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introduced. The method in SPIMAV assume:, that at least one of the attitudes 

in the bloek is elose to the fimil desired block average, i.e. , the block ;iveragc 

wiiich would Im ..I'lained if tlie ambiguities were correctl 3 ' resol\t!d. Clearlj', 

it would be vc-rv unusual for a set of attitudes to 3 'icld an average which did not 

lie elose to an>' aititU'Jc in the set. Therefore, this simplifying assumption is 

reasonable, and in virtually all cases it should yield the same result as an ex- 

N 

hausti\*e scai'ch of all 2 cases. 

The computation ti:nc required for SPINAV is approximately proportional to 
. (The number of trial attitudes is 2N . For each trial attitude, SPNAVl 
requires a computation time proportional to N ,) Therefore, the block- 
averaging method is impractical when N is large. In practice it has been 
found tliat a block with 200 useful output frames requires several minutes of 
processing time on the IBM S/360-95. In most cases such exhaustive proc- 
essii^ of a large block is unnecessary. A small subset of the block can be 
processed, and the average of this block can be used as an a priori attitude. 

If an a priori attitude is available, it is used as a trial attitude, and the search 
for a trial attitude is eliminated. The short iterative procedure described 
abo\*e can be used without searching for a trial attitude. The results in most 
cases will be identical to those which would be obtained using the complete 
method. 

Finally, no method for resolving the ambiguities, not even an ejdiaustive search 
N 

of all 2 combinations, will yield correct results in all cases. If the attitude 
is allowed to change in an arbitrary maimer from one frame to the next, re- 
solving ambiguities is clcaii.v impossible, even if tlic daUi are perfect. Also, 
if the attitiKle is constant, but sj'stcmatic or r.andom errors in the data cause 
the "correct” altilutlo solutions to vary within the block b\’ an amount compa- 
rable to tiic tlifforence between the two attitudes in c:ich pair, ambiguities 
cannot be resolved. 
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3.2.3. 5 Other Dctci*ministic Atliledc Compubitions 

OASYS provides four methods of independently calculating attitude: one using 
the Earth-in triggering, one usirig the Eartli-out triggering, a third requiring 
both triggerings to calculate an Earth nidth angle and then a nadir angle, and 
a fourth using the Earth scan mid-time as computed from the two Earth trig- 
gerings. This latter method is the double-horizon dihedral angle method. (See 
Reference 8 for the mathematical details of tliese methods. ) 

3. 2. 3. 6 Data Weighting 

An option is provided in OASYS to compute data weights to be used in the bloek- 
averagpng procedure. These data weights are computed from user- specified 
values for the uncertainties in the observables, ff the data weighting option is 
selected, five uncertainties are computer for each single-horizon crossing solu- 
tion: the arc-Ie '•.ertainty in attitude, the uncertainty in right ascension 

of the spin vector ihe uncertainty in declination of the ;L'piu vector ^ , the 

uncertainty in the naai.- mgle, and the uncertainty in the dihedral angle. For 
each double horizon solution, the uncertainty in the nadir angle and the uncer- 
tainty in the dihedral angle are con 4 >uted. 

hi the block-averagii^ procedure (SPINAV), the unit spin vectors are averaged 
writh each spin vector assigned a weight equal to l/(arc-length uncertainty in 
degrees) . 

In addition to pro\iding weighting factors, OASYS provides a measure of the 
unccrlaintj' in the block-average attitude, by computing a weighted average of 
the unccrUiintics for all frames. A w'cightcd average is computed for the un- 
certainties in arc-length, a , and 6 . The weights used in computing these 

W'Cightcd averages arc the same weights used in the block-averaging procedure, 

2 

l/(arc-length uncertainty in dcgi'ccs) . 


3-12 


t*w*’ l'**» 



ORIGINAL PAGE W 
OF POOR QUALITY 


The following technique is used to compute unccrlaintics; Let f be a function 
describing some computed quantit 3 ' in terms of observed quantities. For ex- 
ample, f might be the function \vhich c^qpresscs the right ascension of a unit 
vector along the spin ;ixis a in terms of the obscr\'ablcs yS , y , P , A , and t . 

a =fC8, y.p, A. lf(t)) 


where = Sun angle 

Y = sensor mounting angle 
p = effective angular radius of the Earth 
A = rotation ai^le from Sun crossii^ to horizon crossing 
R (t) = spacecraft position at time t 

Strictly speaking, not all of these quantities are directly observed. For ex- 
ample, the Sun angle is computed from a coded field in the telemetry, and 
the rotation angle A is computed from a spin rate and crossii^ times, which 
in turn are computed from clock counts in the telemetry. However, this set 
of observables has the advantage that uncertainties in each observable may be 
conveniently estimated and treated as constant, at least over a single block of 
data. • i 

Standard error analysis yields the following expression for the uncertainty 
in a : 



where a = standaixl dexiation of or 
a . 

Oo = sUindaixi deviation of fi 
p 

and similarh’ for each observable. 
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This analysis is based on the assumption that the errors in the observables are 
normally distributed, uncorrelatcd errors which arc small enough so that the 
function f may be treated as linear, i. e. ,* the partial differentials do/hP , 
etc. , may be treated as constant. 

Given that the function can be treated as linear in the desired region, the de- 
rivatives may be computed numerically, as follows: 

.2 

^ Cf 03, y, ...) - f 08 + o^, y, ...)] 

This technique is easily applied to any function f which can be evaluated nu- 
merically. 

hi actual implementation the method works as follows: the attitude, the nadir 
angle, and the dihedral angle are computed using the observed data- Then the 
specified uncertainty in the Sun angle is added to the observed Sun angle, and . 
the attitude, the nadir angle, and the dihedral angle are re-computed, using 
the perturbed Sun angle ivith all other parameters as before. The changes in 
attitude, nadir angle, and dihedral angle are sa%'ed. Next the specified uncer- 
tainty in some other parameter is added to the observed value, and attitude is 
computed again, using a perturbed value for only one observable at a time. 

After each observable has been perturbed in this manner, the uncertainty in a , 
for example, is computed as 


a 


a 



where Aa. - the change in or resulling from a perturbation to the ith 
observable 
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Note that for the uncertainty in time t the ephemeris routines must be called 
to obtain the position at the pcrtui-bcd time t + o . This ertor in time can be 

If 

used to represent both the effect of an absolute timing error (i.e. , the absolute 
times fi*om the attitude telemetry arc not consistent with the absolute times 
used for orbit determination) and the effect of an in-track orbit error of a given 
time magnitude (i. e. , the spacecraft is 30 seconds ahead of or behind the orbit 
tape prediction). 

This method of uncei'tainty computation breaks down if attitude cannot be com- 
puted from the perturbed data (i. e. , if -the function f (fi + Oot 7* . . .) is unde- 

P 

fine<^. In this case, the function f is certainly not linear over the region 

P to P + Oa * because f is not defined over the entire interval. l\*hen this 
p 

condition occurs, OASYS assigns a very large uncertainty (99999. 0 degrees) 
to the attitude, nadir angle, and dihedral ai^le, resulting in a very low weight 
for the frame. This is a reasonable action in most cases, because the uncer- 
tainties are generally very large in the region close to the point at which attitude 
becomes undefined. That is, for the functions under consideration, the deriv- 
atives go to infinity at the point where the function becomes undefined. (This 
is true, for example, for the square root, arc sine, and arc cosine functions, 
all of which are involved in the attitude computations. These are also the 
only functions which can cause the attitude computation to be undefined. ) The 
problem of perturbed data yielding no solution can always be avoided by using . 
sufficiently' small values for the uncertainties in the observables. For example, 
if the uncertainty in each observable is reduced by a factor of 10, the relative 
weights will remain accurate, and the probability' that the perturbed data will 
yield no solution is reduced. The user must then remember to multiply the 
computed atlilude uncertainties by a factor of 10 to obtain realistic values. 
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3. 3 THE UECUHSr\^E ESTIMATOR APPROACH 

3. 3. 1 Comparison Rctwccn Recursive Processing: and Hatch Processing 
Estimation Methods 

This section compares, nonmathcmatically, the batch processing and recur- 
sive processing methods for estimating the time-invariant state vector X of 
a system using, as input, a block Y of noisy observations y. obtained on that 
system. The discussion of batch processors infers the usual least-squares 
differential correction (DC) algorithm employed, for example, by the GCONES 
program and GTDS. The discussion of recursive processors infers the usual 
nonlinear recursive least-squai'cs algorithm employed by OABIAS. 

With the batch processing algorithms, all observations y^ in Y are, in effect, 
handled simultaneously. Batch processors which take into account observation 
geometry nonlin3arity (to be deHned in Section 3.3. 2) require an a priori esti- 
mate X of X to start the operation. A batch processor does not provide a 
new estimate X until the complete set of observations Y has been processed. 

With the recursive processing algorithms, Y is separated into a number of 
mutually exclusive subsets. For example, in a system such as OABIAS which 
processes the observations one at a time, the subsets are the scalar observa- 
tions y^ . They usually are ari'anged in chronological order and are processed 
sequentially. Updated and improved estimates X are obtained continually 
during the processing operation* In the most common operating mode, X is 
updated 'ter processing each y. , and this updated value is employed in com- 
puting the partial doi*i\'ativcs G required for processing the next observation 

/N 

y^^ j . In other modes, X may be updated- less often, e.g. , as infrequently 
as after the full daUi set Y is processed. \Micn the observation geometry is 

/N 

nonlinear, the X uptialing frequency is of some significance, because Its in- 
flucncc on the Gs affects the dynamics of the filter. Recursive processors 

-A. 

always require an a pi*iori estimate to start the processing operation. 
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An iterative processing technique is one in which the observation set Y is 
reprocessed several times. The state vector estimate X obtained at the end 
of each complete pass of Y through the processor is used as the initial input 
for the following pass. Itcratiw processing is useful in improving the accuracy 
of the final estimate when the i-esult obtained by only a single pass would be 
degraded by the nonlinearity in the observation equations and the error in the 
a priori cstinuite . In practice, batch processors used for attitude deter- 
mination problems usually must bo operated iteratively to yield trustworthy 
results; the procedure is called differential correction (DC) . The nature of 
the recursive processing algorithms makes iterative operation of recursive 
processors less essential in attitude determination problems, because the 
estimate X is continually being updated. However, iterative operation of 
recursive processors is possible, often advantageous, and sometimes nec- 
esrary. 

Recursive estimation algorithms have a number of potential operational ad- 
vantages over batch processing algorithms. The main ones are as follows: 

1. Recursive processors handle only a small number of observations 
at a time. Therefore, because they need to store only a small 
number of observations at any one time, recursive processors 
can require less computer core space than batch processors and 
have no a priori limit on the size of the data arrays they can proc- 
ess. 

2. Recursive processors arc better suited for real-time opci*ations, 

because the operation need not be delayed until a block of obser- 
vations is accumulated. Individual obscrv.ntions can be processed 
immcdialely, thus providing immediate estimates of the attitude 
state of the spacecraft with various measures of the accui'aey of 
that estimate. f 
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3, The series of successive state esTimalcs provided by the recursive 
processor allow the operato' to watch the convergence of the solu- 
tions and develop an intuitive feeling for the quality of the solutions 
that is not possi’de with the single output of the batch processor. 


The potentiality for the real-time operation of a recursive processor is deterred 
when iterative processing is required. Also, both the real-time and core stor- 
age advantages of recursive processors can vanish if preprocessing of the input 
data is required or if storing the input data in large blocks is necessary rather 
than passing the data to the recursive processor on receipt. This has been the 
case wifli MSAD/OABIAS. Hence, the development of an intuitive feeling for 
the qualit}’’ of the solutions has been the main advantage of recursive processing 
over batch processing in the OABIAS application. 


The principal advantage of batch processors is stability. Because recursive 
processors can update their state vector estimate X continually (after proc- 
essit^ each y. , if necessary) , in iterative operation, they tend to converge 

rs 

to a final X estimate faster than batch processors. Because batch processors 

/s 

update X onlj' at the end of each iteration they are more likely to converge to 
a valid solution in difficult problems, i.e. , to be more stable, bi a recursive 
processing operation, however, updating X only at the end of each iteration 
will overcome this difference and yield dynamics virtually identical to that of 
batch processing. 


In addition to stabilit}', batch processors have hra other potential advantages 
over recursive processors The first is naming speed. Because batch proc- 
essors generally i*cqui»'e fewer numerical operations than recursive processors, 
to process a block of dati they tend to run faster. Secondly, batch processors 
can take into account the effects of correlated errors in y^ better than recui’- 
sive processors. This capability, however, is not utilized often and, therefore, 
can rarely be considered a significant advantage. 
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3. 3. 2 TIk; IiMsIc Hociirsivc KyCimnior AltyorHlim 

This section discusses the concept of observation equations and geometry non- 
linearity, defines some terminology to be used in the remainder of Section 3, 
and presents and discusses the fundamental x*ecursive estimation algorithm on 
which OABIAS is based. 

Recursive estimators require a set of observations Y as inputs. Let the in- 
dividual scalar observations . '' be designated as y. where j = l, 2, p. 

Associated wiUi each observation . j is an algebraic expression y . (X , t) 

cj ~ 

which enables the value of y to be predicted, y . (X , t) is a mathematical 

] CJ 

model of the satellite and its pertinent sensors. Thus, for each observation j 
there is an equation of the form 


y, = y„j S. tj) + vj 


(S-l) 


where t^ is the time at which observation j was obtained, and y^ is the actual 

"measured” value of the observation, v. is an error term which must be in- 

J 

eluded to make the two sides of the equation balance, v^ results from the error 
in the y^ measurement and from modeling errors, i. e. , from approximations 
lu the y^j (X , t) expression. The actual value of v^ of course is unknown. 

In some problems Equation (3-1) can be placed in the form 




(3-2) 


whci*c ji. and b. arc constant or time-dependent coefficients, but arc not 
explicit functions of X . 
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In tins case it is said that the observation gcomclry is linear. Wlicn Equa- 
tion (3-1) cannot be placed in the above form it is said that the observation 
geometry is nonlinear. Recursive estimation problems in which the observation 
gometry is linear tend to pi'cscnt fewer difficulties than those in which it is 
nonlinear because there is a fully developed and rigorous body of theory for 
generating mathcmaticallj' optimal results in the linear case. Unfortunately, 
tlie observation gcoractry in attitude dctermir.atioa problems usually is suffi- 
cientli' nonlinear that its effect must be taken into account in the basic algo- 
rithms of the recurs’ve S 5 'stem and/or compensated for by ad hoc procedures. 


It is necessary at this point to define some ..eims which will be used in die 
remainder of Section 3. This terminology is nearly identical to that used in 
Reference 1 . The word "measurements" will oe used when referring to the 
independent inputs m^ supplied, via telemetry and preliminary ground proc- 
essing, to the overall attitude determination system. To be specific, the 
measurements supplied to MSAD/OABIAS are the Sun angles , Sun sighting 
times tg , central body-in horizon crossing times t^^j , and central body-out 
horizon crossing times t„_ . The word "observables" will be used when re- 
ferring to the basic vai .yles which serve as inputs to the recursive estimator 
portion of the attitude determination system. The terms "real observations" 
or merely "observations" y^ will be used to designate the measured values 
of the observables. The term "model observations" will be used to designate 
the predicted values y^, of the observations. The word "model" is used here, 

because y , is obtained using mathematical models of the satellite 'jnd its 

cj 

sensors. In Equation (3-1), y , (X , t) Is the algebraic expression which is 
used to compute 3 ' . 


In most estimation studies, a distinction is not made between measurements 

m and observations 3 ' . Instead, the observables arc considered to be the 

a j 

basic parameters which are measured. A distinction is being made in the 


•. * 
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present analysis, however, because the observables used by OADIAS (and also 
those used by OASYS/GCONES and OASYS/GllECRS) arc not the basic measured 
parameters t^ , t^j^ , and . The direct use of event times, such as t^ , 
SlI * and t^j^ , as observables is difficult because of the difficulty in predicting 
these times via algebraic models as is required by Equation (3-1). The easier 
and more usual approach is to convert t^ , t^^^ , and tj^^ (via preliminary 
processing) into parameters which can be handled more easily as observables 
by the estimator portion of the sj^stem. Fox example, in OASYS/GCONES and 
OASYS/GRECUS the observables are Sun angles, nadir angles, and dihedral 
angles. (The observables used by OABIAS will be discussed in d; tail in Sec- 
tion 3. 6. ) One of the drawbacks of transforming measurements m^ into new 
observables is that it tends to increase the statistical correlation betu'een the 
observables. This difficulty occurs in OASYS and GRECRS which process the 
observations one at a time. 

The equations of the basic recursive processing algorithm used by OABIAS 
follow. (A derivatiop of these equations is given in Appendix B.) 


^cjR " ^cj ®jR^ 


(3-3a) 


%R = ^j V 


(3-3b) 


*jR = ' ^ciR 


(3-3c) 




I 
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SXS 

(3-3c) 

• • • 

AX. = kJ?.., - G^Jx. , -X.„! 

J j( JK “jR|( 

(3-30 

/\ /N 

X = X. . + AX, 
“j-1 ~i 

(3-3g) 

SXl 



where 


ayo,® 

8X1 


(3-41 


The above equations are applicable only to estimators, such as OABIAS, which 
(1) assume that X is constant in time and (2) process the scalar observations 
one at a time. Equatl<;ns (3-3Q and (3-3g) are not identical to the equations 
actually implemented in OABIAS; the modified eqxiations used in OABIAS are 
discussed in Section 3.4.4. The Equations (3-3f) and (3-3g^ aie used in this 
Introductory' discussion because they are more basic and easier to comprehend 
than the corresponding OABIAS equations. 


The computations denoted bj' Equation (3-3) are performed sequentially on each 

observation y. . In otlicr \rords, y is run through the equations, then y and 
j * . “ 

so on until the complete observation vector Y has been processed. Initial 

values of and must be provided for processing observation y^^ , (X^ 

is used as the initial reference vector .) In addition, a wcigliting factor 

Wj must be provided or computed for each y^ . The computations must be 

performed in the order shown in Equation (3-3) except ic innovative 
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rcsidunl can be computed at anj' point prior to the £iX^ computation. 

Also, the covariance matrix P, can be computed at r . point after calculation 

of the gain vector K. . P. is not used in calculr , but is used instead 

~J J . 

in processing the next j + 1 observation. 

The individual equations of Equation (3-3) now will be discussed, hi Equa- 
tion (3-3;'' is the aforementioned model observation. The subscript R 

* cjU 

in Equation (3-3) indicates that the parameter is computed using a reference 
value X. of the slate vector. In most estimation systems, is the pre- 
vious estimate • QAB?\S, however, has the capability of updating the 
reference vector less often than this; this capability will be discussed in Sec- 
tion 3. 4. 

hi Equation (3-3b), is an s x 1 partial derivative vector where s is the 

dimension of the state vector X . Its value is computed by inserting X.^^ into 

the algebraic expression G. (^ which is fornxed by differentiating the y (^ 

J c ■" 

eiqiression with respect to X as indicated in Equation (3-4). Essentially, 
and G^jj are the first two coefficients of the Taylor series expansion of the 
observation Equation (3-1) ; i. e. , 


y = y + G {X - X 1 + v. + higher order terms in (X - X ) (3-5) 

J cjtt JK JK J j.*t 


hi Equation (3-3c), z^j^ is called the residual, or innovative reaidual, of obser- 
vation j . z ^ provides an indication of the deviation of the true state vec- 
jR 

tor X front tlic reference vector . This can be seen more clearly by 
combining Equations (3-3c) and (3-5) to 3 rield 


Zjjj = G.*j^ [X - X^j^) + V. + higher order terms in (X - X.^^) (3-6) 
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In Kqurition (3-30 , it is seen that, to first oi*der, is converted to an ^ 
reference prior to its use in upchiting the state vector. In oOicr words, the 
terms within the outer braces I } on the right side of Equation (3-3f) were ob- 
tained bv truncating the scries 


yv *1/ /V 

s. = 2j {X^ - X^j^) + higher oi-der terms (3-7a) 


where 


GiR = 5,(XjR»=-^(X,/ 




(3-7b) 


Equations (3-3f) and (3-3^ show that the updated state vector estimate X. is 

>N /V ^ 

obtained by adding a correction vector AX* to the previous estimate X. ^ . 

^ -j ^ 

AX. is the product of the modified residual wetor z. (X, ,) and a gain vector 
Kj • K. is computed using Equation (3- 3d) ; the computation also requires a 
matrix P. ^ which was computed, via Equation (3-3e) , when processing the 
previous observation j - 1 . Except for degradation due to observation ge- 
ometry nonlinearity, establishing the gain vector through Equations (3-3<^ 
and (3-3e) is an optimal technique because the resulting estimate is opti- 
mal according to Uic several statistical and nonstatistical criteria. With the 
nonstatistical approach taken in Appendix B, this method is optimal because 
the resulting X. minimizes a generalized least-squares loss function. 

The term w of Equation (3-3d) tells tlic processor how heavily to weight 
observation j in generating the new estimate X^ . The matrix ^ , which 
was compiiUxI \ ia Equation (3- 3c) when processing observation j - 1 , tells 


the processor bow much weight to attach to the preceding estimate X 


j-1 * 
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Statistical considerations indicate that optimally (1) the initial input matrix P 

should be the covariance E(X - X^j - XqJ of the uncertainty in X^ and 

(2) each should be the inverse of the variance Oy^ in y* due to random 

errors in t’ne basic measurements m . If these conditions are satisHcd for 

a 

each observation yj , the resulting matrices Pj will be the covariance 
E{X - Xj} lx - of the unceitainty in the estimate Xj . Analyses which 
lead to these conclusions usually assume tliat (1) the obser\'ation geometry is 
linear, (2) the y. errors are uncorrelatcd, and (3) the input measurements 
enter directly into the y. , not the , computations. 

3. 3. 3 Discussion of Weighting Factors 


The question of obserration errors, optimal weighting factors w. , and the 

conditions under which processing the observations one at a time is an optimal 

technique will be considered at this point. The residts of this development will 

be used in the discussion of the OABIAS observables and weighting factors (see 

Section 3. 6). Let M be the composite n x i measurement vector of a block 

of data and let m^ , where a = 1 . . . n , be the individual scalar measurements; 

i.e. , m is an element of M. Let dm and dM be the scalar and vector 

measurement errors. Let Y be the p x i observation vector. Neglect flie 

degenerate case in which each element y. of Y is identical to a corresponding 

scalar measurement m . It is assumed that each observation y. is computed 

a j 

using one or more ; i.e. , 


yj=yj(M) (j = ltop) 


(3-8) 


Let dv. and dY be the scalar and vector observation vectors. Assume that 
■ J ~ 

they result entirely from the measurement crroi*s dM . 
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dY = H • dM 
pxi pxn nxi 


where 




dm. 


(M) 


(3-9a) 


(3-9b) 


Let R-, and R„ be the covariance matrices of dM and dY . From Equa- 
M Y _ _ T 

tion (3-9a), to first order 


pxp nxn 


(3-10) 


Statistical approaches to the estimation problem indicate that processing the 

observations y. one at a time, as has been assumed in the current section, 

can be optimal only if the dy 's are statistically uncorrelated; i.e. , if R 

J ' Y 

is diagonal. In this case, the inverses of the diagonal elements of R^ are 
the optimum weighting values w* . 

Using Equation (3-10), the neccssarj'^ conditions for a diagonal R can be 


shown to be that (1) each scalar observation y. is computed from its own set 
of measurements M. whose elements m. , w'hcrc a= 1, 2, .,,, n, , are 
not used in (lie comi)ut:i(ions of any other observation y , where k / j and 

Iv 

(2) the measurement subsets M , M^, ..., M are statistically uncorrclated. 

2 1 - P 

llicn, the variance Cw and optimum weighting factor wt of observation j are 

■’J J 


*-l 2 , T „ 

w* =0 -h R h. 
j y^ -j M. -J 
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where 


h. =T— ^(M.) 
ja Sm. ~j 
ja 


In the special case where the elements of M. are uncorrelatcd with variances 
2 

Cijj. , Equation (3-1 la) reduces to 
jot 


^-1 2 ,2 2 
wf = <r = > h. a 
i y Z-* J0£ m. 

^ 0 C =1 jc 


(3-llc) 


J» ft r 



3.4 impu:mkntation of the recursive estimator ai.goritiim ri 

OABIAS 

This section discusses the general features of the implementation of the re- 
cursive estimator algorithm in OABIAS, including the main options which arc 
available to ii\flucncc its performance as an estimator. 

3.4.1 Principal Inputs 

Preprocessed telemetry information is the principal input required by OASYS/ 
OABIAS. Each frame in this data set must contain the following information 
from a single spin period of the satellite: 

• The measured angle jS between the Sun vector and the satellite's 
spin axis 

• The time t at which the Sun was sighted 

s 

• The central body- in crossing time t ; i. e. , the time at which 

HI 

the line of sight of the horizon detector crossed the s!qr-to-central- 
body horizon 

• The central body-out crossing time t^^ 

• The satellite spin rate o) 

P t t *HI * ^IIO elements m^ of the measurement vec- 

tor M noted in Section 3.3.2. Thus, if there arc n useful frames in the 
telemetry data set and none of the measurements in any of the frames is dis- 
carded, in the dimension of M is 4n . In this context, the spin rates cc 
in the data set arc not considered mcasurcnicnts because they norma 11 j' are 

computed algebra icallj- using the t mcasuremcnls and, mere significantly, 

s 

because they arc used only in preliminary pi'occssing in 0.4SYS, not in OABIAS. 

In addition to the telemetry inputs M , the OABIAS recursive estimation sub- 

/nGI 

system requires the unit Sun vectors in geocentric Inertial (GI) coordinates U 
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at each Sun sighting and the ccntral-bod 5 '- to- satellite vectors In GI coordinates 
— ^GI 

U at each horizon crossing. These arc obtrined from cphemcris data or an 

orbit generator. The OABIAS repursive estimation subsystem also requires 

the angular radius p of the ccnti'al bodj' on the unit celestial sphere at the 

c 

horizon crossings, is computed by the program from the orbit data. In 

addition, the nominal wine of the horizon detector mounting angle y and the 

s 

initial estimates X of the state vector X and of the covariance matrix 
^ T 

= E (X - X^) (K - X^) are required. Pro^'am operating instructions and 
error parameters wliich are supplied by the user through the NAMELIST 
are also needed. 

3.4.2 Obscr\*ation Models 

OABIAS uses eight different observables, commonly referred to as "models. " 
(See Section 3. 6 for a detailed discussion of the mathematics of the eight 
models.) When running OABIAS, the user has the option of selecting which 
models are to be used. These models are analogous to and replace the Sun 
angle, nadir angle, and dihedral angle models employed in GCONES (Refer- 
ence 8) and OASYS/GRECRS (Reference 10). 

OABIAS possesses models with error-free real observations — a feature which 

is not commonly found in recursive estimation systems. In the usual estimation 

system, the measurements M and their errors enter into the computations 

through tlie real observations Y ; i.e, , Y = Y (M) and Y = Y (X) In 

_ _ — — c *~c ~ 

OABL\S Models 2, 3, and 4, however, the real observations arc error-free 
and M enters solely through Y^ ; i.e., Y^ = Y^ (M , X) . The Appendix B 
least-squares derivation of the OABIAS recursive estimator algorithm shows 
that this unoidliodoxy docs not alter the validity or optimality of the algorithm 


^Rcfcr to Section 3.3.2 for the distinction between measurements, real obser- 
vations, and model obscrv'ations. 
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because the cssenlial pi'oiierlies of the computed i*esidual vcctoi- Z arc re- 
tained. 


3.4.3 State Vector Elements 


The state wclor X used by OABL\S contains 12 elements x 
ments with their alternative symbols are defined as follows; 


V 


These cle- 


• X (s ) and x (s ) which define tlie attitude of the satellite; i.e. , 

X 1 £m 

the orientation of its angular momentum vector 

* X (0 ) which defines the phase of the satellite in its spin cycle at 

«j o 

the start of the run 


• X (ce) which is the satellite's spin rate 

• x^ (Ay) , Xg , Xg (Cjj) which define the effccove 

mounting alignment of the horizon scanner 

• X- (A^) and x ^ (c) which define the effective mountii^ alignment 

o 10 

of the Sun sensor 


• x^ (AP) which is the effective error in the central body angular 

radius p computed in OABIAS using orbit information 
o 

• x ^2 (At) which is a timing-bias due to an effective error in the 
computed location of the satellite in its orbit 


The precise mathematical definitions of the 12 state vector elements arc given 
in Section 3. 5. Note that and 0^ arc the cffcctiTC azimuth angles of the 
horizon scanner, relative to a body-fixed reference frame, at the Earth-in and 
Earth-out horizon crossings, rcspcctivelj'. OABIAS considers these as distinct 
pai'amctcrs. Also, is an alignment error of PAS and is not required with 
conventional horizon detectors. 


The initial cov. 'ianc . n.atrL\ which the user sui)|)lics to OABL\S is diag- 
onal. \Mien selected diagonal dcinents of this matrix arc set to zero, the 
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corrcsiwndin" olcnienls of X rcmn in constant, at theix* a priori values, 

/S 

throughout the run. Tims, the user can select the components of X to be 
updated. This is a powerful and axeccssary tool for the proper, use of OAIIIAS, 


OABIAS also provides the user with the option of not updating the state vector 

/S 

estimate X after each observation. Equation (3-3) shows that a reference 
vector Is used in the processing operations performed on each observa- 

tion j . X. enters into the mathematics primarily through the partial de- 
rivative vector . Most recursive processors which model the observation 
geometry nonlincax'ity (i.e. , which make each a fimction of ^update X^^^ 

at each observation using X._ = X. , . In OABIAS, hoxvever, X.„ may be 

~jR ~j-l -jR 

updated after eacn observation, after each telemetry frame, or after every 
N telemetry frames, where N is a user- selected integer.. These options 
have been included in OABIAS to provide additional flexibility of operation, 
particularly when observation geometr 3 ' nonlinearity is significant. OABIAS 

/S 

computes an updated state vector estimate X only at those points where the 

reference vector (now to be denoted as X ) is updated. For this reason, pre- 

R 


vious CSC reports have called the operation state vector updating rather than 
reference vector updating. The reference vector is updated by setting it equal 
to the updated state vector. 


To derive the state vector updating equations implemented in OABIAS, combine 

£ 0 iUi‘.kJt;S (3- 3f) and (3-3g) into a single equation. Replace X by X and 

**jR R 

repK'.ec the i subscripts by k*s . The first observation after the most recent 
upd.ati'ig ir, signified by k = 1 . Subtracting X^^ from boUx sides of the equation 
yields 







/N I 


(3-12a) 
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where, by definition 
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/S /V 

Tk ”k ”R 




(3-l2b) 


and 

= 0 (S-12C) 

It should be noted that the of Equation (3-12) is not the same variable 
as the of Equations (3-3f) and (3-3g) beeause it is defined relatiw to a 
different referenee. OABIAS sets AX equal to 0 when the state vector is up- 
dated and then processes the observations sequentially using Equations (3-3a) 
through (3-3e) and (3-12a). Let the number of obser vations between updates 
. be n . When k = n , OABIAS updates X and using 

5c = 5c + AX (3-13a) 

-n -o -“n 

% 

X„ = X (3-13b) 

“R -n ' 

where X is now Uie previous state vector estimate, 

“O 

3,4,4 Iterative Operation 


When the effect of observation geometry nonlinearity is significant, the per- 
formance of a recursive filter can often be improved by iterative operation. 
With this techniciue, the composite block Y of observations, or subsets of 
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this block, arc run through the filler several times. The state vector csti- 

mate X ol)tainec! at the end of any one pass \ is used as the a priori input 

X * X 

X for the next pass. If the procedure is successful, the estimates X, 

— o 

will converge toward a constant value as X increases. 


Recursive filters require, as inputs, not only a priori state vector estimates 

X^ but also an a prioi'i covariance matrix Pq . For iterative operation, a 

technique for establishing the mati'ix input to be used at the start of each 

o 

pass X must be decided upon. The two simplest approaches are; (1) to reset 

P to its original a priori value {i.e., P^ = P , where X = 1, 2, ...) or (2) to 

® ® X X-1 

use the value obtained at the end of the previous pass (i. e. , P = P^ , where 

X = 1, 2, ...). 


Three types of iteration capabilitj' are provided in OABIAS, In the first, the 

complete block Y is run through the processor in each pass. For this method 

of iteration, the user has the option of employing either of the two P^ updating 

schemes noted in the abo\'e paragraph. The number of iterations are selected 

a priori by the user; the program has not been given the capability of using the 
/^X 

Xj convergence as a criteria for automatically ending the operation. 

In the second typo of iteration, the iterations are performed on subsets 
of y . The Y^'s are composed of the observations between state vector up- 
dates. Thus, if X is being updated every seven frames, Y^^ will contain the 
ob .arvations obkuned from framas 1 to 7. The program will continue to re- 
process Y^^ , until the estimalc converges or until the specified limit 
on the number of passes is reached. It then will move on to set Y. which is 
composetl of the observations obtained from frames 8 to 14. For this mode of 
iteration, OAblAS has not been given the P^ - P^ covariance matrix resetting 
capability; it employs only the P^ ; method. 
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3.4.5 Siny:lc Ot)scrv;ition Ituralion 

In nddition to the two iteration techniques discussed above, OABtAS has a third 
iterative method called the linearity fix (Ilcference 11). Although eonccptually 
and mathematically similar to the other two iteration methods, the li'iearity 
fix should be distinguished from them. With this thii'd teehnique, the program 
iterates the scalar obscrv'ations one at a time. That is, any single obsein’ation, 
e.g, , observation j , is re-run thi'ough the filter (X = 1, 2, , ..) until the 
state vector estimate X. converges to a constant value or until the Uoer- 
specified lim... on the number of passes is reached. The program then moves 
to observation j + 1 , etc. After these operations on observation j have been 
completed, it is never necessary to recall observation j for further proc- 
essing. As a result, the method provides the advantages of block iteration 
with less degx'adation to the capability for real-time operation. The algorithm 
employed in OABIAS for single observation iteration was obtained directly from 
Reference 1. 


' The mathematics of the single observation iteration method can be delineated 

using the basic recursive estimator equations (Equation 3-3) as a starting point. 

Assume that the processing of observation j - 1 has been completed to yield 

Xj ^ and Pj ^ . Observation j is to be processed next. Let superscript X 

signify the Xth pass (of observation j) through the filter. In pass X , the single 

/\X-i 

observation iteration algorilnm uses the state vector estimate X. from 
pass X - 1 as the reference vector. (In other words, the X of Equa- 
tion (3-3) is now X. .) Equation (3-3) now can be rewTitten to encompass 
the single observation iteration option. 





(3-14a) 
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/so /s 

-J " -i-1 


While the filter is operating on observation j , the subscript j in the above 
equations is constant; X takes on values 1, 2, 3, etc. The covariance matrix 
P is not updated during the X passes. Instead, it Is updated only after X. 
has converged or the limit on the number of passes has been reached. Letting 
n be tlic total number of passes of observation j through the filter, the P 
updating equation is 




(3-15) 
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3.5 BASIC gi:omi:tky 

This section delineates the coordinate frames, direction cosine matrices, and 
ir.iin geometric variables used in the OABIAS subsystem and/or emplo}'ed in 
the discussions in Section 3.G of the eight OABIAS obscn'ation mcdcls. Exact 
definition.^ of the 12 elements of the OABIAS state vector X also are included. 

3. 5. 1 Coordinate Frame Flow Diagram 

The Cartesian coordinate frames and the main geometric variables to be used 
in the remainder of this section at'e shown in Figure 3-1. This figure uses a 
standard technique for displaying the relationships between coordinate frames. 
The circles in the figure . . present the coordinate frames. A straight line 
between any two circles defines tlie rotation by which the coordinate frame on 
the left is transformed -to the frame i the right. For examp e, rotating 
frame spacecraft inertial (SI) about its s-axis through the angle 0 yields 
frame spacecraft reference (SC), In the interest of simplicity, only the coor- 
dinate frames which are fundamentally significant have been given names and 
symbols: the remaining frames are intermediate ones of lesser importance. 

The rotation angles shown in Figure 3-1 are to be regarded as the true angles 

of file satellite being studied rather than as OABIAS time-varying estimates of 

these angles. Thus, the angles are time-invariant except for the rotation 

angle <t> and (in RAE-B problems) the detector mounting angle y . 

s 

The squares in Figure 3-1 indicate direction cosine matrices. For example, 

A is the direction cosine matrbe which transforms vector components from 
GI resolution to SC resolution. 


-«^C 
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Similarly, ii'^nsforms vector comiioncnt” from korizon-in crossing 

prime (HI') [horizon-out crossing prime (HO*)] resolution to SC resolution. 



= B,V 


!ir 


Table 3-1, which supplements Figure 3-1, defines the directions of the most 
significant axes of the main coordinate frames. In this section, coordinate 
frame axes will normally be indicated by x , y , and z with appropriate sub- 
scripts. WTicn it is necessary to specify unit vectors aloi^ coordinate frame 
axes, those shown in Table 3-1 will be used. Unit vectors along axes of par- 
ticular significance, such as the one along the satellite's spin axis, have been 
given special symbols as shown in the table. 

3.5.2 Coordinate Frames Gl* and Gl 

Frame Gl' is the conventional geocentric inercal frame defined with respect 
to the celestial equator and poles. Let O' and 6* (not shown in Figure 3-1) 
be the conventional right ascension and declination of the satellite's spin axis 
relative to frame Gl'. When 6' is close to 0 deg^e or ±90 degrees, the co- 
ordinates which O.A.BIAS uses to specify satellite attitude encounter disconti- 
nuities (see Section 3. 5. 8). OABIAS avoids these difiiculties by automatically 
performing a 45-degree rotation of frai.ie Gl' ic frame Gl in runs where 6' is 
close to (i. c. , within a user-spccificd tolerance of) 0 degree or ±90 degrees. 
Frame Gl in Figure 3-1 is the geocentric inertial frame, rotated or not, which 
is used in the OABIAS calculations. The 45-degree rotation entails the trans- 

/N 

formation of the unit Sun vector U and the orbit radius vector R from frame 
Gl' to frame Cl. 

The 45-dcgrc( rotation of frame Gl' normally is taken about the x -axis, 

Gl 

However, rotation about tliis axis produces an insufficient change in declination 
when 6' is close to 0 degree and ct' is close to 0 degree or ISO degrees. 
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Tabic 3-1, Coordinate Frames, Axes, and Unit Vectors . 


NAME OP 
COORDINATE 
FRAME 


STANDARD 

geocentric 

inertial' 

(SEE NOTE II 


rotated 

GEOCENTRIC 

inertial' 

(SEE NOTE II 



DIRECTION OF 
SIGNIFICANT 
AXES 


TOWARD VERNAL EQUINOX 


TOWARD NORTH CELESTIAL 
POLE 


DESIGNATION ALTERNATE 
OF UNIT DESIGNATION 

VECTORS ALONG OF SIGNIFICANT 
AXES UNIT VECTORS 



SPACECRAFT 

inertial' 

SI 

*S1 

''si 

*SI 

SEE NOTE 2 
SEE NOTE 2 

ALONG SATELLITE SPIN AXIS 

SPACECRAFT 

REFERENCE^ 

SC 

ml 

1^9 

SEE NOTE ? 

SEE NOTE 2 

ALONG SATELLITE SPIN AXIS 

SUN SENSOR^ 

SS 

“SS 

Vss 

*ss 

perpendicular TO SUN SENSOR 
FACE 

SEE NOTE 2 
SEE NOTE 2 

HORIZON IN 
CROSSING^ 

HI 

“hi 

Vhi 

»hi 

SEE NOTE 2 
SEE NOTE 2 

ALONG HORIZON IN CROSSlNG LOS 

HORiZON-OUT 

CROSSING^ 

HO 

“mo 

'^ho 

*MO 

SEE NOTE 2 
SEE NOTE 2 

ALONG HORIZON OUT CROSSING 
LOS 



'fixed in INERTIAL SPACE 

^FIXED IN THE BODY OF THE SPACECRAFT 

^FIXEO IN THE BODY OF THE SPACECRAFT FOR NORMAL HORIZON DETECTORS; ROTATES RELATIVE TO THE 

spacecraft body on rae b 


NOTE 1: THE DISTINCTION BETWEEN tRAMCS GV AND Gl IS DISCUSSED IN SECTION 3 5.?. 

NOTE 2. the OmCCTtQNS OF THESE AXES FOLLOW FHOM THE BASIC COORDINATE FRAME GEOMETRY 
SUMMARIZED IN FIGURE 4 I AND DISCUSSED IN TEXT. 
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Tlius, whenever a rotation is necessary and ex* is witliin ±15 degrees of 

0 degree or ISO degrees, the 45-degree rotation is taken about y , ratlior 

Cl 

than X . Figure 3-2 summai;izes the Cl* to GI rotation criteria. It should 
GI 

be noted that the transformation from GI’ to GI coordinates affects only tlic 
internal worlungs of OADIAS; it docs not affect the inputs provided by the user 
or the outputs generated by the program. 

3.5.3 Coordinate Fram es SI and SC and State Vector Elements x fiJ) and 

Frame SI is Hxed in inertial space in any single OABIAS run. Frame SC is 
file basic body-fL\ed reference frame of the spacecraft. Figure 3-3 shows the 
relation among frames GI, SI, and SC. The z - and z -axes both lie along 
the unit spin vector S of die spacecraft. Since OABIAS does not include nu- 
tation effects, frame SC is obta. om frame SI by a sit^le rotation angle 
0 which defines the instantaneous of the spacecraft in its spin cycle. 

The X - and y -axes are oriented such that rhe aximuth (measured in the 

SCr OC.* 

• plane from the +x -axis) of the center of the spacecraft’s Sun 

SC sc sc 

sensor is zero. 

OABIAS assumes that 0 is a linear function of time 

0 = 0^ + W {t - t^) (S-16) 

where 0 and to are constant, and t is the time at the start of the run. 
o o 

The value of is normally of little or no direct interest. It must be included 
in X , however, because it is an essential component in three of the eight 
OABIAS observation models. Similarly, to is an essential component of six 
of Uie eight models. 




.» NoavNnoso ivmoiuo 


3-11 


ROTLIM- USSRSPECiflED TOLERANCE 


>1 

z 

o 

*5 

z 


Q 


(A 


< 


z 

o 

QC 

< 

z 

§ 

£ 

o 


o o 

o m 
z < < 
O o o 

L in 

5 » 

w til uj 
5 I- K 
ff < < 

tt H H 

0 o o 

z cc a: 

1 I I 



^ m CO 






OF POOR QUAL-mr 




3.5.4 Sun Sensor Gooinclrv and State Vector Klcnumts x (A/3) and x (c) 

8 ~^ 10 "* — 

Figure 3-4 shows Uie Sun sensor geometry at a Sun sighting time. (3^^ on this 

figure is the angle measured bj'-tlie Sun sensor, and jS is the true Sun angle; 

y\ ^ _ 

i.e. , tlie. angle between U and S. c and Ap are bias errors which cause 

the variation between )3 and )3 . 

M 

This document employs a coordinate frame Sun sensor (SS) whose axes are 
attached to the Sun sensor. The x - and z -axes form the sensor's reference 

OO 

slit plane. Idcall}', the Sun sensor is mounted such that the spacecraft's spin 
vector S' lies in Uiis reference slit plane. The sensor sights the Sun and 
measures the angle /3 at those instants when (1) the spinning Xcc“*oo 

Ai Oo oO 

/N 

crosses the unit Sun vector U and (2) Uie angle between U and x is within 

So 

the rat^e limits of the sensor. If the Sun does lie within the sensor's FOV, it 
is sighted once per spin cycle. 

OABIAS models tlie Sun sensor biases as actual or effective alignment errors 
C and A^3 . € and A/3 are included in the state vector X as e'ements 
X - and X , respectively. C is an alignment tilt of tlie reference slit plane, 
and AjS is a rotation of the instrument in this plane. A/3 is, in effect, an 
alignment error of die main slit plane of the ^nsor. 

3. 5. 5 Horizon Detector Geomotn' and State Vector Elements x (Ay) , 

Ki) ■ Ph) • ’•u 

Parameters which are related to horizon crossings — eidier ccntral-bcd}' in or 
central-body out — are signified by Uie subscript II . Indices I and O are 
used, as necessary, to indicate iu-crossings and out-crossings, respcctivclj'. 

The horizon scanner representation used in OABIAS employs separate coordi- 
nate frames for in-crossing and out-crossing conditions. In Figure 3-1, 
horizon-in crossing (HI) and horizon-out crossing (IIO) arc the main coordinate 
frames of the horizon scanner. HI' and HO' are Uie intermediate frames and 
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are clesigaa(ccl only because later sections utilize the matrices and 
which transform vector resolution from 111' and llO' back to SC coordinates. 

The geometry of a hori7.on-in crossing is shown in Figure 3-5. and 
designate unit vectors along tlie and respectively* A horizon- 

y\ 

in crossing occurs when crosses the s’:y-to-ccntral-body boundary. A 
horizon-out crossing occurs when crosses the central-bod 3 '-to-sl^ bound- 
ary. 

The orientation of frame III relative to the satellite body-fixed reference 

frame SC is specified by angles ^ € Ay , and y . The orientation of 

n n o 

HO relative to SC is specified by 0^ , € , Ay , and y . y is the nominal 

mounting angle of tlic horizon detector's 'line of sight relative to the satellite's 
spin axis. For normal horizon detectors y is constant. Thus, frames HI 

O 

and HO and the intermediate frames HI' and HO' are invariant relative to 
frame SC. For the PAF-B PAS, however, y is stepped in 0. 7-degree in- 

O 

crements. Hence, in PAS problems, frames HI' and HO' are fixed relative 
to frame SC, but HI and HO are time-varying relative to SC. 

OABIAS models the horizon scanner bias errors through the rotation aisles 

€ , Ay , 0 !, , and 0^ . These are included in X as elements x , x , 

H 11 11 11 4 

X , and X , respectively, f and Ay are regarded as alignment errors 
5 6 H 

of the unit and arc considered to be identical for both in-crossii^s and out- 

crossings. Cjj , however, is used only witli PAS. As can be seen in Figure 3-5, 

when Ay , y , and y are constant (as they are wiUi normal horizon detec- 
s 

tors), a non-zero affects the horizon triggering times in the same way as 

a bias A0„ on the azimuth angle 0 , Hence, it is not possible for OABIAS 

to distinguish between and A0jj , and onl}' one of tlie pair is observable 

from the available horizon crossing time data. Wlicn running OABIAS, tlie 

usual teclmique for circumventing this difficult}’ is to constrain €„ to zero. 

” I 

In PAS runs, however, e , in principle, is distinguishable from A0„ and 

11 H 
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NOTE: THE RELATION BETWEEN THE SC AND HI FRAMES CAN BE SEEN MORE CLEARLY IN 

FIGURE 4-1. 


Figure 3-5. Horizon nolcclor Goonncir}" at 
Central-13ocly-lii Crossing 
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A0 , bc-causo of the lime variation of y ; hence C., is not constrained to 
tl s II 

zero in KAK-B problems. 


anv^ 0^ arc azimuth angles of tlic horizon detector and ca.n be regarded 
as the sum of tlie nominal azimuth angle 0,,vTn\T smaller perturbation 
angles A0^j and A0^j ; i. e. , 


0^ = 0 + A0* 

0^ = 0 + A0® 


The biases A0l^ and A0^ result not only from physical instrument misalign- 
meat, but also from improperly predicted electronic phenomena associated 
with horizon detector triggering. The latter tj'pe of bias is not necessarily 
identical for Earth-in and Earth-out crossings, and for this reason OABIAS 
uses distinct azimuth angles 0j^ and 0^ for in-crossings and out-crossings. 

‘ s indicated above, horizon detector azimuth bias errors A0!^ and A0^ 

H H 

are included in X . However, an analogous azimuth bias A0 of the Sun 

s 

sensor is not included. Instead, the Xg^-axis has been defined in such a man- 
ner that A0 = 0 , As shown in Figure 3-4, flie x -axis employed internally 

3 oV./ 

In OABIAS lies along the intersection of the Sun sensor reference slit (Xj,.,- 

Do 

z ) plane with the plane perpendicular to the spin vector S , The OABIAS 

DO 

state vector has been set up in this mannr ' because it is not possible to deter- 
mine absolute azimuth misaligiunents (relative to an arbitrary body-fixed 


reference fi-ame with an axis along S ) of both the Sun sensor and the horizon 
detector from the available Sun sighting and horizon crossing time data. 
Instead, only the azimuth misalignments of one instrument relative to the oUter 
can be determined. The OABIAS solution has been to constrain the azimuth 
alignment angle of the Sun sensor to zero. 
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3.5.6 State Vector Klcnicat (AP) 

Using horizon crossing time data in attl*udc determination computations re- 
quires that the values of tlic angular radius p of the central body on the unit 
celestial sphere as seen from the satellite be known. (The technique used in 
OASYS/OABIAS to compute p is discussed in Section 3.7. ) State vector 

element x is a constant angular bias AP on the computed values P of P . 
# c 

The assumption that a p^ might have a constant bias is realistic because a 
constant bias on the triggering threshold of the horizon scanner is the main 
error source which generates a AP . Figure 3-6 demonstrates this phenom- 
ena. The attitude determination algorithms assume that the sensor triggers 
whenever the center point of its FOV crosses the boundary of the central body 
disk. Figure 2-6 shows the result when the triggering occurs significantly 
early on the in-crossing and late on the out-crossing. The figure indicates 
that this triggering time error increases the apparent p of the Earth by a 
constant value, i. e. , AP , which is independent of the path of the sensor 
across the Earth's disk. This, of course, is a simplified view of the trigger- 
ing phenomena and is based on a number of assumptions and approximations 
which include the following: 

1. The FOV is circular 

2. The oblatcaess of the central body is negligible 

3. The change in tlie true angular radius p of the central body 
during file rpn is negligible 

4. The sensor triggers, on the average, when a given fraction dA 
of its FOV is illuminated 

5. dA is identical for boUi in-crossings and out-crossings^ 

6. dA docs not change significantly during the run 

^Tlus is usuallj* the most important assumption. 
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Figure 3-6. Generation of Earth Radius Bias Ap 
by Sensor Triggering Level Bias 
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From the above paragraph, it should be evident that AO is la* ^elj' a result of 
sensor behavior rather tlian of an actual uncertainty in the true Farth radius. 


3.5.7 State Vector Klemcnt (At) 

State vector element x is a time bias At in the location of the satellite in 
its orbit. It represents either a lag or a lead of tlie true in-track location of 
the satellite in its orbit, relative to either the location predicted by the OAP’AS 
orbit geaerator or tlie location specified by the orbit data supplied to OABIAS. 
The use of such a bias is convenient {' r approximating simuL'icd orbit data 
errors which result from aerodj'namic drag near the perigees of highly ellip- 
tical orbits. At has been included in X mainly to h?* *'’' > tills type of problem. 


It should be emphasized that At is not an error in the telen. ed sensor 
event times t , t , or t . At enters the OABIAS n'athematlcs only 

S tit xiCl 

dirough the orbit radius vector R 

3.5.8 Transformation Matrix A and State Vector Elements x (s , ) and 

The OABIAS algorithms require computation of the transformation matrix A 
shown in Figure 3-1. This section discusses thi • computation and includes the 
parameters s and s. which OABIAS uses as the x and x co.v.ponenls 
of 2^ to specify spin vector attitude. 

Using Figure 3-1 and elementarj' rotation matrix techniques, A may be deter- 
mined as a function of the angles O' , 6 , and 0 . The result is 


A = 


c0s6ca - s0sa 
-s0s6ca ~ c0sof 
c6ca 


c0s6sa + 60c a -c0c6 
-s0s6sO' + c0ca s0c6 
efsO! sfi 


1 - 

where s - sine and c s cosine. 


(3-17) 
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A can also be expressed in the following form; 


A = 


^.GI.T 

1 

^,GI.T 


^GI.T 


(3-18) 


where 'e* , , and “S' are unit vectors along the x-, y-, and z-axes of 

i. ^ «5 

frame SC, and superscript GI signifies resolution along the axes of frame GI, 
Superscript T signifier Ui.tt the elements are 1x3 row vectors rather than 
3X1 column vect'' lies along Oie spin axis of the satellite. Thus 


Let the scalars f , s. , and s signify the components of S , Then 

1 ^ «5 


Equations (3-V^l and (3-18) show that 



= c6ta 


(3— 19.x^ 

8^ “ efisa 

m 

(3-19b) 

S3’s6 = i/-s^-s| 


(3-19C) 


The plus sign is used with the radical sign in Equation (3-19c) when 6^0, 
and tlic minus sign is used when 6^0., 

In place of the usual or and 6 angles, OABIAS user s and s as the spin 

vr 'tor altitude coordi ates x and x of X . 'fherefore it is necessary 

1 . ^ 

to rewrite Equation (3-17) to yield A as a function of , and 0 . 
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This can be done, element by element, usings Equations (3-17) and (3-19). The 
result is 



Equations (3-20a) and (3-20b) are the ones actually used in OABIAS to compute 
the elc 'enls a^^ of A . 0 is computed using Equation (3-16). 

Eq>iaiions (3-19a) and (3-19b) show that s and s approach zero as 6 ap- 
proaches ±9C degrees. Equations (3-20a) and (3-20b) show that the elements 
a , a , a. , a^^ of A then approach a singular condition (zero divided 
by zero). The net result is tliat Uie use of coordinates s and s_ does not 
eliminaia the well known singularity at 6 = 90 degrees which is encountered 
with a and 6 coordirates. The ; irpose of the 45-<’ijgree coordiiiafe rotation 

when the user-supplied declination estimate O' is close tc *90 degrees is to 

o 

avoid lilts d;;!icully. 

s ■ ■ ' 
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As noted previously, the correct sijjn of the square root term in Equation (3~20b) 
is positi’.e when 6 > 0 and negative when 6 < 0 . OABIAS determines the 
correct sign of s at the start.of each run and maintains this sign throughout 
the run. Thus, erroneous computations would result if the filter's transient 
response were such that the sign for s had to change during the run. 

•5 

OABIAS prevents f is from occurring by performing tlie aforementioned 
45-degree rotation w. enever the user-supplied declination angle estimate 6^ 
is close to zero. The proper sign of s. is determined from the sign of the 
initial declination angle 6^ of the rotated coordinate frame GT. As a result 
of the rotation, s is sufficiently large that it wi'l not pass through zero during 
the run; hence the potential change difficulty is avoided. 

3,5.9 Transformation Matrices and 

The derivation of the equations for the direction cosines matrices Bj and 
of Figure 3-1 requires onlj* straightforward rotation matrix techniques. Be- 
cause the form of the equations for B^ and B^ is identical, only a single 

matrix, to be designated as B , is shown. The azimuth angle is s.ginfied 

il 

by 0 . The result is 

11 




'*b 

b .. 

b 

11 

12 

13 

b 

b 

b 

21 

22 

23 

b 

b 

b „ 

L 31 

32 

33 J 


■ \ \ 


cAyc,jjj - sAysCjjS0jj ^ 

cAys0jj + sAysC^^cf - cAysf^jCd^j 


(3-21) 


-sAyce 


II 


S( 


H 


cAyCf 


H 


where T ( ) signilics llic matrix for Oio rotation about the a.xis within the 
parentheses and superscript T signiftes matrix transpose. • 
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3.6 OnSKltVATION MODELS 

Th» } section discusses each of the eight OABIAS observation models. Also 

included are the derivations of the y and y algoritlims for each of tlie 

c 

models. First, however, the pertinent features of tlie general observation 
model equation used by OABIAS will be recapitulated. 

The OABIAS observation model equations are in the form 

y(ffi) = >’ (IS* S + V (3-22) 


where y * the real obsen^ation 

y = the model observation 
c 

m s the 4X1 veetor of the basic measurements for the frame being 
processed 

V = the net error 'hie to the error in m and to modeling approxima- 
tions and inaccuracy 

ffi’’ = l“i- "“2- “3- “4^ ' ‘hi- ‘ho’ 

V 

OABIAS processes the telemetry frames individually. Assuming no Sun angle 
smoothing, no observation model utilizes measurements m^ ; com more than 
one telemetry trame in any single processing operation. For this reason, in 
the remainder of this section, single frame measurement vectors m will be 
referred to rather than tlie composite block measurement vector M used 
earlier in this report. Also, the form oi Liquation (3-22) will be maintained 
throughout Uiis section with the real observation on the left and the model 
observation on the right. 

For each of the eight models there is an otiscrvation equation of th-* general 
form shown above. In Models 1, 5, 6, and 8, y is calculated using one or 
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more elements of m . In Model 7, y is the telemetered horizon detector 

angle. In Models 2, 3, and 4, y is zero bj’ definition, y is computed in each 

c 

model from an algebraic expression involving one or more elements of X 

and, directlj' or indirectly, one or more elements of m . In Models 2, 3, 

and 4, titc m dependency enters solely through y . 

c 

When running OABIAS, the user selects tl\e elements of X to be used and tlie 
models to be employed. These must be compatible with each other and with 
the attitude determination problem being pursued. Because knowledge of the 
elements of X which appear in each of the models is important, they v/ill be 
noted explicitly in the discussion of each model. (See Table 3-2 for a sum- 
mary of these elements. ) 

The OABIAS recursive estimation algorithm also requires computations of 
partial derivatives o'* the form 


hy 

= 1 to 12 

• 

The 12 g^*s are comprised in the G vector used elsewhere in the report. 
Because there are 8 models and 12 state vector elements, there are 96 partial 
derivatives including some which are zero. (Becau -.o of the number and com- 
plexi^' of the partial derivative equations, they are contained in Appendix A.) 

3. 6. 1 Model 1 — Sun Angle IVIodcl 

The angle which is measured by the Sun sensor is the observable for 
Model 1 and is clearly the most significant clement of m used by Model 1. 
Model 1 is similar to the Sun ai^lc models used in OASYS/GCONtUS and 
OASYS/GRKCRS except Model 1 includes the Sun sensor misalignment angles 
€ and A/3 . 
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The basic geometry at a Sun sighting is shown in Figure 3-4, The features 
of this figure essential to Model 1 are repeated in Figure 3-7. 

The Sun angle observation equation is obtained by applying the law of cosines 
of sides to the spncrical triangle shov/n in Figure 3-7. Thus, 

~ ^ ^ (3-24) 

A suitable equation for cos P now must be derived. This equation should 

include elements of X but not the measured Sun angle . The derivation 

- ^ 

starts by expressing cos p as the dot product of U and S 


cos j8 = 0 • § (3-25) 

Let U and S be resolved aloi^ the axes of frame GI and let their components 

be designated as u , u , u and s , s , s_ , respectively. In OABIAS, 

X ^ 3 X «5 

an ephemeris file or orbit generator provides flie unit Sun vector components 

u' , u* , u' (at Sun sighting time t ) along flie axes of frame GI’. These 
X M u s 

are identical to u , u^ , u in runs in which the 45-degree frame GI rotation 

is not made. In runs in which this rotation is i.. ^de, the OABIAS transforms 

u' , u* , u' into u , u„ , u using standard rotation equation methods. 

X A 3 X 3 3 

Hence u , u , and u are not functions of X . s* is computed using 
X M 3 3 

Equation (3 -19c) given previously. 

In sumniaiy, the Model 1 equations are 


y - (3-20a) 

y = + arc cos C(u s + u s + u s ) sec e] (3-2Gb) 

C 1 1 A 2 o 3 
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Figure 3-7. Geometry for Model 1 — Sun Angle Model 



where 
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Sg = ±/l - Sj - Sg (3-26C) 

• ^ 

The correct siga in Equation (3-26c) is determined at the start of each OABIAS 
nm. 

Equation (3-26) shows that Model 1 contains the following X elements: 

s^, Sg, and e 

Model 1 is the only model which contains ^ , and it is the onl;. model that uses 
the Sun angle measurement . Hence, Model 1 normally is included in all 
OABIAS runs. 

3.6.2 Model 2 — Sun Sighting Time ^lodel 

Model 2 employs only the Sun time measurement t^ . The essential features 
of the geometry are shown in Figure 3-4. In this figure, the Sun sensor refer- 
ence plane x _-z rotates in inertial space due to the spin of the satellite. 

Do Do 

y\ 

The Sun is sighted when this plane crosses the unit Sun vector, U . 

N is the unit vector along the Sun sensor axis y_-, which is perpendicular 
s OD 

/S /V. 

to X and z . N is perpendicular to U at Sun sighting times t . 

So oo S S 

Therefore, 


l/N 


(3-27) 


The obscr^•ablc for Model 2 is U • N = cost , where T is the ancle between 

8 

U and N as sliown in Figure 3-4. The value of cost varies as a function 
s 
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of the angle 0 of the satellite in its spin cycle. Assuming the orientation of 

U and S to be constant, cost Is a periodic function of time; its fundamental 

period is the satellite's spin period 2ir/o) . Cost is obsei*ved only at those 

times t when the Sun is sighted. However, cost is zero at t = t as was 
S- s 

indicated in Eejuation (3-27). Therefore, the real observation y and model 
obser\'ation y for Model 2 are 


(3-28a) 
(3-28b) 

where y is evaluated at the measured Sun sighting time t . 
c s 

It should be apparent that there is a fundamental difference between Models 1 
and 2. In Model 1, the significant measured parameter entered into tlie 
mathematics tlirough the real observation y . Thus, y was not known per- 
fectly. In Model 2, the measured parameter t enters the mathematics 

s 

solely through the model observation y rather than through y . With 

c 

Model 2, y is zero by definition. Hence the name "obse* vation" as applied 

to y in Model 2, in a sense, is a misnomer. However, this difference be- 

tween the methods of Models 1 and 2 does not affect the performance of the 

system as a least-squares recursive estimator; in botli cases the residuals 

y - y have equivalent properties, 
c 

A suitable equation for y containing t and elements of X now must be 

c s ~ 

derived. Using frame GI resolution. Equation (3-2 8b) becomes 

= "l "l * "2 "2 * “3 "3 


y = 0 


/V /\ 

y = U • N 


s 
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As was noted in the Model 1 discussion, the components u , u , u of 

are provided by cphcmcris data. Hence, the remaining problem is to develop 

/\GI 

equations for the components n* , n , n. of N . 

1^3 s 

Because N lies along y , 

S Oo 


j.j 

The Uj , Og , equations now can be developed with the aid of Figure 3-1 
and standard coordinate frame transformation methods. The result is 


Siss 


S 


^^3 


T T T 

= A T^^) T^(y) 




(3-31) 


where superscript T signifies a matrix transpose and 






cA/8 0 -sA8 
0 10 
sAj8 0 cAj8 

10 0 
0 C€ sc 
0 -sc cc 


(3-32a) 


(3-32b) 


A is the direction cosine matrix defined in Figure 3-1 and discussed in Sec- 
tion 3. 5. S. It is computed in OABIAS by Equation (3-20). 
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The final equation for n , n , n is obtained by substitutinj; Eouations (3-32a) 

i. ^ O 

and (3-32b) into Equation (3-31) and performing the multi plieations. By indi- 

- • - 

eating tlic elements of A by lower case letters, as shown in Equation (3-20), 
the result is 




cc + a. 


31 


cc + a, 


32 


cc + a, 


33 



(3-33) 


A A docs not appear in the result. 

The final equations for Model 2 are Equations (3-28a), (3-29), and (3-33), 

Equations (3-33) and (3-20) show that Model 2 includes the following elements 

of X : s, y s , , w , and c . When computing 0 and, hence, A for 

Model 2, the measured Sun sighting time t is used. This is the mechanism 

s 

through which m enters the Model 2 mathematics. 

3,6,3 Model 3 — Nadir Vector Projection Mode l 

The central body horizon-in crossing time t^^ or the horizon-out crossing 

time t are the elements of m which are used in Model 3, When tlie user 
HO ^ 

requests Model 3, separate calculations using first t and then t are 
performed. In the interest of simplicity, the present discussion and notation 
will not distinguish explicitly between the liorizon-ln crossing and horizon-out 
crossing cases. 

Figure 3-S summarizes Uic esse. .dal geometry fo” Model 3. ui tliis figure, 

^ is the unit vector along the lino of siglit of the horizon scanner and r.)l:!tcs 
in inertial space due to the satellite’s rotation. Figure 3-8 shows the locus 

y\ 

traced out on the unit celestial .sphere by Uic tip of L . -11 is '.’»e unit nadir 
vector of the central body; it moves relatively slowly on the celestial sphere 




Figure 3-8. Geometry for Model 3— Nadir Vector 
Projection Model 
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due to the in-track motion of the satellite around its orbit. P is the apparent 
angular radius of the central body on the unit celestial sphere as seen fron; the 
satellite. ■ 

/S 

The observable for Model 3 is the parameter cos p + R • L . This parameter 
varies c 3 'clically at the orbit pex'iod 2ir/w . There is also a slow variation due 
to the motion of R on the celestial sphere and (in elliptic orbits) the chatting 
magnitude of P , The measurements of this observable are contained in the 
horizon crossing times t , A check of the geometry in Figure 3-8 will show 
that the observable is zero at these times. In this respect. Model 3 is analo- 
gous to Model 2, The equations for y and y , therefore, are 

c 

y = 0 (3-34a) 

= cos P + ^ • i) (3-34b) 

where y is computed at the horizon crossing times t . 

C IT 

OABIAS uses the following equation to obtain cos P for Uie y computation: 

c 

cos P = cos (p + AP) = cos P cos AP - sin P sib AP (3-35) 


where p is tlie computed ''aluc of P obtained using orbit data. (.-Jee Sec- 
c 

tion 4.G for a discussion oi Uiis computation. ) AP , which is element of 

X , is a bias on P . 

** c 

To derive the equation for R • L , frame GI resolution is used. Let R be 
the distance vector from the Earth's ccnlci to Uic satellite. Then 


gOI 



(3-36) 
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Letting' the components of H ■* be designated as R , R , R , 

X M «5 

(u| = /Rj + Rj * Rj (3-37) 

The Rj^ , It^ , R^ values are generated for the OABIAS reeursive estimator 
via an ephemeris file or orbit generator; these data are transformed as nec- 
essaiy in those applications where the 45-dcgree GI' to GI traasf».^*.*ation is 
made. 

Let the components of L^* be designated as I, , l„_. . Use of 

IGl oGl 

Eijiatiops (3-a5) and (3-36) now enables the y expression in Equation (3-34u) 

0 

to i?:"’ ■ ‘ten. as follows: 

y « cos p cos Ap - sin p sin Ap 
c c • c 

(3-38) 

jW[ !^1 'iGI ^2 *2GI * ^3 ^3GI) 

Suitable equations for the components 1, __ , , l„__ of are derived 

IGI 2G1 3G1 

using Figure (3-1) and the usual rotation matrix techniques. The intermediate 
result is 


, Ty?H ^ , T 
L =T^(y) L =Ty(y) 



-^GI 



- A 


T^SC 

JLi 
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(3-39b) 

(3-39c) 
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f» •'Si 
\^icre 


TyOO 


% ® -“Vs 


By^ 0 


C7 


(3-40) 


Th" aecoad part of Equation (3-39a) follows from the fact that £ lies along 

the Zjj-axis. The equatioas for the elements and of matrices A 

and B are Ecpiations (3-20) and (3-21). The final scalar equations for I 

% 

1_ , l„_. are obtained by performing the matrix multiplications of Equa- 

2G1 oGI 

tions (3-38), (3-39), and (3-40). The result is 


IGI ’ 


(*lSCi 




^ *^33 % * ** 


3SC- 


(3-4 la) 


IGl 


t°'=n 


2G1 


3G1 



1 ^ 

+ a^ 

1 ^ 

+ a 

1 

11 

ISC 

21 

2SC 

31 

3SC 


1 

+ fl 

1 

+ a 

1 

12 

ISC 

22 

2SC 

32 

3SC 


1 ^ 

+ 9 

1 

+ a 

1 _ 

13 

ISC 

23 

2SC 

33 

3SC 


(3-4 lb) 


The final equations for Model 3 are Equations (3-34), (3-37), (3-38), and 
(3-41). Equations (3-lG), (3-20), and (3-21) also are needed for computation 
. of the elements of the A and B matrices. 

In using Equation (3-lG) to cominUc the angle 0 which appears in A , the 
time t is the measured horizon crossing t 

* 0 I 1 o “ 1 1.. “ t ! (3*42) 

model 3 o I II 01 ' 
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The elements of It , however, are evaluated at t^^ + At where At is 
element x of X. Those are the two mechanisms through which m enters 
into Model 3. 


A check of the equations will show that Model 3 includes the following elements 
of X : * 

3.6.4 Model 4 — Horizon Crossing Time Model 


As with Model 3, Model 4 employs the central body horizon crossing times t„. 

til 

or t„^ . Likewise, separate calculations using t„, and t„_ are made 
HU HI HU 

with Model 4. Once again, the present discussion and notation will not dis- 
tinguish explicitly between horizon-in crossings and horizon-out crossii^s. 


Figure 3-9 shows the essential geometry for Model 4. N is a unit vector 

H ^ 

along the j'-axis of frame II and the horizon detector line of sight vector L 
lies along the z-axis of frame H . Hence, N and L are perpendicular to 

^ /s 

one another and N • L = 0 . Both vectors rotate in inertial space due to the 
H 

spin of the satellite. 


Figure 3-9 shows the locus swept out on the unit celestial sphere by the tip of 
. 6 is a unit vector whose tip is at the intersection of the perimeter of the 
central bod 5 ' disk and the locus of L . D moves (relatively slowly) due to the 
motion of the satellite in its orbit which produces motion of the central body 
disk on Uic unit celestial splu 'c. It is important to realize that the D motion 

/S 

is not a function of tiic satellite's spin rate w . 

The observable for Model 4 is cos C where C is the angle between and 
D . Thus cos C = Njj • D . Cos C ”" ries as a hmetion of tlie ai^le of the 

satellite in its spin cycle. The motion is near-periodic witli the fundamental 
licriod being tlic satellite's spin peri«»(l 2it/oi . There also is a slow variation 
due to the motion of D . Cos C is observed at the horizon crossing times t^^ 
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SATELLITE SPIN AXIS 



Figure 3-9. Geometry for Model 4 — Horizon 
Crossing Time Model 
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when L and D coincide. Because L is ortliogonal to , cos C - 0 at 
t = tjj . Therofoce, the y and equations are 


y = 0 , (3*>43a) 

y =N„* D (3-43b) 

•'c H • 


where y is evaluated at the horizon crossing times t . 

C H 

Resolution on frame GI is employed in the y computation. Thus, 

c 


aGI aGI 


(3-44) 


Referrit^ to Figure 3-1, 


(2)^T (x)^T 

H 0JJ fjj Av H 


.TaH* 


(3-45a) 




(3-45b) 


where 0„(z) , c (x) , and Ay(y) are the usual coordinate frame transforma- 
li ii 

tion matrices. From the definition of frame H' and the vector N , 

H 



(3-4G) 


J ( 


•* 
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The sc.alar equations for are obtained by substituting Equation (3-46) ^to 
Equation (3-45) and performing the niatrl ’ multiplications. The result is 


«h| 

C0R (3-47a) 


1 

yvGI 
N = 
H 

( "hIGI \ 1 *^11 "hISC hi “h2SC ^31 "h3SC^ 

I “h2GI / " S *^12 “hISC *^22 "h2SC ®32 "h3SC . 
(”h3GI ) V13 ”h2SC ®23 "h2SC * ®33 “hSSCJ 


^Gi 

The scalar equations involving D can be obtained with the aid of Figure 3-9 


>^GI ^ gGI ^ ^ 

(3-48a) 


. 6^^ = -c(p + Ao) 

C 

(3-48b) 


A 

^Cl ^Gl , 

D • D =1 

(3-48C) 



At any stage of the OABIAS computations, numerical values (or estimates) for 
all parameters in Equations (3-48), except B , are available. Hence, Equa- 
tions (3-48a), (3-48b), and (3-48c) constitute a set of scalar equations which 
can be solved for the throe components of . Because the derivation of the 

algorithm for is lengUiy and not essential for the discussion of Model 4, 

/sGI 

it will be deferred to Section 3.8, However, it should be noted that D is 
a function of the following elements of X ; s^ , s^ , Ay , AP , Cjj , and At . 
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The ccfuations for Model 4 are (3-43a), (3-44), and (3-47), plus the solution of 

GI 

Equation (3-48) for D . In addition. Equation (3-20) is required to calculate 

the elements of A , and a supplementary calculation of P iS needed. The 

— ^GI ^ 

algoritlim also requires 11 and the nominal mounting angle y as inputs. 

s 

As in Model 4, 0 (and hence A) is evaluated at the measured horizon crossing 
“^GI 

time t.. . R (and hence p ) is evaluated at t„ -f At where At is element 
H e H 

X - of X . These are the two places where m enters into the Model 4 

mathematics. A check of the equations should show that the X elements 

which appear in Model 4 are the same as in Model 3: s , s. , , Ay , 0^^ 

O , 1 2 o H 

or 0JJ , AP , w , Cjj , and At . 

3.6.5 Model 5 — Sun to Earth-In and Sm to Earth-Out Dihedral Angle Model 

Model 5 uses two elements of m : the Sun sighting time t and either the 

s 

horizon-in crossing time t^^j or flie horizon-out crossing time t^^^ . When 

the user requests Model 5, separate calculations using t^ and t^^ and then 

t and t are performed. In the interest of simplicity, the present d? jus- 
S nv/ 

Sion and notation will not distinguish explicitly between horizon-in crossings 
and horizon-out crossings. 

The geometry for Model 5 is shovra in Figure 3-10. The left portion of the 
figure shows the geometry at a Sun sighting time t^ ; the right portion shows 
the geometry at a horizon crossing tine t^^ . Since the spin rate to is con- 
stant, the total spin angle change between t and t will be tO(t - t ) . 

s H H S 

The observable for Model 5 is (t„ - t ) . Thus 

H s 

y = - t^ (3-49) 
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SATELLITE SPIN AXIS 




UNIT CELESTIAL SPHERE 


Figure 3-10. Geometry for Model 5 — Sun to Earlli-bi and 
Sun to Earth-Out Dihedral Angle h jdel 
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The basic equation for the model observation y can be developed with the 
aid of Figure 3-- 10. 


" W ^ - 0JJ - 2 irn} (3-50) 

, as shown in Figure 3-10, is the dihedral angle between the U-S plane 

and the ^-6 plane. As/w is the change in the Sun sighting time t due to 

s 

Sun sensor misalignments € and A/3 . Similarly, Ali/o} is the change in 
t„ due to horizon scanner misalignments € and A7 . 0„/w is the time 
required for the satellite to spin through the azimuth angle 0 between the 
Sun sensor and the horizon scanner. 2an , where n ^ 1 , 0 , or <1-1 , is an 
additional term which is required to make the angle inside the braces fall 
within the proper range under all conditions. 

A in Equation (3-50) is computed using the well known dihederal angle ex- 
do 

pression 


A ... *= arc tan 
d5 



5- [uxB3 

D-Cfi. g^][u. 


(8-51) 


where all ve tors are resolved on frame GI. 

As noted previously, U is provided by ephemeris data and is evaluated at 

/^Gl 

t = t . The elements of S are s. and 


, which are the state vector 

com- 


t = I . cic;iiu:ut9 vi o ui-e s aw , which ure uie siaic vei 

S 12 g. 

elements, and s wind is obtained via s = ±yi - s - S- . d'^ is 
3 3^12 

puted using the technique summarized in the Model -1 discussion and delineated 
in Section 3. 8. When computing , the vector is evaluated at 
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The term As in Equation (3-50) is computed using the two spherical triangles 
associated willi the Sun sensor in Figure 3-10, From the upper triangle and 

the law of sines, • ' 


sin E = sin c/sin $ 

(3-S2a) 

From the lower triangle and the law of sines, 


. * sin E sin (90 - /5) 

sin As ; — \ ^ 

sin A 

(3-52b) 

Combining Equations (3-52a) and (3-52b), 


. ^ sin t cos 

sin As = -■ o . — ^ 
sin p sin A 

(3-52c) 

From the lower triangle and the law of cosines of angles, 


cos A = sin c sin $ 

(3-52d) 

Comparing Equations (4-52a) and (4-32d) shows that 


A = 90 - (T 

(3-52e) 

The final cquatioci for As is obtained by substituting Equation (3- 

-52e) into 

(3-52C). 


As B arc sin (tan c ctn $) 

V 

(3-53) 

♦ 
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When using Kquntion (3-53), P is computed as follows: 


P = arc cos (U « S ) 


(3-54) 


The remaining task is to develop the equation for AH to be used in Equa- 
tion (3-50). Figure 3-10 shows that this task is geometrically identical to 
tJhe development just completed for computing As , Hence, Equation (3-53) 
can be employed directly wiUi appropriate changes In variables. The result 
is 


AH = arc sin Ctan e ctn y] 

XI 


(3-55) 


where y is calcubted Croxr. 


y = arc cos [ • 6^^] (3-56) 

In summary, tlie equations for Model 5 are Eqiiati^'ns (3-49) through (3-51) 
and (3-53) through (3-56). A check of these equations shows that Model 5 uses 
the following elements of X s ®i * ®2 * * ^H ^ ^ * 

and At . 

3.6,6 Model 6 — Earth Width 

Model 6 uses two elements of m ; the horision-in ct'ossing time t^^j and the 

horixon-out crossing time t . The observable is the difference between 

HO 

these two times. Thus 


y 




(3-57) 

r 
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Figure 3-11 shows Uic geometry for Model 6. Model 6 is mathematically 
identical to Model 5 except the Sun time and geometric variables of Model 5 
are replaced by tlio Sun time and geometric variables associated with the other 
horizon crossing. 

The basic y equation for Model 6 can be written with tlie aid of Figure 3-11. 
c 

The result is 


y 


c 



+0jj + 27rn 


■] 


(3-58) 


Equation (3-58) is analogous to Equation .(3-50) for Model 5. Equation (3-55) 

shows that AH and AH are identical in value. Hence, they cancel one 

^ I O 

another and do not appear in Equation (3-58). 0 and 0 are elements x_ 

ti n 5 

of X . 

6 

in the y^ equation is computed from Dj , ^ , and using the 
usual dihedral angle expression 



= arc tan 




(3-59) 


where all vectors are resolved on frame GI. 

/sQl 

Dj in the above equation is computed by tlie usual method (see Section 3. 8), 
with and Pc evaluated at t + At , 6^^ is computed by the same 
technique witli u ' and P evaluated at t + At . 

C il V/ 

A check of the equations shows tliat Model 6 contains the following elements 
of X ; Sj , Sg , Ay , 0[j and 0JJ , AP , O) , , and At . 



^-75 

f 




3-76 


page w 
quality 

3.6.7 Model 7— Small Tnr«;ct Model 

Model 7 was developed for Oie ti’ansluuar phase of the UAE-B mission. In 
this phase, the angular radius of the eentral body (the Moon) was sufficiently 
small that it could be approximated as a point source for attitude determination 
purposes. The PAS horizon detector angle y was stepped in a predetermined 
manner. The telemetered data included the relative times (t - t ) when the 
central body horizon was detected and the scanner angles y^ at these times. 

The geometry for Model 7, shown in Figure 3-12, is mathematically similar 

to Model 1. The observable for Model 7 is the horizon scanner angle y . 

8 

The observation equation can be derived using Figure 3-12 and the law of 
cosines of sides The result is 


OWGINAL 
OF POOR 


y *= -Ay + arc cos 


[ cos y *1 

cSi-c~J 


(3-60) 


Cos y can be computed follows; 


cosy = - 


(3-61) 


Therefore, the Model 7 equations for y and y are 

c 


y «=y. 


y *= -Ay + arc cos 
c 


I .! 


(3-62a) 


(3-62b) 


'•‘Gl 

In Utc R-\K-B application, H was computed at frame times which were 
accurate to ±0.5 spi \ period. 
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BODY 



Figure 3-12. Geometry for Model 7— Small Target Model 
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Equalioa (3-G2) shows lliat the following elcnicals of X appeared in Model 7: 
®1 * ®2 * ^11 * 

3.6.8 Model S — Sun to Earth Mid-Scan Dihedral Angle Model 

Model 8 originally was designed for the RAE-B mission and, like Model 7, 
approximated tlie central body as a point. However, Model 8 has recently 
been modified to eliminate the restriction that the central body must be small. 
The present discussion covers only the new version of Model 8. 


Model 8 uses three elements of m : the Sun sighting .ime t and the two 

s 

horizon crossing times t and t . The geometry for Model 8 is shown 

HI HO 

in Figure 3-13. The Model 8 observable is the time interval between t and 

s 

the midtime between the horizon-in crossii^ and the horizon-out crossing 


y.0.5tt^*tHo3-t, P-63) 

Model 8 is mathematically similar to Model 5 except t and t of Model 5 

/s. HO 

are replaced in Model 8 by 0. 5 tt + t 3 . Also D , D of Model 5 are 

Hi HO 1 O 

used for die Model 8 dihedral angle (A^) computation. The equation for die 

model observation y can be obtained with the aid of Figure 3-13 or by 

o 

analogy to the corresponding Model 5 result— Equation (3-50), It is 

J'c ' S hs - '’•® K * oj,! - AH * 2 ™] (3-«4) 

In Equation (3-64), 0^ and 0^^ are the horizon detector azimuth angle ele- 
ments X and X . Motlcl S thus cniplo)'S the average azlnii.'h aiwlc of the 
5 0 

horizon detector. 
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Fipjiirc 3-13. Geometry for Model 8 — Sun to Earth 
Mid-Sc:ui Dihedral Anglic Model 
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The dihedral angle A is shown in Figure 3-13. Us equation can be obtained 

dS 

by using the Model 5 result — Equation (3-51). 

In fact 



0.5 



(in) + A 


d5 



and 


= 0.5^y^ (Model 5, in) + y^ (Monel 5, out)^ 




arc tan 


s • l6 x6j] 


u- 6j-cej.§Ke. 9]| 


+ arc tan 


s- [0x6^3 


U. 6q-C6q.§]C0.§]| 


(3-65) 


where all vectors arc resolved on frame GI. U is computed at t ® t . 

6 

As and AH were defined previously (see Figure 3-10 and Equations (3-53) 
and (3-55)). For As , 


As = arc sin (tan ( ctn $) 


(3-66a) 


where /3 is computed from 


fi ~ arc cos (U^^ . §^^) 


(3-6Gb) 


'I 


; 
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- • 

AH = arn sin [tan Cjj ctn y] 


(3-67a) 


vhere y is computed from 

y = arc cos [cos cos (v^ + Av)] (3-67b) 

In summar 5 % Model 8 employs Equations (3-63) throi^h (3-67). A check of the 
equations shows tliat Model 8 uses the following elements of X : * 

. AP , to , Ay , C , and At . 
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3.7 COMPUTATION OT CUNTIIAL BODY ANGULAR RADIUS P 

c 

OABIAS uses Uie following cquntioii to compute the angular radius P of the 
central body ou the unit celestial sphere as seen from tlie satellite: 


P = P + AP 
c 


(3-68) 


where P is computed deterministically using orbit information and AP is 
c 

a bias element which is included in the state vector X . 

OABIAS provides two methods for computing P . The first assumes that the 

c 

central body is spherical; i. e. , central body oblatcness is omitted. Fig- 
ure 3-14 shows the geometry for this case. For notational simplicity, the 
magnitudes of vectors will be referred to in this discussion as they are desig- 
nated in Figure 3-14. B 3 * definition, P is the angle between the vector to the 

c 

center of the central body -R and the horizon sensor line-of-sight vector L 

at a horizon crossing. At horizon crossings, L is tangent to the central body 

surface. Hence, for a spherical central body, L is perpendicular to the 

radius vector ^ shown in the figure. By simple trigonometry, 

6 

* arc sin 

where 0 degree ^ P ^90 degrees , R is the nominal radius of the central 
** c ” e 

bedy, and h is a user-supplied correction. R is obtained from orbital data. 

The second p computat'.on method allows oblatcness to be taken into account 
c 

when the central body is the Earth. The mathematical model is an approximate 
one which uses Equation (3-69) willi R redcQned to be the geocentric radius 
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Figure 3-14. p Computation Geometry— Spherical 
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of Earth at the horizon crossing point. R is computed by the following 
equation: 


O^ilomcters) = 6378. 16 


2 ^ 

s X + c. 


sxj 


+ h 


(3-70) 


where X is the geocentric latitude and and c^ are constants which 
default to 0.0033528 and 0.0, respectively, c sX permits the user to include 
seasonal variations in the effective thickness of the atmosphere. 

In order to compute with the oblate Earth model, the location of the 
horizon crossing on the Earth's surface must be known. However, the com- 
putation of this location requires knowledge of the value of R^ . Because the 
equations of the model are hi^ly nonlinear, an iterative procedure is needed. 
The main steps employed by OABIAS are summarized below utilizii^ the nota- 
tion in Figure 3-14. 

1. The latitude X of the subsatellite point is computed as follows: 

o 

. , r ftGI' 

/s 

-R is obtained from orbital data, and because Zq^i is the 

unit vector along the z-axis of frame GI*. = [ooi]"^ . It 

should be noted that if a 45-dcgree rotation has been made, the 

vectors must be retransformed to frame GI. X serves as the 

o 

initial estimate of the latitude of the horizon crossit^ point. 

2. R is computed using Equation (3-70). The latitude X obtained 

c 

. in step 1 is used in the first pass. 

3. P is computed using Equation (3-69) with R from step 2. R is 

c e 

’ obtained from orbital data and h is supplied by -the user. 
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/SQI* 

4, The unit horizon crossing vector L is computed in subroutine 

CONES using inputs (from orbital data), and y (from 

X and y ), iuid P (from step 3). CONES returns two solutions 
s c 

corresponcing to Ear(h-in and Earth-out crossings. The correct 
one is established with the aid of subroutine PHASE. 


5. L is computed using L R cos P . 

c 

-fc-GI' 

6. The Earth radius vector to the horizon crossing point is 

computed using 


e 

where R is obtained from orbit data, L fr^m the computations 
in step 5, and L from step 4. 

7. is now normalized to yield fiO"’. 

e e 

8. The latitude X of the horizon crossing is computed using 

X,* arc sin ; 

9. The program returns to step 2 and continues in a loop until the 

change in P computed in step 3 is less than 0. 0001 r.ndians or 
c 

until the maximum number of iterations (10) is reached. 

In system testing, it has been found that convergence usually occurs in one or 
two iterations. 
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3. 8 COMPUTATION OF HORIZON CROSSING VECTOR B 

Models 4, 5, and 6 require that the following set of nonlinear algebraic equa- 
tions be solved for : 


= cos (y + Ay)’ cos c„ 
''s H 

(3-71a) 

.^GI ^ gGI _ ^ ^ 

c 

(3-7 lb) 


(3-71C) 


The purpose of this section is to discuss the algorithm employed to solve this 
set. The variables in the above equations were defined in Figure 3-10. Super- 
script GI signifies that OABIAS uses GI coordinates for the calculation. Let 
the componenU >f S , B and R be designated by lower case letters. 
Then Equation (3-71) can be rearranged into the following form: 

pi ' 2 ] pd . p“ 0's * ‘h • “s ‘*3 

h ’^ 2 ) Kl ' fo * ^0) - >^3 dj 

+ d^ = 1 
1 2 3 


(3-72a) 

(3-72b) 


Assuming that the determinant (Del) of the 2 x 2 coefficient matrix in Equa- 
tion (3-72a) is not zero, Equation (3-72a) can be rearranged into the following 
form: 


Plj _ _1_ p2 -'‘'2I 0's * '11 ■ '3 S 

(d^j " Dot [-i j Sj (-cos (P^ + AP) - Tg dj 

2x2 2x1 


(3-73) 
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where 
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Det'« s r 

X u 


®2 


(3-74) 


The condition s r - s r = 0 occurs when the two vectors formed by proj- 
X M X 

ecting S and li onto plane X - Y of frame GI are colinear. The problem 
can be circumvented (except in the unscivable degenerate case of S = by 
permuting the subscripts of the gOI _ gGI 

, and components. 

Equation (3-74) can be put in the form 


= Pi ^3 + qj (3-75a) 


where 


•k 

n S' — [s r - s r 1 
^2 Det *■ 3 1 1 3-' 

Qj = ^ tr^ cos (y^ + Ay) cos + s^ cos (p^ + AP)] 


^ [r^ cos (y^ + Ay) cos f + s^ cos (P^ + AP)] 


(3-75C) 

(3-75d) 

(3-75e) 

(3-75f) 


Substituting Equations (3-75a) and (3-75b) into (3-72b) and rearranging yields 



where 
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Solvit^ Equation (3-76), 


a^ = l^Pl^P2 




2 2 , 

a3 = qi^q2-l 


(3-77a) 

(3-77b) 

(3-77C) 


**3 " [“ ^2 *'^2 “ 


(3-78) 


Equations (3-75), (3-77), and (3-78) constitute the basic algorithm for solving 

Equation (3-71) for the components d , d , a:id d of , However, the 

12 3 

problem of resolving the sign ambiguity in Equation (3-78) still remains. 

Geometric considerations show that one sign signifies a horizon-in crossing 

and the other, a horizon-out crossing t^^^ . However, mere knowledge 

of whether a t„, condition or a t condition is being processed is not, by 
HI HO 

itself, sufficient to resolve the ambiguity. Therefore, OABIAS first computes 

^G1 

the horizon crossing vector (D ) using the plus sign In Equation (3-78). The 

/S /N 

subroutine PHASE then computes the dihedral angle from the * S plane 
to the - [-R] plane, i/ is shown for botli the in-crossing and the out- 
crossing cases in Figure 3-15. The dihedral angle computed by PHASE is 
measured In the conventional right hand sense and will be in the range from 
0 degree to 3G0 degrees. Since the angular radius of the Earth is always less 
than a right angle, 1/ always lie in the range from 0 dcgi ee to 180 degrees for 
an out-crossing. For an in-crossing, u always will be in the range from 
180 degrees to 360 degrees. Ilencc, the sign ambiguity in Equation (3-78) can 
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HOaiZON-IN CROSSING 



Figure 3-15, Resolution of Sign Ambigui^ rf 
Horizon Crossing Vector D 


3-90 


ORIGINAL PAGE Itt 
OF POOR QUiU.nY 


be resolved bj' checking whether v is within the correct range. The exact 

* 4 * 

criteria is as follows: 


Value 

Horizon 


(degrees) 

Crossing 

Correct Sign 

0 5 i; s 180 

in 


+ 

out 

+ 

180^1/ S3G0 

in 



out 

- 





A check ox the D computation equations shows that the following elements 

of X are used: s , s , Ay , € , and Ar . The element At of X also 

12 H ■* 

appears implicitly Uirough R . 
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3.9 WlilGIITING FACTOUS OF THE OBSERVATION MODELS 


Each time one of the eight obscn'ation models in OABIAS processes an obser- 
vation j , the prelim iniiiy operations include Oie caleulation of a weightijig 
factor . Each model has its own distinct equation for computing . 
These calculations use the following user-supplied inputs: 


1. a - -The square root of the variance of the error in the Sun angle 
measurements (degrees) 

2. a — The square root of the variance of the error in the Sun sighting 

s 

time measuremcats (seconds) 

3. ffjj — The s^are root of the variance of the error in the horizon 
crossing time measurements (seconds) 


4, Oy — The square root of the variance of the error in the RAE-B PAS 

angle (degrees) 


The optimum technique for establishing the observation weighting factors w^ 
v'^s discussed, in general terms, in Section 3.3.3. Ideally, each w^ should 
be the inverse of the variance of the error in the observation which is to be 
processed.^ This conclusion assumes that the errors in the observations are 
statistically uncorrelatcd. The final result of Section 3. 3.3 was Equa- 
tion (3-llc). 


With modification in notation to tailor it to the present discussion, this equa- 
tion can be written as 


CSC rt'porfs on atliludo estimation have been inconsistent in the use of the 
t:)rm "weighting factor" and the corresponding symbols W or w . The 
actual OABIAS coding uses 'Jie iuv. cs of the weighting factors of Uic pres- 
ent discussion and dcsigiiates Uiem by W( ) . 
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(3-79) 


where z. signifies the residual (y. - y .) of observation i and w* signifies 

j J cj' j 

the optimum weighting factor. The third term on the right side involves the 
assumption tliat the statistics of the error in the horizon crossing time meas- 
urements are identical for in-crossings and out-crossit^s. 


The values of the partial derivatives in Equation (3-79) depend upon the model 
being processed. For example, with Model 1 (the Sun angle model) &e basic 
measurement .8 enters z throu^ y;y = 8. Hence, Sz/b8 = 5y/S8 = 1 
and the remainir^ pairtials are zero. With Model 8 (the Sun to Earth mid-scan 
dihederal angle model), y = 0.5 Ct^j + " *-g * Hence, oz/d8 = 0 , 

b z/d t = d z/a t = 0. 5 , and a z/a t = -1 . 

Hi HU S 

A modification to Equation (3-79) is convenient for those models (2, 3, and 4) 

in which the measurement-dependency enters the residual z^ via die rotation 

angle 0 in the y computation. For these three models, the following is 
c 

adequate: 


bz 

at 


() 




at 


() 



where ( ) sigjiifics s , III , or 110 . 



Then, Kquation (3-79) becomes 
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€ 

2 22 /^^e\ 

a/w*) = 0 . + 0 . Oy +« 


(3-80) 


+ w 


\oi& ]. , ^ 

' o/ in-crossing 




© 


out-crossing 


II 


The equations for the partial derivatives dy /ht are given in Appendix A. 

c o 

OABIAS has used two methods to compute die wei^ting factors; both are avail- 
able in the current version (see T^dile 3-2). In principle, because Mediod 2 
yields wei^ting factors which are in closer agreement with the theoretical 
ideal w* , it should be superior to Method 1. This can be seen by a careful 
comparison of Table 3-2 with Fixation (3-79) or, for Models 2, 3, and 4, 
Equation (3-80). Testify also has shown that Method 2 yields better results. 

The question of resolution effects in the wei^ting factor computation now will 
be discussed. The conventional statistical filterimr theory which was implicitlj' 
usee! in the preceding development assumes that the errors in the input meas- 
urements can be modeled statistically as white noise. With digital systems, 
the resolution error due to finite word length of the sensor or processing tech- 
niques is sometimes significant with respect to this type of noise; in some 
cases, may completely dominate it. Ou attitude determination sj'stems, the 
resolution of digital Sun sensors, in particular, can be sufficiently great that 
difficult questions arc raised on how to handle or weight ‘he data. It is well 
known that errors due to resolution can be treated as white noise when tlie 
signal is sufficiently variable tliat it rarely stays in the same resolution cell 
on successive measurements. The near-constancy of tlie Sun angle seen by 
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Table 3-2. WcigUtii^ Factor Equatioa Used in OABIAS 


Model 

1 

2 


Observable 


5 

6 

7 

8 


3 

COS 


D .5^, 


H 


H s 


Sio -' 


HI 


Weigbtit>g Factor (ce ^ ) Equation 


Method I 


.-1 


.-1 




' LJjLfil 

. PS7.296/ 8J 

[k)T 

/ [‘"Z " 

flvf • 

[&)'f 


Method n 


-1 


[feldT 

l-B-JT 

[(■fefJ-l 

K>' * 

[^f] ‘ 


NOTE 1; See AppcndL'c A for tiic equations for the partial derivatives 
NOTE 2; The terms within the brackets have tlic dimensions of variance. 
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spin stabilised satellites, however, makes this assumption difficult or impos- 
sible to justify in processing tlic data from tlie digital Sun sensor. 

The question of how, in theory, the resolution phenomena can best be handled 

has not been addressed in this document. The equations in Table 3-2 assume 

that resolution is negligible and, therefore, that the errors in the Sun angle. 

Sun sitting times, and horizon crossirig ..mes can be modeled as white noises. 

The square roots of the variances of these noises are » ct,. . and a . 

P s H y 

The single exception is a. in Method I which is modeled as tn«. resolution 

p 

ceil width. 
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3. 10 SUMMARY OF SECTION 3 


The elements of the state vector X used by OABIAS arc as follows: 


Xj(Si> 


=‘2<V 


*3»o> 


x^(Ay) 

4 *\i) 


x^(AP) 

Xg(Aft 

V"* 

Xio«) 

X,J«„) 

Xj2(At) 


component of spacecraft's unit spin vector (S) 
along the x-axis of frame GI 


component of S along the y-axis of frame GI 

initial phase angle of the spacecraft in its spin 
cycle 

bias on horizon sensor mounting ai^le 

azimuth of the horizon sensor relative to the Sun 
sensor at borizon-in crossing 

azimuth of the horizon sensor relative to the Sun 
sensor at horizon-out crossing 

bias on the angular radius of the central body 

bias on the Sun ai^le measurement 

spin rate 

tilt of Sun sensor reference slit plane 
PAS horizon detector plane tilt (RAE-B) 
timing bias in spacecraft orbit data 
A summary of each of the eight OABIAS observation models follows. 


Model 1 — Sun Anglo Model 


y 


M 


y = -AS + arc cos 


IpG^t 


I 


COS c 


) 
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ft'hore P , = measured Sun angle 
M 

U = unit Sun vector 

A- 

S = unit spin axis vector 
GI = resolution on frftme GI 


Model 2 — Sun Sighting Time Model 

y»0 

^GI qGI 

-Ng 

|COS €} 

|sin c i 




a a 
21 31 


^22 %2 


S3 ^saj 


where Ng = unit vector perpendicular to &in sensor reference slit plane 

a = elements of transformation matrix A from frame GI to frame SC 

ij 

Model 3 — Nadir Vector Projection Model 

y = 0 


ftGI ^GI 

y = cos (P + AP) + R * L 
•'c c 

l” . [OOlf 


where P = computed angular radius of the central body 

^ ' 
U = unit vector from central body toward spacecraft 

L = unit vector along horizon detector line of sight 
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B 


TyO') 


11 

T 


= ti'ansformatioa matrix from frame H (Figure 3-1) to frame SC 


cos y 

0 

siny 

s 


8 

0 

1 

0 

-siny 

0 

cos y 


s sJ 


= nominal mounting angle of horizon detector relative to spin axis 


Model 4 — Horizon Crossing Time Model 


y c 0 


/nGI 


B N 

H H H 

= Cooif 

ll 




where N„ = unit vector perpendicular to PAS plane 
H 

fiGI = unit vector on central body perimeter at horizon crossing point 
Model 5 — Sun to Earth-In and Sun to Earth-Out Dihedral Angle Model 


y = -i- U - As + All - 0„ + 2irn) 
”'c .U) d5 H 

1 S.rOxn] ( 

Ajj ’ jffTir7iirrsinrr¥]t 

As = arc sin [tan C cln /5] 

= arc cos [u *sj 
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AH = arc sin [Um ctn y] 
y = arc cos [?. 6] 

where tjj = measured horizon crossing time 

t = measured Sun crossing time 
s ® 

= dihedral ai^le between spin axis/Sunline plane and spin axis/ 
horizon vector plane (Frame GI resolution is used with ail vectors) 

As = correction due to tilt of Sun sensor reference slit plane 

AH = correction due to tilt of PAS horizon detector plane 

NOTE: t and t are measured by different sensors separated by dihedral 
a^le . 

Model 6— Earth Widtli Model 


y = t - t 
^ HO HI 




A = arc tan 
do 


✓V /V /\ 

S . (Dj X Dq) 


®o"^‘ ^I> 


where t.„ = measured horizon-in crossing time 
" HI 

SlO ~ horizon-out crossing time 

A^g = dihedral angle between spin axis/horizon-in plane and spin axis/ 
horizon-out plane (Frame GI resolution is used with all vectors) 

Dj = horizon-in crossing unit vector 

Cq = horizon-out crossing unit vector 
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Model 7 — Small Target Model 


y =y 


8 


y = -Av + arc cos 



where = measured PAS horizon detector angle 
Model 8 — Sun to Earth Mid-Scan Dihedral Angle Model 


y = O.o [tjjj + - tg 


■it 




”] 




arc tan 


j § ♦ [0 X Dj] I 


S . Cu X Dq] 


( ‘ 


As = arc sin (tan C ctn /3) 

fi - arc cos (U • S) 

All = arc sin (tan c ctn y) 

H 

y = arc cos (cos c cos (y + Ay)) 
H 8 


where frame Cl resolution is used with all vectors. 
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Models 2, 3, and 4 require tlie transformation mati*L\ A whose components are: 


A = 


1 


(2 2 


S S CM - S 80 S S C0 -f* s 80 
1 3 2 2 3 1 


-S^ Sg S0 - Sg 00 -Sg Sg 80 + Sj C0 


'-K-a 

K-a 


C0 

80 




2 . 2 
2 


/■2' ■ “2 
®2 A ®2 


/ 2 ^ 2 
® 3'^®1 ^ ®2 


where s. 




0 = 0^ + W ft - 

and c and s signify cosine and sine, respectively. 

The equation for the transformation matrix B shown in Models 3 and 4 is 

ii 




cAy C 0 JJ - sAy spjj S 0 JJ -cr^j 80^^ sAy 00 ^^ + cAy 80^^ 

cAy S0JJ + sAy se^i c<t>^ ce^^ c0^ sAy 80^^ - cAy sCjj ©0^^ 


-sAy Cf 


H 


B€ 


H 


cAy Of 


H 


/'GI 

D (used for Models 4, 5, 6, and 8) is computed by solving the following set of 
equations: 


-gGI ^ ^GI _ ^ ^ 


H 




. 6°' = . cos (p +A0) 
c 




The clcpciidoncc of Uic observation models on the stale vector components is 
shown in Table 1-1. 
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The equations of tlie basic recursive processing algorithm implemented in 
OADIAS arc: 


K. = fw/ + p. , G. J P, , G 

i Li j-1- 

sxl 

T 1 


JR 


S5d I I 


'’J = 

sxs 


s;d 

Sy 

= reference vector used in processing observation J . 

The equations f r the wei^ting factors (w to w ) used in processing the 

(1) (o) 

eight observation models are in Tabic 3-2. 
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APPENDIX A - Ot^SKUV ATION 

PAltTlAT. DEPWVATIVES 


This apjx'ndix presents equations for the derivatives dy ./3x. which are used 

0 1 

in the OAlilAS observation models. Sections A. 1 thi*ough A. 8 present these 
eo' .ons for each of the eight models in numerical sequence. Section A. 9 
presents eq. ations for the derivatives dA/dx. , dB/dx^ , and dD/dx^ of 
matrices A and B and unit vector O which arc used in the preceding eight 
sec ons. 

In the interest of conciseness, the derivatives by /bx.. which are zero are 

c i 

not shown explicitly; all others are. 

Tu'O equations Tor each derivative are included with Models 1, 2, 3, 4, and 7. 
The first equation in each case is in a form which is intended to provide max- 
imum understanding of the geometrical far ors (mainly the relationship between 
significant unit vectors) which influence the derivative's value. The second 
equation in each case is in a form which is analogous to that used in the OABIAS 
coding. Single equations which serve both functions are given for Models 5, 6, 
and 8. The second equations are not idr>r.tical to the coding in all cases. In 
particular, the appendix employs a fuller use of rotation matrix notation with 
Models 3 and 4 than does the coding. 

The main notation and notntional techniques used in this appendix were defined 
previously in Section 3. Attention is called, in particular, to Figure 3-1 which 
shows the relationships bebveen the various coox*dinatc frames and to Table 3-1 
which defines these cooi'dinalc frames and the significant unit vectors. Lowur 
case letters (Uj , s^ , etc.) arc used to signify the components of unit vectors 
(U, S. etc.) along the axes of Frame GI. 'I’he definitions of the 12 slate vector 
elements x^ ai.d their alternate symbols are given in Section 3.10. 
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A. 1 MODKL 1— SUN ANGLE ’ODEL 
Subroutines : SAN FUN 
Observation Equations: 




+ arc cos X 


where X = U • S sec c 
Partial Derivatives: 


dye 

II 

1 

II 

^ /N 

S • (eg X U) 



• ej /cos^ e - (U • S)^ 
o 


11 

/O 1 
II 

-§.(^^x0) 

^^2 


(S' /cos^ c - (0-§)2 

U 



dXg ~ dAjB 




/N /\ 

U • S tan € 


/ 


cos'* C - (U • S)" 


X tan c 
\/ 1 - X ^ 
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A.2 MODEL 2— SUN SIGHTING TIME MODEL 
Subroutines: STM FUN, AMATKX, APARTS 
Observation Equations: 


z = 



^ /*v 

y = U • N = > u. (a„. cos € + a^. sin C) 
c S ^ 1 ' 2i 3i 


Partial Derivatives; 


dy dy dN_ 3 / da„. 

c •'c ^ S ^ I ^ 2i . ^ 3i 

T— = = U • -r-— = > u. I cos € — — + sin € 


da 

ds 


dy dy dN„ 3 

-^c ^ S I 2i . 

r-- = -r— = U • = > u. (cos e — — + sin c 

*'‘2 *=2 *®2 ‘\ ^®2 


da 

ds 


") 


dy 3 / da^. da A 

T~ = -r— = -U • cos e = ^ u. [cos e -rj-^ + sin e -r^) 

dy dy 3 / da da \ 

V 1=^1 


dy dy 
dx 


y ay ^ ^ \ 

^ ^ , G • (cos c S - si„ "i fr-os < - Sl„ c 
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whore c,' - T cos 0 + ^ sin 0 
e ' = -T sin 0 +7 cos 0 

/V >\ xs. ’ 

NOTE: dNg/ds^ (also dL/ds^ and dNjj/ds^ which appear in Sections A. 3 
and A. 4) arc the partial derivatives which would be detected bj' an 
observer whose orientation was invariant I'clativc to Frame GF. These 
partials do not appear explicitly in the OABIAS coding. Algebraic ex- 
pressions for them were not developed in the present stud^'. 

A. 3 MODEL 3— NADIR VECTOR PROJECTION MODEL 

Subroutines: LRFUN, EPHEMV, VECROT, AMATRX, APARTS 

Observation Equations: 


z = -y^ 


✓V /N 

y = cos p + R • L 
c 


Partial Derivatives: 


dye 

dye 


»L 

dXj 

-ds^ 

= R* 



dye 

/N 


^*2 


= R • 



R 


ds. 




S /$yfx 1 ^GI.TdA'^fSC 
- — r: — — =R*(SxL) = -:nrr H "n — l 

dx dv |u| 


^ , fi . S;. . ^ ^ J.H- 

dXj dAy 1 |ji| dAy 
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"'c ^ ^ 1 ^GI.T T >^Iir 

>r — = — r = R • (S X L ) ^ — sr- R A — r- L (in-crossing) 

* IHI 


5 


H 


d0 


H 


= 0 (out-crossing) 


dye dy^ 

r — = — — = 0 (in-crossing) 
dx O 

6 dSy 


y\ ^ 

= R • (S X 


SB 

1 ^I,T *T O ^HO* , ^ 

Tr=r R A — ^ (out-crossm^ 
R O 


d0 


H 


dy dy 
^ c c 

■r — = -c: — = -sinp = -cos p sin Ap - sinp cos Ap 
dx^ dAp c c 


1 bA^ ^SC 

dXg doj j^l 


db? 


®^c ^^c /N 1 -^LT T dB ^ 

T-^ = = -R • N cos (Ay + y J = ^ R ^ A 3 ;;— L’ 

dx,, dc„ H ® R 


11 H 


H 


A V yN • ^ • 

^c ^c dR R dcosp cv dR R 

L • + — r r^; — = L • T — + — r— tan p sin p 

dXi2 dAt dAt I II I ^ju[ dAt (r j c 


where 


- 13L»’ 
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R = 



dR 

dAt 


-^[V- (V • R) R] 
R 


p=p^+ Ap 


siL.* =(R^+h)/lRl 

w 

cosp^ = (l- (R^+h)VlR|^) 


where R = Earth radius at horizon crossing 
h = effcctiw height of atmosphere 

NOTE: I signifies an in-crossing and O signifies an out-crossing; neither 

symbol is used with equations or variables which are applicable to both 
crossings. 
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A. 4 MODEL 4— HORIZON CROSSING TIME MODEL 
Subroutines: LNFUN, AMATRX, LPARTS, APARTS 
Observation Equations: 


z = 



y^ = D.NH 


Partial Derivatives: 


dS ^ ^ ijGI ^>GI,T dA^ ^SC 


H 


SD ^ H ft /dD\ ftGI ftGl.Td^ftSC 

^ ds„ H ds^ Ids^/ H ds^ H 


dN 


axj, dSj C,Sjj H Ss^ 


if 


as^ H 


/N ^ /nGI.T dA*^ ^SC 

T-^ = T7^ = S«(N„xD) = D * ^7-N„ 

dx_ ' H ' H 

3 o 


by by R • (N„ x D) sin (y_ + Ay) cos c 
C _ 11 b 

dx. ” dAy 

4 


H / dfi ftGI 

■“ = W/ 


ft • (I> 


^•'^c A /s /^ /sGI.T T -^III 

A ^N,„ (in-crossing) 


II 


b<t> 


II 


= 0 (out-crossing) 
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/\ /s yv xvGI T T '^IIO 

= s . (Njjq X D) = A* nJJq (out-crossing) 

d0„ 


^ ^ S .(N^XD) Stop .^vGI 

dx^ dAp- g.(6xS) * H 


dx„ dw~' bu) H 

9 . O 


^ X D) cos (y„ + Ay) sin e 

C C H O H ^ 

^ " R.(Dx§) a 


/dD\®' OGI ^^GI,T .T bB 

= \^) •''h 


een H 


dye dy^ Njj • (S X D) 
dXj2 R • (O X §) 




OGI 

H 



= defined in Figure 3-9 
V = spacecraft velocity vector 
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NOTE: In the models in which D appears, the dc rivatives arc cvaluMed at 

the measured horizon crossing times. As a result, it would be possible 
to replace D by L in the dy^/dx^ equations. dD/dxj however can- 
not be replaced by dL/Sxj , 

A. 5 MODEL 5— SUN TO EAItTH-IN AND SUN TO EARTII-OUT DIHEDRAL 
ANGLE MODEL 

Subroutines; DIAFUN, LCOMP, LPARTS 

Obser\'ation Equations: 


A, =arctanfl2E\ 

"d5 \B0T/ 

TOP = §.(U x6) 

BOT = U» D- (U • S) (D • S)* 
As = arc sin [tan e ctnjS] 
fi = arc cos { U • S^) = f (s , s ) 

X M 

All = arc sin [tan f ctny] 
y = arc cos (D • S)* 

*D • s - cos fjj cos (y^ + Ay) = f (Cjj. Ay) 
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) 


dy dv 
c ■ c 



(C I ds^ 


sec As CSC 


3 


J3 tan e U 





1 


sec As CSC 


3 


iStan € U 



!>£ 

^^4 


dAy 




(0 

dAy 


sec AH CSC y sin € 


H 


sin (yg + 


Ay) 


I 


% 1 

f = - ~ (in-crossing) 

5 a*„ 

= 0 (out-crossing) 


dy 

_ = — ^ = 0 (in-crossing) 

6 

= - ~ (out-crossing 




by bv 
c _ c 

bx^ ~ bAp 



O) dAp 



dXg do; w 
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}_ 

(jC 


sec As ctn jS see 


2 


c 


dA 


^ 1 i ^"d5 


+ sec AH ctn y sec c 


H 


- tan c„ CSC y sin c„ cos (y_ + Ay) 

li H D 


-]| 


where 


dx^ dAt to dAt 


dA 


d5 


dx. 


_1 / 1 dTOP _ TOP dBOT\ 


/TOP \^ 
^ * \BOT/ 


BOT^ / 




SBOT 5 D ^ A. /y\ SS\ *^ / 

557 • ^>(“ • Is;) * • s>f " <>- 


+ cos (y + Ay) sin € 

9 


II dx, I 


+ AX) cos c„ 4^ 
s H dXj 
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(i= 1) 



. 1 _ _ 

-a i 

(1 = 2) 

- S.ii' 2 


= 0“ 


(i>2) 

dA7 

■ 

(i = 4) 


= 0 



dXi 

(i=ll) 


= 0 




NOTE: 3 D/ax^ ^ 0 in any of the above expressions. 
A. 6 MODEL 6— EARTH WIDTH MODEL 
Subroutines: DHFUN, LCOMP, LPARTS 

« 

“ w (^dO ■^^H “^II) 
A-12 
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= arc lan 



TOP = § • ( 6j X D^) 


DOT = D, 


Do-(S 


✓N /S 

Dj)(S 


“o> ' “l 


ft 2 2 

D - - cos €„COS 
U H 


Partial Derivatives: 







" 


dx “ dAy ~ 0) dAy 
4 



by 


bt> 


c 

r 

H 


iti 


(in-crossing) 


= 0 (out-crossing) 




0 (in-crossing) 


= — (out-crossing) 
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_ 1 

“ dAp U dAp 


&Xg “ dto u> 


*’‘n ^'h “ 


dXj^ dAt U dAt 


where 


dA 


d6 


Sx, 


_1 l_i__ STOP _ TOP dBOT\ 

BOT^ ' 


2 \BOT 

^ \BOT/ 


dD. 


SD, 






dBOT /N ^ 2 

^ + 2 cos c„ Sin c„ cos^ (V + Ay) 

dx O ox. I dx H H '^S dx 


2 dAv 

+ 2 c'^s €jj cos {y^ + Ay) sin (y^ + Ay) 
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(1 = 4) 
(i^4) 

( 1 = 11 ) 

(i?^ll) 
d S' 

For - — , see Model 5. 

A. 7 MODEL 7— SMALL TARGET MODEL 
Subroutines: £?BFUN, EPHEMV 
Observation Equations: 


dAy 

dx. 


= 1 


= 0 


d€ 


H 


dx. 


= 1 


= 0 




y = “Ay + arc cos a 
c 


where a = -R • S sec c„ 

H 

Partial Derivatives: 
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by^ 

*•'‘2 ^“2 (S • Sj) ycos2 - (U • S)2 



dx^ ~ b Ay 


ax,, " ac 


R • S tail c 


H 


-a tan 


H 


11 ycos2c^-(g- §)2 /l -o^ 




SR 


sec c 


**12 y^s?73-(ft-§)2 yrr? 


H gOI, 


.GI 



For 


ys 

bR 

bAt 


see Section A. 3. 


NOTE: The mathematical similarity between Models 1 and 7 can be seen by 
comparing the above equations with those in Section A. 1. 


A. 8 MODEL 8— SUN TO EARTH MID-SCAN DIHEDRAL ANGLE MODEL 
Subroutines: DIAFUN, LCOIVIP, LPARTS 


Observation Equations: 


2=y-yc 


y “ 0* + 1 jjq) t^ 
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TOPj = S • ( U X Dj) 

/V /\ y\ y\ 

BOTj= U* Dj-(U* S)(Dj* S) 

TOP^ = $*(U xDq) 

/S ^ /V /V /\ /V. 

BOTq-U*Dq-(U*S)(D^. S) 

As = arc sin (tan € ctn /3) 

P = arc cos ( U • §^) = f(Sj^, s^) 

AH = arc sin (tan c„ ctn y) 

H 

V 

/N - 

y = ar^ cos ( D • S) = arc cos (cos c„ cos (y_ + Ay)) = f(e„. Ay) 

Ho H 


Partial Derivatives: 


ds^ OJ \ 


3 „ ^ 

sec As esc p tan € 


»'ii) 




ir \ 

an c U • T-r- I 
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bx. 

4 


1 3 

Sec All CSC y sin Cjj sin (Vg + Ay) 


^yp ^yp -0.5 

~ ~ Ui 


5 


H 


^ K :0^ 

H 


dXg " hit} ~ w 


^yp ^yp 1 2 

r -- — = — sep As ctn B sec € 

ae w ^ 


!Ll-^ 1 


1/2 2 \ 
= — sec AH fctny sec c„ - tan c„ esc y sin c„ cos (y_ + Ay)) 

CO \ H H Ho/ 


dXj 2 ~ 03 dAt 


dA 


ax 


^ = 0.5r^(A_jj(i„)tAj.(out))l 


l_»x, 


where a A ,_/oX. is defined in Section A. 5. 
do 1 


A-IS 



ORiGUm PACE i9 
OF POOR QUALmr 


A. 9 SUPPLL'MEXTAUY PAIITUL DL’UIVATIVES 
A. 9.1 Matrix A Partial Doriv'afive 

The equations for Matrix A were discussed and presented previously in Sec- 
tion 3. 10. For the reader’s convenience, they are repeated below. 


A = 


^11 

*^12 

“l3 

a_ 

a 

a 

21 

22 

23 


a„_ 

i£l 

31 

32 

33. 


1 

cos 6 


s^Sg eos 0 - Sg sin 0 
-s s sin 0 - s eos 0 

X <5 A 

COS 6 


®2®3 0 + Sj sin 0 

-s s sin0 + s COS0 

Xi o X 

s cos 6 


-cos 6 cos 0 
2 * . 

cos 0 sin0 


8^ COS 6 

3 


where cos 5 




+ s^ 
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The derivatives 6 A/dx. can be obtained by direct differentiation. The results 

i 

are as follows; 


hA 

dx. 


ds. 


cos 6 


cos 0 (s^ - 
-s^ a^j/cos 6 


-sin 0 (Sg - s^/Sg) 

-s. a„ /cos 6 
1 21 


cos 6 


sin 0 


-cos 0 (s^ 


cos 0 


+sin 0 (s^ 

-Sj ajj/cos 6 


cos 0 


sin 0 


-(s /s ) cos 6 

X u 


aA ^ ^ 
dx„ as„ 


cos 6 


-sin 0 

-cos 0 (s^ ® 2 ^® 3 ^ 
s^/cos 6 

-cos 0 

+sin 0 (s^ ®2'^®3^ 
"^21 


cos 0 (Sg - Sg/Sg) 
-“i2 ' 


-Sin 0 (Sg - Sg/Sg) 

-“22 ®2'''“ * 


COS 6 


-Sg cos 0 


Sg sin 0 


-(Sg/Sg) cos 6 


♦The use of the partial derivative symbol a here is a misnomer, because de- 
rivatives ds./ds , d cos 6/ds , and d(cos 6)’Vds (a-~l. 2) arc included. 
3 a a a . 
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dA dA 




a 


22 





dA 


dA 


dto o d|p 


A. 9. 2 Matrix B Partial Derivatives 
The equations for Matrix B are: 


B = 





where 



COS0JJ 

sin0jj 


-sin0jj 

COS0JJ 


0 0 


0 

0 

1 




1 0 
0 cos 

0 sin Cjj 

COS Ay 0 
0 1 
-sin Ay 0 


0 

-sin 

cos Cjj 

sin Ay 
0 

cos 
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The partial derivatives of B can be obtained by direct differentiation; these are 


T T 

dx, “ dAy 
4 


-sin Ay 0 cos Ay 
0 0 0 
-cos Ay 0 -sin Ay 


dB dB dB 


axj.- ixj a0„ 


-sin (t>^ 

-cos 

-sin0 

M 

0 


...T 

COS<ijj 

0 


T 

Ay 

0 

0 

0 




as as T 


0 -sine„ -cose 

H H 

0 cos € -sin €„ 

11 u 


Ay 


A. 9. 3 D Partial Derivatives 

The computation of D in subroutine LCOMP was discussed in Section 3. 8. The 
three fundamental equations used in the computation are repeated below. 


§•6 = cos (y + Ay) cos e 
D H 


(A-la) 


• D = -cos (p^ + Ap) 


(A-lb) 


6 * 6 = 1 


(A-lc) 


The partial derivatives of D with respect to the state vector elements x. are 
computed in subroutine LPARTS. The derivation of the LPARTS algorithms 
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starts by differcnlialing liquation A-1 to obtain the following set of perturbation 
equations 


t . dD + D • (IS -■- -sin (y^ + A7) cos Cjj dAT 


- cos (Tg + Ay) sin dc^ 


R • dD + D • dAt = sin (p + Ap) dAp 

c 


, d cos (P^ + Ap) 


dR 


dAt 


dAt 


6 . d6= 0 


where dS is constrained to be 


(A-2a) 


(A-2b) 


(A-2c) 




(A-2d) 


Equation (A-2) defines the perturbation dD in 6 which is produced by per- 
turbations dSj^ , dSg , dAy , dAp » dCj^ . and dAt . 

The present discussion will use the simple symbol R to signify the magnitude 
of R , As before, R will signify the rate of change of R , and 'V will signify 
the sixicccrafl's velocity vector. diR/dAt of Equation (A-2) is the velocity vec- 
tor of U ; later in the derivation, it will be specified by the following equation 



where 
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R 


(A-3b) 


The second term on the right side of Equation (A-2b) is the change in cos p 

due to the change dR in R which is produced, in non-circular orbits, by dAt . 

The following cciuation for this term was developed with the aid of the equations 

for sin p and cos p listed in Section A. 3. 
c c 

d cos (p + Ap) • 

« ' R (A--*' 

/\ 

The problem now is to solve Equation (A-2) for d D . As a tool for accom- 
plishing this, the following 3x3 matrix M will be introduced: 




(A-5) 


The left sides of Equations (A-2a) and (A-2b) do not necessarily imply that the 
vectors arc resolved upon a coordinate frame. With the introduction of M into 
the development, however, it will become convenient (although not absolutely 
essential) to assume resolution upon a specified frame; namely, frame GI. For 
nolational simplicity, however, the superscripts Gl used elsewhere in the 
report (o designate GI frame resolution will be omitted from the remainder of 
this discussion. 

r 
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Arranging Equation (A-2) in a mati'ix format, introducing Equations (A-4) and 
(A-5), and solving for d 8 no»v yields 


-D • d S - sin (y_ + Ay) cos « „ dAy - cos (7 + Ay) sin c.. df 
S H S il II 


dD = { sin to + Ap) dAp - 6 • ~ sin (p + Ap) tanp dAt ^(A-6) 

c oAt iv c c 


Equation (A-6) is not meaningful when M is singular. The condition of 
singular M occurs if and onlj’ if the three vectors D , S , and ^ are in a 
common plane. 

A general equation relating d 8 to the state variable perturbations dx is as 
follows: 


12 


a6= V 


✓N 


(A-7) 


The desired equations for the partials d8/dx. now can be obtained by com- 
paring Ecjuation (A-6) with Equation (A-7). Equation (A-3) is employed where 
necessr ry. I or this work, it is convenient to define the columns of M ^ ex- 
plicitly: i.c. , 




(A-8) 


The resulting dD/dx. equations arc as fol lews: 


d8 

d8 

[-Q a -c. 

L'* ' 

dx^ 

-ds^ 
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-5i 



(A-9b) 



/N 

dP 

dAy 


= -Q^ sin (yg + Ay) cos 


(A-9c) 


/\ /\ 
SD ^ bP 

dx^ dAp 


= 92 sta (p^ 


+ Ap) 


(A-9d) 



d6 


= -g^ cos (yg + Ay ) sin Cjj 


s6 ^ ^ 
&x^2 dAt 



dAt 


.|pin(p^ 


+ Ap) tanp 


% 


(A-9e) 


(A-9f) 


Equations (A-9a) through (A-9f) are similar to those implemented in the OABIAS 
coding. An alternate form which better shows the geometric factors that influ- 
ence the dD/dx. values, however, is possible. To derive these equations, 

-1 ^ 

M is written in the following explicit form: 


m‘^ = (det M)‘^ [rxB §x^] 


(A-lOa) 


where 


det M 6* (^x fi) = §• (It X 6) =1^ • (B X §) (A-lOb) 


Equation (A-10) can be verified by preimiltiplying it by Equation (A-5) to 
produce the identity matrix. ,• 
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The alternate form of the d equations is obtained by substituting Equa- 

tion (A-10) into (A-G) and proceeding as before. The dS/dAt equation can be 
simplified using Equation (A-2d) and (A-lb). The results are 





fixfi/e. ii-1 



dB _ 

aB 

1 





B • (Bx §) 





1 


dB _ 

aB 

\ 

f 




B*(Bx§) 


. s6 . 

Bx 

B sin (y + AT) cos e 
S ri 


SAT ■ 


B • ( B X B) 


dB 

_ 

_ D X S sin p 



“ dAjC 

> B*(Bx§) 



Bx 

B cos (yg + Ay) 

sin^H 




B » (Bx B) 


(A-lla) 


(A-llb) 


(A-llc) 


(A-lld) 


(A-lle) 


_ S X D ‘/^ ^ g cos Ap \ ^ 

^^12 It • ( D X §) [\ Pj. / J 


(A-llf) 


where the dS/ds^ equations are given in Section A. 5. Unlike Equation (A-9), 
the above cqualions can be regarded as vector equations per so, and not merely 
as equations for vector components along the GI fiaime axes. 
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APPKX DLX D - pKim'ATIpN OF U KCt'KSP/ K PUOCESSIKG 
AI^GOIUTIIM USFD IX OABIAS 

B.l INTRODUCTION 

The purpose of this appendix is to derive the equations of the basic recursive 
processing algoidthni used by OABIAS. The resulting equations are applicable 
to any sj'stem which uses subroutine RECURS (Reference 21) or RECUR? now 
employed in O.ABIAS. A generalized least squares approach is taken in the 
derivation. 

Section B.2 discusses the generalized least-squares loss function and then 
derives the batch processing differential corrector equations. This is a 
necessarj' preliminarj' step in the OABIAS equation derivation. Section B. 3 
presents the derivation of the nonrecursive least-squares algorithm. Sec- 
tion B.4 derives a preliminary set of recursive processing equations. These 
equations are jentially identical to the ones in Section B. 3, except they are 
directly applicable to recursive processing. From the resulting equations in 
Section B.4, Section B. 5 derives ttie basic algorithm used in OABIAS. 

References 11, 20, and 22 through 24 provide the bacl^ound for the derivation 
methods used in this appendix. * 

The main notational techniques used in the appendix are as follows. All sym- 
bols with an underbar, but no superscript T, are column vectors. Super- 
script T signifies the transpose of a column vector or of a matrix. Matrices 
are designated by uppercase symbols with no underbar. Lowercase sjanbols 
with no underbar are scalars. The derivative of a scalar with respect to a 
vector is considered to be a row vector. The dimensions of matrices and 
column vectors are indicated the first time they apiiear in equations. Tlie 
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notation dislin<;uishcs between variables and algebraic expressions which are 
used to compute tlic values of these variables. For example, in i4ie equation 

Y = Y (M, X, t) 
c e 

Xq • ii t U is a laiown algobraie exijressioa. The numerical values of the 
elements of Y^ are eomputod by inserting values of ^ , 2L » ^ i*^t° i^* 

Parentheses are used only in this manner. 

B.2 DISCUSSION OF THE LEAST-SQUARES LOSS FUNCTION 

The purpose of the algorithms derived in this appendix is to compute optimal 
estimates of the state of a sx'stcm using, as primary input, a sequence of 
observations performed on this system. The system state is assumed to be 
constant and is modeled as an sxl vector. The symbol X will denote the 
unknown true state vector. A distinction will be made in the early stages 
of this discussion between "arbitrary" estimates of X , which will be denoted 
by X , and the optimal estimate X* . The symbol X^ signifies an a prio,ri 
estimate of X . 

A generalized least-squares criteria will be used to develop the algorithms 

/\ /\ /V 

for computing X** 21* » by definition, is the X which minimizes the 
following loss fimction: 


W. 


/j = 0.5Z, ^ 

pxp pxl 


-1 


0.5 jx, - X 


» s ^ X - J 

-o| 0 j-i -oj 

sxs sxl 


(B-la) 


where 




!LV 


nxl 


t) 


(B-lb) 
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and 

Yj = Y^ (M^, t) (B-lc) 

--cl " -cl ^1^ 

Tho 1 subscripts in the above equations signify that the estimate is to be 
obtained from a set of data designated as set 1. For notational simplicity, 
these 1 subscripts henceforth will be omitted until the development reaches 
the point where they can i ve a useful purpose. 

Z Is the residual vector. Y is called the real observation vector or just tlie 

observation vector. Y is the model i >servation vector. It if- important to 

o 

distinguish between Y and the measurement vector M . The elements m^ 
of M are the basic scalar measurements generated by the onboard sensors 
and preliminary ground processing operations, and are the primary inputs to 
the composite altitude determination system. The elements of Y are 
the inputs to the estimation subsystem of the attituue determination system. 

The elements y^. of the model observation vector Y^ are the predicted 
values of y^ using a mathematical model of the satellite and its sensor sys- 
tem. Under ideal conditions, each y^. would be identical to its corresponding 

y, . The residuals z. between eac^ v. and y , provide an indication of the 
t 1 * i ci 

error in tlic esiimate of X . 

In most estimation studies, a distinction is not made between Y and M . 

Also, Y is not normally considered to be a functioii of M as shown in 

OT ^ -I 

Equation (D-ld). The OABIAS estimation system is unorthodox in several 
respects, howcv?r, and the purpose of the apparently general format of Equa- 
tions (li-lb) through (B-ld) is to make the forthcoming derivation completely 
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applicable to OADIAS. The most unusual feature of OABIAS is that in several 
of the sensor models, y. is known exactly and the M-dependenej' enters soicly 
through y . 

The loss ‘unciion I of Equation (B-1) is a blend of the new measurement data 
M and the a priori state vector estimate X . W and S are weighting 
matrices whose function is to establish the relative weights to be assigned to 
the individual residuals and to ti.w ^cate vector deviations (x^ - . 

Hie term "least squares" often is used to signify loss functions in which W 
is diagonal (sometimes with all diagonal elements identical) and/cr is 
zero. The present development will not make either of these restrictions in 
its early stages and, thus, can be considered to be a generalizea least-squares 
approach. 

W and S are assumed to be symmetric. Tor X to be meaningful, both W 
o 

and S must be nonnegative definite. The further restriction that W be 
o 

positive defini‘t; usually is iustified. A positive semidefinite W would imply 
measurement data to which no weight a{ all is to be attached. A positive 
semidefinite is meaningful and acceptable, unless an X* algorithm which 
requires its inverse is developed. A positive semidefinite implies that 
the state vector contains components whose values are completely unknown 
a priori. 

Estimation algorithm derivations which employ a least-squares approach some- 
times assume a priori that W is the inver.se of the covariance matrix R 
of the errors in tlic obscxwations or residuals and that S is the inverse of the 
covariance m; irix of the imccrtainty in 2S^ • The presen. ierivation, however, 
does not al'a a anv ncccssaiw stat al significance to S and W . They are 
rcgaixlcd only as weighting matrices :uid arc .subject only to the mathentatical 
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restrictions noted in the above paragraph. W = R ^ and usually arc 

considci’cd, without verification, the optimal weighting conditions which make 

X* , which satisfies the generalized least-squares criteria, truly the optimal 

estimate; However, such statistical restrictions on W and S are considered 

o 

here to fall fundamentally outside of the basic least-squares mathematics. 

In principle, a distinction should be made between the sensor event times t , 

s 

t , t and the time variable t which is included in Equation (B-1). t is 

intended here to signify the time or times at which orbital or ephemeris data 

is evaluated to generate the Sun and central body reference information reqaireti 

by the estimation system. These times are not necessarily identical to tlie 

telemetered sensor event times. For OABIAS application, however, the point 

is not a significant one, because the program does use the telemetered sensor 

event times to generate the Sun and central body reference information. For 

this reason, the model obseiwations y . of all the OABIAS models are func- 

Cl 

tions, to at least some extent, of the elements i , t , or t _ of M . In 

s HI HO 

the interest of notational simplicity, the t- and M-dependence will not be 
shown explicitly in the remaining equations of this appendix. 

B.3 DERIVATION OF THE NONRECURSIVE LEAST-SQUARES ALGORITHM 

A necessary condition for X * to be the estimate which minimizes the loss 
function jt of Equation (B-1) is 


/s T 

=0 


ax 


sxl 


(B-2) 
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/V ^ /V 

Differentiating Kquation (B-1) with I'espect to X » setting X to X * » 
using Kquation (B-2) yields 


^(X^) W Z (X*) + S [x* - X 1 = 0 

“ ' " " SXl 

sxs 
sxp 


(B-3) 


The problem now is to derive a method for solving Equation (B-3) for 

/V 

The usual Newton-Raphson procedure expands Z(X*) in a Taylor series 

about X . A slightlv more general approach will be employed here, however , 
~o 

in order to obtain equations which can be used in the recursive estimator deri- 
vation. Instead, each element z. of Z will be expanded in a Taylor series 
about an arbitrary reference to be denoted as X . The p Xo-*® 

assumed to be identical to X nor to each other. They are assumed only to 

~o 

be sufficiently close to X* to partially justify the forthcoming series tnmeation. 
The Taylor series result can be placed in the following form: 


Zi2*) = Y - - Gj^. X* + H.O.T.(X* - X^) 

pxl pxl pxs 


(B-4a) 


where 




cp Up I 



(B-4b) 


(B-4c) 
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(B-4d) 


and 


Sm - al' S r/ = 

8 X 1 

The last e>q 3 rcssion on the right side of Equation (B-4a) signifies higher order 
terms in the series. In addition to Equr*’’''n (B-4a), the Taylor series expansion 

of the partial derivative of 7^ with respc X * needed. This can 

% 

' be obtained by differentiating Equation (B-4a). 


^1 (S*) = + H. O. T. (B-5) 

pxs 


Substituting Equation (B-4) and (B-5) into (B-3), dropping the higT ‘ order 
terms (II.O.T.), and performing some minor algebraic manipulat yields 


h * "'i °ir] 1 % - 1 = "'i 1 ^ ciR 


*®1R"®1R^o! 


(B-6) 


The previously noted 1 subscripts which signify that the estimate is obtained 
from data set number 1 have been inserted into Equation (B-6) for convenience 
in future references to the c(iuation. Also, the oplimal estimate will now be 
signified merely by 2 rather than by X*. 1 u.a simplification is being made 
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because the * no lonp;or serves a useful purpose. Also, it brings the current 
notation into closer agreement with that commonly used in Computer Sciences 
Corporation reports. 

For the remainder of this section, the work can be restricted to the case where 

all reference vectors X_. arc identical; X„. = X where i = 1 to p , Equa- 

-Ri “Ri ~1R 

tion (B^) shows that can be replaced by X^^^ . 

Let be the sxs coefficient mati^x on the left side of Equation (B-6). 

Equation (B-6) possesses a ui\ique solution if, and only if, Q is nonsingular. 

T ^ 

Assuming that both S and G,„ W G,_ are at least positive semidefinite, a 

sufficient condition for nonsingular is that one or both of the pair be posi- 

i T 

tive definite. Positive definiteness of G^j^ signifies that X is observ- 

able from the set 1 obse n ations alone. Nonsingularity of can be assured 
by choosing to be positive definite. 

y\ 

In the nonsingular Qj^ case. Equation (B-6) can be solved for , and the 
result written as follows: 


where 


A. 

X =X + 
”1 -o 


■'i 

sxp 


W-Y +G 

I 1 clR 1R|-1R -ojj 


(B-7a) 


K = P G^ W 
1 1 IR 1 

(B-7b) 

sxp 


P = 
1 

(B-7c) 

sxs 


Q = S + G^ W G 
0 IR 1 IR 

(B-7d) 

B-3 






original 

OF POOR QUALfTY. 


At this point in the development, no special significance can be attached to 
and ; they are merely convenient matrices for use in the equations. 

Equation (B-7) constitutes a slightly generaliired form of the usual nonrecursive 
least-squares algorithm. \Vlicn employing this equation for batch processing, 
the reference vector normally is set equal to the a priori estimate . 

This enables the matrix on the right side of Equation (B-7a) to be dis- 

carded. It will be recalled, however, that when deriving Equation (B-7), it 
was necessary to ti*uncate the higher order terms of Equations (B-4a) and (B-5). 
Thus, the result of Equation (B-7a) will minimize the loss function I 
defined by Equation (B-1) only when is sufficiently close to Xj^. (now 
= X^) that tlie er-ror due to discarding the higher order terms is negli- 
gibly small. An iterative processing operation, called differential correction, 
commonly is performed to overcome this difficulty. With this technique, data 
set 1 is passed through the processor several times. The state vector esti- 
mate obtained in an}' given pass becomes the a priori estimate for the following 
pass. In order to delineate the operation mathematically, let superscript X 
denote the Xth pass through the processor. Equation (B-7) then can be con- 
verted to the following form: 


/^X 1 ■^x 

^ 


(B-8a) 


nl 

(B-8b) 


''’i 

(B-8c) 



(B-8d) 

<5^ * s 

rw, 0, (s«) 

(B-8c) 
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5 ^ 



(B-80 


The process is said to have converged when the elements of become 

negligibly small. The iterative operations tend to degrade the meaning and 

usefulness of S . Equations (B-8a) Uirough (B-SQ are equivalent to those used 
o 

by GCONES, except GCONES does not include . 

B.4 DERIVATION OF A PRELIMINARY RECURSIVE LEAST-SQU.ARES 
ALGORITHM 

The preliminary recursive least-squares algorithm derived in this section is 
not the algorithm used in OABLAS. However, it is a necessary b 5 '-product in 
the present derivation of the OABLAS equation and is of some significance in its 
own right. 

Assume that data set 1 has been processed, iteratively or noniteratively, and 
that the final results and P^^ of have been saved. Let data set 2 now 
be received. Set 2 contains q scalar observations (q ^ 1) . Let the data set 
which Is comprised of both sets 1 and 2 be den''*^ed as set . Set contains 
r = p + q scalar observations where r s 2 . The problem now is to obtain a 
new estimate utilizing all the observations in set H ; the subscript 2 
signifies that the estimate employs all data sets up to and including set 2. The 
batch processing approach to the problem would process the composite set 53 
in the same way that set 1 was processed. With appropriate changes in sub- 
scripts, the equations developed in Section B. 3 arc applicable to this method. 

In essence, tlic methods minimizes a loss function I defined bv 


I =0.5 7^ W Z +0.S|X -$!^S jXo"^! (B-9a) 

_ I -2 -oj o j_2 -oj A I 


T '5: 

rxr rxl 
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where 


„T 1 „T I 

h ~ 1-1 -2 1 






(B-9b) 


(B-9c) 


The major drawback of this approach is that the resulting algorithm [Equa- 
tion (B-7) or (B-8)l with appropriate changes in subscripts) includes vectors 
and matrices with dimensions equal to the total number of scalar observations r 
in set ]C , This increases the computation and storage requireme>'*s. If the 
approach is continued when additional data sets come in, these requirements 
will increase withoi i bound. 


The recursive processing approach to the problem uses only the new observa- 
tions Y and the results $ and P or Q obtained from processing set 1. 

X XX 

With these approaches X« estimate which minimizes a loss function l„ 

defined by 


/2 = 0.5 


7*^ 


W. 


^2 


+ 0.5 


( ^ /s 
1^2-^ 


I 


-II 


qxq qxl 


1 

sxs 


(A. ) 


(B-IO) 


As will be shown Inter in this section, the recursive processing appi'oach 

requires that = 0. Comparison of Equation (B-10) with the loss function 

jfc of Equation (B-1) indicates that the batch processing equations developed 

in Section B.3 arc applicable to this case with appropriate changes in subscripts 

and thus constitute a valid recursive processing algorithm. However, tlie 

problem of selecting remains. Arbitrary selection of , independent 

of S , is undesirable. A better approach is to establish S such that the 
o 1 
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X cor.ipiitcd by the recursive method is identical (at least under certain con- 

A 

ditions) to the estimate obtained using the batch processing technique which 
minimizes . An algorithm for computing Sj^ as a function of is 
needed. The remainder of this section is devoted to this problem. 

The derivation starts with Equation (B-G). Altc ring the 1 subscripts to malie 
the equation applicable to set H yields 




r X s 


‘X -X f = 
(“2 -ol 


= G* W W -Y +R 

-cZR -zr 


r X 1 




"^ER^ol 


where 



(B-12a) 




(B-12b) 


T 

Y 

-cER 


yT j 

-clR-c2R( 


gT ^ I T T I 
®ER (-1R -2Rj 


W, 


w 


12 






(B-12c) 


(B-12d) 


- (B-12e) 
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For use below, let the corresponding set 1 equation (B-G) be x*ewritten in the 
following form: 


t 




3, 1 ^ = S $ + 

IRj ~1' o "o 


T 


K' -Y + 
|-1 -clR 


^ir1 


(B-13) 


Equation (B-11) now is written in its full form by inserting Equation (B-12). 



Performing the matrix multiplications in Equation (B-14), canceling identical 
terms where possible, and rearranging the result slightly yields 


k * “IrVir * <=Jr "iSr * <=m " iAr * “ -2 


r]*2 


■ ^ “Ir*! 1^1 - Vr • Sir') • KR"Ii K -2 cir * 5.r!]<B-‘5) 

* t'L '■'zIlz-XcaR'esRl] 


Development of a recursive algorithm requires elimination of terms involving 

X , Y, , Y , „ , and B.^ from Equation (B-15). The terms within the 
-0 -1 -clR ~1R 

first set of brackets on the right side thus are unacceptable. These, however, 
are identical to the right side of Equation (B-13) and hence can be replaced by 
the left side of (B-13). The remaining unacceptable terms in Equation (B-15) 
are those inside the second set of brackets on the right side. These involve 
the coupling matrix W . To obtain the desired recursive processing algo- 

X M 

rithm, must be restricted to zero. 
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A restriction that all reference vectors in set 2 are identical 

— 2Ri 

(X-,^. - , where i = 1 to q) now is acceptable. This enables B to be 

2R 

replaced by • The concept of nonidentical reference vectors was a 

mathematical tool used to avoid the restriction that the same Taylor series 
expansion reference vector be employed for both set 1 and set 2 . The tool now 
has served its purpose in the derivation and is no longer necessary.- 

The desired recursive processing algorithm is obtained by implementing the 
above-noted operations into Equation (B-15) and performing a few minor 
additional algebraic manipulations. The result is 


where 


^2'^" S k-Xc2R*°2R !^R-Slj| 


sxq 


qxs 


•= 2 = *’2 Sr ^ 

sxs 


^2 “S 


-1 


Q = Q + G W G 
^2 ^1 2R 2 2R 

SXS 


(B-16a) 


(B-16b) 


(B-16c) 

(B-16d) 


Q = S + G. W G 
^1 o in 1 IR 


(B-17) 


Equations (B-7) and (B-16) constitute a valid method for processing data sets 
1 and 2 sequentially. They' can, in fact, be gcncrali^,cd immediately for 
sequential processing of an arbitrary number of data sets. Before doing this, 
however, it Is desirable to answer the previous question concerning the 
correct weighting factor for Equation (B-10). It should be evident that 
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direct use of ICquation (B-10) to derive an optimal estimate with X^^ and 
as the a priori inputs would yield results identical to Equation (B-7) obtained 
from the I equation, except the 1 and 0 subscripts would be replaced by 2 
and 1, respectively. Equation (B-7d) would now be 




(B-18) 


The r'>sults obtained from Equation (B-10) should be identical to those of Equa- 
tions (B-16) and (B-17). Comparing Equation (B-18) with Equations (B-16d) 
and (B-17), it is concluded that 

Equation (B-1 9) is the desired equation for updating S to obtain S, . It is 

o 1 

evident that this result can be extended to yield S == Q , S = Q , etc. 

A 2 3 3 

Because the variables Q can be replaced by S , their use has become super- 
fluous. Also, superscripts 2 and 1 of Equation (B-16) can be replaced by more 
general ones j and j - 1 . Therefore, the final recursive processing equations 
of this section are 


As /N ( 

, + K, Y. -Y . 
“j j-1 j -J "cj 


+ G Jx 

cjR jR |-jR -j-ll( 


K. = P. G^_ W. 
) j jR 3 


(B-20a) 

(B-20b) 


P^ = S"^ (B-20C) 

T ' 

S =S. , + G.„ W. G.„ (B-20d) 

j j-1 }ll J jR 
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B.5 DL'RIVATION OF THE RECURSIVi: IJ^iAST-SQUARES ALGORITHM USED 
IN O A BIAS 

Equation (B-20) oi* tlie equivalent usually is considered to constitute a poor 
recursive processing algorithm because the inversion of an sxs matrix is 
required at each step. They are not used in OABIAS and hence constitute only 
an intermediate result in the present development. 

In deriving the OAElAS algorithm, Equations (B-20c) and (B-20d) first are 
combined into 
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Now assume that W. is nonsingailar and manipulate 
J 


using the matrix identity 


the term within the brackets 



W = 


CW"^ D] 


(B-24) 


The result is 


r -1 tI”'’ 

P=P -P G Iw +G P G G P (i5"i 

t j-1 j-i VL j V j-1 jR] jR i-1 ' 


Equation (B-25) now is substituted into Equation (fi-20b> to give 


25) 


i ' V j* - b"' * °)B ^,-X «Jh] ' -,-x| °l 


G P . G‘ JW. 
JR j-1 jR j 


To finish the derivalion, Equation (B-2G) is substituted into Equation (B-25) to 
Simplify the cc. 
algorithm are 


Simplify the cciuation for P. . The final equations for the new processing 

J 


•^1 



Vlll 

(B-2'^a' 

T P "“1 T 



“ °JR n °jR ^j-l ^juj 


(B~27b) 

u 

C 1 

1 

1 


(n~''7c) 
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Assuniins that \V. is prusscd lo the recursive processor, a processor which 

employs liquation (B-27) must imert two p.xp. matrices in each processing 

J J 

step j, wl.rre p. is the number of scalar observations on set j. Only one 

^ -1 
p.xp matrix inversion is needed il ..\e processor is given . However, 

the observations can be px'ocessed one at a time in applications whevc all 

off-diagonal elements . where fi ^ v , of the weighting matrix' can 

be made zero. This method avoids matrix ii Jrsions and minimizes computer 

storage requirements. Assuming the observations to be processed ot e at a 

time enables Equation (B-27) to be converted into the following form 




(B-28a) 


°iR ' <^iR> 

sxl 


r -1 T 

: = w. +G.„P. G.„ p. , G.„ 

j [ J “jR -j-l -jRj -3-1 -]R 


sxl 

sxs 


I -K. G.^ 
-j-jR 

sxs 


j-l 


AX, = K. y. - y \St. , - X. J ' 

-j "j I J cjR jR l-j-1 

✓S . 


y\ 


(B-28b) 

(B-28C) 

(B-28d) 

(B-28e) 

(B-28f) 


Unlike the other equations of this appendix, Equ.itions (B-28a) through (B-280 
have been arranged into a possible sequence for utiliz.ation in a computer 
roulir''. Tb-'se equations arc the end result of this appendix and ait: the basic 
recursive esiimalion equations used in OAIJIAS except for implement modifica- 
tions dir..' u<sed in c ion 3.4. 
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The symbols most commonly used in this document arc listed below. Refer- 
ences arc provided to figures or equations in which symbols are defined when- 
ever such figures or equations occur. The most common superscripts, 
subscripts, a» ’ overhead symbols arc listed at the end of tlic glossary.' 

S\Tnbol Definition 


A 



A 

liV 

B 


B, 



DOT 

b 

tiv 


A 3 X 3 matrix which transforms vectors from 
frame GI resolution to frame SC resolution 
(reference Figure 3-1) 

A dihedral angle 

The Model 5 dihedral raigle (reference Figure 3-10) 
The Model 6 dihedral angle (reference Figure 3-11) 
The Model 8 dihedral angle (reference Figure T -13) 
The elements of A 

A symbol used to indicate either B^ or B^ in equa- 
tions that are applicable to both 

A 3 X 3 matrix which transforms vectors from 
frame HI' resolution to frame SC resolution (refer- 
ence Figure 3-1) 

A 3 X 3 matrix which '^ransforms vectors from 
frame HO* resolution to frame SC ’•esolution (refci.- 
ence Figure 3-1) 

Defined in Equation (B-4d) 

The denominators in the equations for tan A^^ , 
tn.. . or t;m A^^ 

Th" elements cf B 
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Svmbol 


D 


/K 

D. 


D. 


Det 
Det M 


“l' "z- ^ 


dM 


dm 


dY 




/V 

e 


C 

a 


-Ak /s 

®S1’ ‘‘s2’ ®S3 


■§, 6,6 

1’ ?/ 3 

“1- '2- 


M pit '^11 

II* 12’ 13 

e’* €^' 

Ol* 02’ O', 


Definition 

Cosine 

A symbol used to indicate either Dj or Dq m equa- 
tions tliat arc applicable to ’x)tli (reference Figure 3-9) 

The unit vector directed from the spacecraft toward the 
horizon-in crossing point on the central bodj"^ surface 

The unit vector directed from the sp:icecrafl toward the 
horizon-out crossing point on the central body surface 

The determinant of the 2x2 matrix in Equation (3-72a) 

The determinant of M 

y\ 

The components of D along the axes of frame GI 

T1 . .'ector of the errors in a set of measurements M 

The error in element m of M 

Of - 

The vector of the errors in a set of observations Y 

The error in element y. of Y 

) 

The angle generated by Sun sensor misalignment c 
(reference Figure 3-10); also, expectation operator 

A sj'mbol used to signify a unit vector along an axis of 
a reference coordinate frame 

See Figure 3-9 

The unit vectors along the axes of frame SS 
The unit vectors along the axes of frame GI 
The unit vectors along the axes of fram^ SC 
The ’.nit vectors along the axes of frame HI 
'1 he unit vectors along the axes of frame HO 
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Symbol 


Definition 


FOV 


Field of view 


G 

G. 


GI’ 


GI 

S 

H 

HI 

HI* 

HO 

HO’ 

H. O.T. 
h 

h. 

“J 

h. 

JOf 

I 


The partial derivative matrix of computed observation 
vector with respect to state vector X 

The partial derivative vector of computed observation 
y^js witli respect to state vector X ; the jth column 
of gT 

The geocentric inertial frame (usually the true-of-date 
frame); Xqj directed toward the vernal equinox of epoch, 
Zgi directed toward the celestial north pole of epoch 

The geocentric inertial frame used internally in OABIAS 

The partial derivative of computed observation y^j with 
respect to x^ ; the 77 th element of G 

The partial derivative matrix of observation vector Y 
with respect to measurement vector M 

The horizon detector line-of-sight frame for in-crossiugs 
z along L_ (reference Figures 3 -1 and 3-5) 

ill 1 

The horizon detector alignment frame for in-crossings 
(reference Figure 3-1) 

The horizon detector line-of-sight frame tor out- 

crossings; z along L (reference Figure 3-1) 

HU O 

The horizon detector alignment frame for out-crossings 
(reference Figure 3-1) 

Higher order terms 

The effective height of the atmosphere (reference Fig- 
ure 3-14) 

T 

The jth column of 1! 

’ r 

The at!i clement of h. 

Identity matrix 
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Symbol 

/N /V 

i, 3. k 

K, K. 

3 

K. 

~3 

L 

L 

/N 

L 

I 

^IGI’ ^2GI’ ^3GI 

^ISC’ *2SC’ ^3SC 
M 


M 

M. 

“3 

m 


XT. 

a 


m. 

30r 





N 


HO 



Definition 


The unit vectors along the axes of frame SC 


The gain matrix of an estimation algorithm 


The gain vector of the estimation algorithm used in 
OABIAS 


The vector from the spacecraft to the horizon crossing 
point (reference Figure 3-14) 

The magnitude of L 

The unit vector directed aloi^ L 

The least-squares loss function freference Equa- 
tion (B-la)] 

The components of L along the axes of frame GI 

The components of L along the axes of frame SC 

A matrix whose three rows are the components of the 
three vectors S , ^ , and 6 , respectively, along the 
axes of frame GI 


A vector composed of a set of scalar measurements 

A* 

A subvector of M 

A vector composcc^ of the four n.tasuremcnts , 

’ SlI ’ ^I!0 ^ single telemetry frame 

An element of M 


An element of 

The unit vectors along the y axes of frames HI and 
HO (reference Figure 3-9) 


The unit vector along the y axis of frame SS (refer- 
ence Figure 3—1) 
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Symbol 


Definition 


”lIlGr '''ll’Gl’ ^II‘'GI components of along the axes of frame GI 

’^HISC* "iPSC’ "lI3SC components of along tlie axes of frame SC 


"l* ”2* "3 


The components of N rUong the axes of frame GI 

s 

A matrix generated in the estimation computations; usu- 
ally considered to be the estimated covariance matrix 
of the error in 5v 


o 

Pr P2 
Q 


^ 2 * ^3 

V^2 


R, |Rl 


It 


R 


/s 

R 


The a priori P matrix 

Defined in Equations (3-75c and 3-75d) 

The inverse P ^ of the estimated covariance 
matrix P of a state vector estimate X 

The columns of the inverse M ^ of the 3x3 matrix M 

Defined in Equations (3~75e) and (3-75f) 

The vector from the spacecraft to the center of the cen- 
tral body 

The magnitude of R 

The unit vector oriented in the direction of R 

Radius of the Earth (reference Figure 3-14) 

The vector from the center of the central body to the 
effective horizon crossing point (reference Figure 3-14) 

The unit vector in the direction of R (reference Fig- 
ure 3-14) ® 


R 


M 


The covariance matrix of a set of measurements M 
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"y 


^ 2 * ^3 


‘'1 

SC 


SI 

SS, ss 


8 


8 . S . S 

1* 2* 3 


T (X.) 
a I 


TOP 
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Definition 

The covariance matrix of a set of observations Y 

The jUi element on the diagonal of R. , 

M 

The components of U along the axes of frame GI 

The components of U along the axes of frame GI 

The unit vector along the spacecraft spin axis 

The weighting matrix in least-squares loss function 
[reference Equation (B-la)] 

The weighting matrix in least-squares loss function 
[reference Equation (B-IO)] 


The spacecraft-fixed reference frame; Zg^ along S 
(reference Figux'e 3-1) 

The inertial frame with Zgj along § (reference Fig- 
ure 3-1) 

The Sun sensor-fixed reference frame (reference Fig- 
ure 3-1) 

Sine 

The components of § along the axes of frame GI; 

8, and s are elements x and x of X 

A 3 X 3 matrix which transforms vecter components 
fro -n the axes of an initial frame, a , to those of a new 
frame, b , where frame b is generated by rotating 
frame a about one of its axes (x. = x , y , or z^) 
through an angle a 

The numerators in the equations for tan , tan A^^ , 


Time 
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Definition . 

The time at which a horizon crossing, in or out, is 
encountered 

The time at which a horizon-in crossing is encountered 

The time at which a horizon-out crossing is encountered 

The time at which an observation j was obtained 

The initial or reference time 

The time at which a Sun sighting is encounter d 

The unit vector along the Sunline (reference Figure 3-4) 

The components of U along the axes of frame GI 

The components of U along the axes oi frame GI 

The velocity vector of the spacecraft center of mass 

The error in observation j 

The observation weighting factor matrix 

The weighting factor for observation j v 

The optin.um value of [reference Equation (3-11)] 

The true state vector 

Thf kimatc of X 

The optimal estimate of X 

The reference value of X 

The subvcclor of X composed of the Sun sensor bias 
parameters 

i* 

The a priori estimate of X 
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Definition 


A vector of "I'eal" observations computed algebraically 
from a vector M ‘of measurements 

A vector of "computed" observations obtained using a 
mathematical model and an estimate of the system's 
state vector X 

The ith element of Y 

The ith element of Y ' 

* -c 

The vector of residuals; Z = Y - Y 

’ - - -c 

The jth element of Z 

✓S 

The right ascension angle of S relative to frame GI 
The right ascension angle of S relative to frame GI 

/s A 

The angle between S and U (refei*ence Figure 3-4) 

The angle measured by the Sun sensor; ideally identical 
to P (reference Figure 3-4) 

The angle between L and S (reference Figure 2-3) 

The angle of horizon detector li '-of-sight vector L 
relative to the detector's reference axis “z ; ideally 
identical to y (reference Figure 2-3) 

The perturbation in dihedral angle A^j due to horizon 
detector misalignment angle e (reference Figure 3-10) 

The pci'turbation in dihedral angle A^j due to Sun sen- 
sor misalignment angle c (reference Figure 3-10) 

A time bins in the location of the spacecraft in its orbit; 
clcmcMt of X 

The correction to slate vector estimate Xj_j or to 
refcreiKC .state vector generated by processing 
observation j [rcfex’unco Ko’iations (3-3) or (3-12)j 
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A)5 


The bias in Sun angle measurement; clement of X 
(reference Figure 3-4) 


Ay 


Ap 


A0 


H 


A0 


O 

H 


6 

6» 




c 

"l 

c 

n 

A 

X 


The horizon detector alignment error; clement of 
X (reference Figure 2-3) 

A bias in th*, effective angular radius p of the central 
body; element x^ of X 

The bias on horizon detector azimuth angle (t>^ at 
horizon-in crossings 

The bias on horizon detector azimuth angle 0^ at 
horizon-out crossings 

/\ 

The declination angle of S relative to frame GI 

/\ 

The declination angle of S relative to frame GV 
The a priori value of 6 ' 

The alignment error of the Sun sensor; element x 
of X (reference Figure 3-4) 

The alignment error of the horizon detector; element 
of X (reference Figure 2-3) 

The angle between horizon crossing vector D and hori- 
zon detector reference axis N„ (reference Figure 3-9) 

The Sun sensor FOV value (reference Figure 2-1) 

90-c (reference Figure 3-10) 

The latitude of the horizon crossing location on the cen- 
tral body surface 


X 

o 


The geocentric latitude of the spacecraft; initial esti- 
mate of X 


V 


The angle between S-D plane and D-11 plane (refer- 
ence Figure 3-13) 
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Snnbol 


Definition 





0 


^IINOM 


The elevation mounting angle of the Sun sensor relative 
to the plane perpendicular to the spacecraft's spin xxis 
(rcfci'cnce Figure 2-1) 

The elevation angle of the Sun relative to the Sun sensor 
(reference Figure 2-1) • 

The angular radius of the central body at the horizon 
crossing point as seen from the spacecraft (reference 
Figure 3-8) 

The computed value of p 

/S ^ 

- R • S cecant € 

li 

The variance of the error in the horizon crossing time 
measurements 

The variance of the error in measurement m. 

ja 

The variance of the error in the Sun sighting time meas- 
urements 

Tlie variance of the error in an observation yj 

The variance of the error in the Sun angle 6,, meas- 
. M 

urements 

Tlie variance of the error in the measurements of the 
panoramic horizon scanner angle y 

s 

/N 

The instantaneous angle between U and the Sun sensor 
axis Ygg (reference Figure 3-4) 

The instantaneous rotation angle of the spacecraft in its 
spin cycle (reference Figure 3-3) 

The azimuth alignment ang! between t .c Sun sensor and 
the horizon detector 

The nominal value of 
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Svml)ol 


0 


H 



X 


The effective value of 0^^ at horizon-in crossings; cle- 
ment X of X (rcfcx'cnce Figure 3-5) 

O 

The effective value of 0^^ at horizon-out crossings; 

element x. of X 
6 — 

U*S secant € 


0 

^0 

CO 


The rotation angle of the spacecraft in its spin cycle at 

reference time t ; element x„ of X 
o 3 - 

The spacecraft's angular rate about its center of mass; 

elenisnt of X 
9 


Superscript 


Definition 


GI, GI’, etc. 


Signifies components of vectors along the axes of coor- 
dinate frames GI, GI', etc. 


I 


Signifies a parameter associated with a horizon-in 
crossing of the horizon detector 


0 

T 


Signifies a parameter associated with a horizon-out 
crossing of the horizon detector 

Signifies the transpose of a matrix or of a colmnn vector 


X 


Signifies a parameter associated with the Xth iteration 
in r 1 iterative processing operation 


Signifies an optimal quantity or estimate 


Subscript 


Definition 


C, c 


Signifies computed 


f Signifies final 

GI, GI*, ct Used with x , y , z , or 1,2,3 to signify the axes 

of the indicated coordinate frain j 
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Subscript 


I, HI 
H 


j 

M, m 
O, HO 


0 
R 

1 


Definition 


Signifies a parameter associated with a horizon-in 
crossing 

Signifies a parameter associated with a horizon cross- 
ing; usually used with symbols that arc applicable to 
both horizon-in and horizon-out crossings 


Usually signifies a parameter associated with the proc- 
essing of observation data set Y or scalar observation 

Signifies measured 


Signifies a parameter associated with a horizon-out 
crossing 


Signifies an a priori or initial value 


Signifies a parameter computed using X = X 

— — K 

Signifies a parameter associated with data set 1 in 
Appendix B 


2 Signifies a parameter associated v'ith data set 2 in 

AppendL'; B 

Signifies the combination of data sets 1 and 2 and 
parameters associated with this set in Appendix B 


+ 


Signifies a pai’ameter computed using the plur sign in 
Equation (3-78) 


Overhead Symbols 


Definition 


“*• Signifies a Cartesian vector of arbilrai'y length 

Signifies a Cartesian vector of unit length; also, signi- 
fies a computed or estimated variable 

r 
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